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ABSTRACT 

A  technique  for  predicting  steady  and  oscillatory  aerodynamic  loads  on 
general  configurations  lias  been  developed  which  is  based  on  the  Doublet-Lattice 
Method  and  the  method  of  images.  Chord-  and  spanwise  loading  on  lifting 
surfaces  and  longitudinal  body  load  distributions  are  determined.  Configura¬ 
tions  may  be  composed  of  an  assemblage  of  bodies  (elliptic  cross,  sections  and 
a  distribution  of  width  or  radius)  and  lifting  surfaces  (arbitrary  planform 
and  dihedral,  with  or  without  control  surfaces).  Loadings  predicted  by  this 
method  are  required  for  flutter,  gust,  frequency  response  and  static  aero- 
elastic  analyses  and  may  be  used  to  determine  static  and  dynamic  stability 
derivatives.  Volume  I  presents  the  theory  and  calculated  results  while 
Volume  II  presents  the  details  of  the  computer  program  used  to  implement 
the  theory. 
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1.0  INTRODUCTION 


Until  several  years  ago  the  kernel  function  procedure  was  the  best  known 
and  most  widely  used  lifting  surface  theory.  The  classic  report  on  this  pro¬ 
cedure  is  by  Watkins,  Runyan  and  Cunningham.1  There  have  been  many  variations  of 
this  procedure;  however,  it  is  not  our  purpose  to  present  a  survey  of  them. 
Ashley,  Widnall  and  Landahl2  have  already  prepared  an  excellent  survey  of  the 
lifting  surface  theory. 

The  use  of  lattice  methods  in  steady  flow  goes  back  to  Fslkner3  and  even 
further.  Recent  developments  and  improvements  made  by  Rubbert4,  Dulmovits5, 
Hedman6,  Belotserkovskii7,  Giesing9  and  others  have  rejuvenated  and  popularized 
this  method.  The  lattice  method  produces  very  accurate  results  even  though 
the  numerical  technique  is  relatively  simple.  James9  has  analyzed  a  two- 
dimensional  model  of  the  steady  lattice  theory.  He  has  proven  that  when  the 
vortex  is  placed  at  the  1 /4-chord  point  of  each  element  and  the  control  point 
is  placed  at  the  3/4-chord  point  of  each  element  (the  lifting  surface  is 
divided  into  a  number  of  elements),  then  the  following  is  true:  1)  the  lift 
and  moment  are  exactly  correct,  2)  the  Kutta  condition  is  satisfied  auto¬ 
matically  without  the  use  of  loading  functions,  3)  in  the  limit  of  a  large 
number  of  elements  the  correct  leading  edge  singularity  and  correct  trailing 
edge  zero  are  obtained,  and  4)  the  method  works  just  as  wall  for  cambered 
surfaces  whose  distributions  are  continuous. 

The  lattice  method  is  simple,  versatile  and  accurate.  The  accuracy  has 
already  been  discussed.  The  simplicity  arises  from  the  fact  that  no  loading 
functions  are  required.  The  method  is  versatile  because  no  prior  knowledge 
of  the  solution  is  required  as  with  the  kernel  function  procedure.  Loading 
functions  appropriate  to  the  particular  problem  must  be  chosen  ahead  of  time 
for  the  kernel  function  procedure.  It  has  been  demonstrated  that  the  lattice 
method  can  handle  a  very  wide  variety  of  configurations,  including:  1)  wings 
of  arbitrary  planform,  2)  wings  with  partial  span  control  surfaces,  3)  T-tails, 
wing-pylons,  wing-tails,  4)  wings  in  ground  effect,  and  5)  annular  wings,  wings 
with  arbitrary  dihedral. 
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Since  the  steady  lattice  method  has  met  with  such  success,  it  seems 
only  natural  that  it  should  be  extended  to  unsteady  (oscillatory)  flow. 

Albano  and  Rodden10  have  done  just  that.  (A  completely  independent  unsteady 
lattice  procedure  was  developed  by  Stark11.) 

The  new  unsteady  method  has  been  termed  the  Doublet-Lattice  Method  (DLM). 
Further  extensions,  applications  and  refinements  of  the  DLM  are  found  in 
references  12  through  18.  The  DLM  has  proven  as  versatile  as  the  Vortex 
Lattice  Method  and  can  handle  the  same  wide  variety  of  configurations  (see 
references  13,  18  and  Part  I  of  this  report).  The  attributes  of  the  Vortex 
Lattice  Method,  i.e.,  simplicity,  accuracy  and  versatility  can  be  applied 
equally  well  to  the  DLM.  Because  of  these  attributes  it  was  selected  as 
the  basis  of  the  present  method. 

Of  major  interest  in  this  report  is  the  interaction  of  bodies  (e.g., 
fuselage,  nacelle,  store)  with  lifting  surfaces  (e.g.,  wing,  tail,  pylon, etc.). 
There  are  two  basic  approaches  to  this  problem:  1)  one  in  which  elements  are 
placed  on  the  body  surfaces,  and  2)  one  in  which  images  are  placed  within  the 
bodies.  Recent  advances  in  the  first  of  these  approaches  have  been  made  by 
Woodward19,  Labrujere20,  and  Bradley  and  Miller21  for  steady  flow.  Part  I  of 
this  report  presents  an  extension  of  these  methods  of  oscillatory  flow.  Recent 
advances  in  the  second  of  these  approaches  have  been  made  by  Giesing8,  Spangler 
and  Mendenhall22,  3orland23,  and  Chou24  for  steady  Flow.  This  report  (Part  II) 
presents  an  extension  of  these  methods  to  oscillatory  flow  and  very  general 
configurations. 

The  approach  of  this  report  is  to  use  the  method  of  images  directly 
replacing  the  steady  vortex  lattice  on  the  lifting  surfaces  and  image  surfaces 
with  an  unsteady  vortex  lattice  DLM.  The  advantage  of  using  images  over  other 
methods  is  the  fact  that  images  do  not  introduce  any  new  unknowns  into  the 
problem. 

This  approach  furnishes  a  practical  method  for  handling  general  configu¬ 
rations  efficiently.  Specifically,  the  configurations  considered  may  include 
a  collection  of  bodies  (e.g.,  fuselage,  nacelles,  stores)  and  lifting  surfaces 
(e.g.,  wing,  tail,  pylon,  etc.).  The  configuration  may  oscillate  in  any  mode, 
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rigid  or  flexible,  and  may  operate  in  or  out  of  ground  effect.  Outputs  from 
the  method  are:  generalized  forces  and  aerodynamic  parameters  such  as  span 
loads,  center  of  pressure,  lift,  moment,  etc.  The  method  can  also  be  used  to 
obtain  dynamic  stability  derivatives  (reference  14). 
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2.0  THEORETICAL  DEVELOPMENT 


2.1  Outline  of  the  Problem  and  Solution 

The  velocity  normal  to  an  oscillating  surface  or  body,  W  =  UooRe(wse1  *) 
is  related  to  the  lifting  pressure  Ap  =  qReUCpe1^)  by  the  following  inte¬ 
gral  equation  derived  in  Appendix  A: 

ws(x,y,z,Yr)  =  §-  fj  K(x  -  y  -  n,  z  -  y$.  Yr»  M)ACp(f  ,n,c)ds 

L.S.  (2.1-1) 

where  L.S.  indicates  integration  over  all  lifting  surfaces.  This  is  the 
familiar  integral  equation  of  lifting  surface  theory  for  surfaces  alone. 

Here  y$>  Yr  are  the  sending  and  receiving  surface  dihedral  anqles.  If  a 
body  is  introduced,  there  are  additional  contributions  to  the  normal  velocity. 
The  first  contribution  may  be  called  the  slender  body  term  and  represents  the 
flow  field  generated  by  bodies  without  considerations  of  interference. 

AW(xry,z,Yr)  =  §7  f  L(x  -  c,  y  -  na,  z  -  v$,  Yr*  M)^(?)dc 

P. 

+  M(x  -  C,  y  -  na,  z  ~  r,a>  Yr»  M)aU)d?  (2.1-2) 

B. 

The  limit  B.  indicates  integration  over  all  bodies;  these  integrals  are  also 
derived  in  Appendix  A.  The  subcript  a  on  n  and  x,  indicate  the  location 
of  body  ayes. 

Here  the  term  us(c)*  represents  an  axial  multipole  distribution  (dipole, 
quadrupole,  etc.)  whose  orientation  is  given  by  y$-  The  second  integral 
exists  only  in  steady  flow  and  is  a  source  distribution  used  to  represent 
the  body  volume  effects.  The  sler.der  body  terms  are  known  since  ns(?)  ancl 
r(?)** are  determined  using  an  appropriate  slender  body  theory. 


*Actual ly  pq (s )  is  twice  the  classic  doublet  strength  uq  Co¬ 
ncurrently  in  the  Present  Method  sources  are  excluded  since  they  do  not 
exist  for  unsteady  flow  and  have  a  '.mail  effect  even  in  steady  flow. 
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A  second  contribution  to  the  normalwash  flow  field  caused  by  the  intro¬ 
duction  of  bodies  into  the  flow  field  arises  from  "image"  lifting  surface 
elements.  These  are  placed  within  the  body  to  help  divert  the  flow  around 
a  body  in  the  presence  of  a  lifting  surface.  The  strength  of  the  image 
elements  is  the  same  as  the  external  elements. 


'r(x>y>z»Yp)  =  J'J'  £>  y  Oj»  2  £j»  Yg i »  Yp>  k^,  M)aC  (?>n»?)ds 


L.S. 


(2.1-3) 


The  subscript!  on  ni?  and  y$  indicates  the  image  position  on  the  image 
surface. 

A  third  contribution  to  the  normalwash  How  field  caused  by  the  bodies 
is  generated  by  an  interference  doublet  distribution  y„(s). 


wn(x,y,z,Yr)  =  h  J  L(x  -  5,  y  ~  na,  2  -  ca,  y$»  Yr>  kr>  M)vTn(c)dc  (2.1-4) 
B. 

The  form  of  this  equation  is  identical  to  the  first  integral  of  equation 

(2.1-2).  The  reason  it  is  not  combined  with  us  is  the  fact  that  us  is 

known  while  £  ,  like  AC  .  is  unknown, 
n  p 

In  the  direct  problem,  the  normalwash  boundary  conditions  are  specified 
and  the  lifting  pressure,  AC  ,  and  body  interference  distribution,  pj,  are 
solved  for: 


w  =  w  +  aw  +  wT  +  vir 


(2.1-5) 


Here  w  is  the  prescribed  normalwash  on  lifting  surfaces  and  bodies.  If  the 
known  quantities  are  placed  on  the  left-hand  side  while  the  unknown  quantities 
are  placed  on  the  right,  then  the  following  equation  results: 


w  - 


(L 

8tt 


/  LIrsd?  +  h  f  Madc)=  bt  //  KACPd5  +  h  If  KiVs 

B.  B.  /  L.S.  L.S. 


8tt 


J 


d5 


(2.1-6) 


B. 
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Here  the  subscript  I  on  K  indicates  K(x  y-nI#  z  YsI,  yr»  K>  M)- 

"  - 4W  ■  h  If  <K  +  Ki>Vs  +  ar /  L«nd«  <2J-7) 

L.S.  B. 

There  are  still  further  contributions  to  the  normalwash  and  these  arise 
from  planes  of  symmetry  and  ground  effect.  If  the  assumption  is  made  that  the 
right  half  of  the  aircraft  Tier  in  the  upper  right-hand  quadrant  of  the  z-y 
plane,  then  a  subscript  UR  may  be  applied  to  contributions  made  from  lifting 
surfaces  in  this  quadrant.  Similarly,  UL  indicates  upper  left  which  contains 
the  contribution  from  the  left  side  of  the  aircraft.  The  subscript  LR 
indicates  lower  right  and, in  tcis  quadrant,  the  ground  effectof  the  right  side 
is  contained.  The  subscript  '„L  indicates  lower  left  and  this  quadrant  * 
contains  the  contribution  of  the  ground  effect  of  the  left  half  of  the  air¬ 
craft.  F.quation  (2.1-7)  may  be  expanded  to  include  these  contributions  as 
follows: 

W  —  (AWjjp  +  £AWy^  +  eMyp  +  e<$AWp> 

=  +  KI)UR  +  +  KI>UL  "  e^K  +  KI}LR  +  ei^K  +  KI>Ll{  ACpds 

+  87  f  J  Lur  +  6LUL  +  eLLR  +  e<5LU.  ^nd?  (2.1-8) 

B. 

The  quantities  6  and  c  are  the  symmetry  and  ground  effect  indicators. 

For  instance,  a  =  1 ,0,-1  indicates  s>mnnetry,  no  symmetry,  and  antisymmetry, 
respectively.  Similarly,  e  =  -1,0,1  ndicates  ground  effect,  no  ground 
effect,  and  anti-ground  effect,  respe:ti„'o  .  The  changes  to  the  argument  lists 
denoted  by  these  subscripts  are  as  follow:, 

UR:  n  =  n,  C  =  ?,  ys  =  y$ 

UL:  n  =  -n>  s  =  Ys  =  "Ys 

LR--  n  =  n,  t,  =  ~C,  Ys  =  _YS 

LL  :  rt  =  -r'j  C  =  -c,  Ys  =  Y$ 


6 


The  bade  method  of  solution  of  Eq.  (2.1-8)  is  to  discretize  the  lifting 
surfaces  into  small  boxes  and  the  bodies  into  sma  axial  elements.  The 
unknowns  are  assumed  constant  over  these  elements  and  the  normal wash  boundary 
condition  is  applied  to  each  box  and  element.  This  forms  ->  many  equations  as 
unknowns  and  the  system  may  be  solved.  Eq.  (2.1-8)  becomes: 


N1  N2 

"  S  ACp  ff  KT  ds  +  S  un  /  LT 
s=l  S  ELEMENTS  rs  s=l  S  ELEMENTg  rs 


(2.1-9) 


where  s  and  r  indicate  sending  and  receiving  points,  respectively,  and 


Wy  =  w  —  aWy 

AWy  =  (awur  +  6AWul  +  sAV.'lR  +  e6AWLL) 

Ky  =  |  (K  +  ^j)yR  +  5(K  +  Ky  )yj_  +  £(^  +  )^R  +  e<5(^  +  Kj )yy  j 

4  =  81T  { LUR  +  6LUL  +  eLLR  +  £5LLl} 

N1  =  number  of  lifting  surface  boxes  for  all  surfaces 
N2  =  number  of  axial  body  elements  for  al i  bodies 

In  matrix  notation: 

/  i  (aC  ) 

jwyj  =  [[)_]  (2.1-10) 

where 

[Dy]  =  [D  \E]  (2.1-11) 

in  which 

D=  ((D+D,)UR 

+6(0+  0j)yy 

4  e(D  +  °I)LR 

+  e6(D  +  Dt)  s  (2.1-12) 

1  LLj- 
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Here 


°rs  ■  //  ^  K  ds 
ELEMENT 

DI  =  ffh  Ki  ds  =  IJ  ET  K<x  ~  C.  y  -  nr  z  -  Cj,  yr»  Ys  .  kr,M)ds 


ELEMENT 


ELEMENT 


(2.1-13) 


The  matrix  elements  D^s  have  the  subscript  r  on  the  receiving  point  quanti¬ 
ties,  x,y,z,Yr  and  the  subscript  s  on  the  sending  element  quantities 
C,n,?,Ys>  ELEMENT.  The  matrix  partition  E  is: 


Ers  =  471  EUR  +  6EUL  +  eELR  +  f6ELL; 


(2.1-14) 


BODY 

ELEMENT 

Once  wT  is  known,  AC^  and  can  be  found,  and  these  can  be  used  to  find 
the  loads  on  the  lifting  surfaces  ard  bodies. 

The  calculation  of  the  flow  fic  d  due  to  the  slender  body  terms,  i.e., 
us  and  a t  is  performed  using  the  same  discretization  technique. 

N3  N3 

-hi  “s  f  /  MrsdE  f2-1-15* 


BODY 

ELEMENT. 


k=l  BODY 
ELEMENT. 


where  N3  is  the  number  of  slender  body  elements.  If  symmetry  planes  and 
ground  effect  are  accounted  for  and  matrix  notation  is  introduced,  then  Eq. 
(2.1-15)  becomes: 

{AwtJ  =  [Lt]  {ys)  +  [MT]{a}  (2.1-16) 


where 


MT  =  47  {MUR  +  6MUL  +  eMLR  +  e6ML’l} 
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Slender  body  theory  states  that  *y  is  directly  proportional  to  the 
local  velocity  normal  to  the  body  axis  (the  direction  of  us  is  parallel  to 
this  velocity). 


vs  =  wD2D  (2.1-17) 

where  D2D  is  the  proportionality  constant  which  is  dependent  on  local  Dody 
cross  section.  The  values  of  w  which  act  normal  to  the  body  axis 

(Eq. (2.1-17))  are  part  of  the  larger  set  that  acts  normal  to  all  surfaces 
and  bodies. 

2.2  The  Normalwash  Boundary  Conditions 

The  normalwash  w  must  be  determined  at  each  lifting  surface  element 
(or  box)  and  at  each  axial  body  element  in  both  the  z-  and  y-di recti  one. 

The  normalwash  boundary  conditions  are  obtained  by  taking  the  substantial 
derivative  of  the  modal  deflections.  There  are  various  methods  of  describing 
these  modes  and  several  of  these  will  be  discussed.  Only  the  first  of  these 
methods  has  been  incorporated  into  the  present  method  (the  polynomial  approach). 
The  polynomial  approach  lends  itself  to  scientific  investigation  where  the 
modes  are  simple.  When  the  modes  become  complicated,  however,  it  may  be 
desirable  to  incorporate  other  more  practical  modal  input  methods. 

2.2.1  Polynomial  Modes 

The  total  deflection  distribution  of  a  lifting  surface  normal  to  itself 
is  made  up  of  a  set  of  modes,  fj. 

NM 

"■'I  Vi  (2-z-1) 

ie 


where  are  the  generalized  coordinates  and  NM  is  the  number  of  modes. 
The  total  normalwash  is  likewise 


co 


I  Vi 

i=l 


(2.2-2) 
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where 


or 


Here 


(  df.  -  v 

w,  =  -  I — —  +  i  f-  f. 


d (x/c) 


dfi 

w,  =  -< - —  +  i  k  2f . 

!  I  j  /  /T\  r  i 


,  d  ( x/'c) 


(oC 

r  "  2U 


For  lifting  surfaces  the  modes  may  be  approximated  by 


(2.2-3) 


(2.2-4) 


where  -  is  the  lateral  distance  in  the  plane  of  the  lifting  surface. 
Eqs .  (2.2-3)  and  (2.2-4)  represent  fifth-degree  polynomials  in  both  the 
lateral  and  longitudinal  directions. 


When  dealing  with  bodies,  two  separate  directions  of  motion  are  possible: 
z-motion  and  y-motion.  However,  bodies  have  no  lateral  coordinate,  thus  the 
mode  shapes  are  as  follows: 


(2.2-5) 


(2.2-6) 
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(2.2-7) 


(2.2-8) 


The  w.  arrays  have  the  following  order:  first,  all  of  the  lifting  surface 
normal  wash  values,  w,  are  determined,  then  the  z-  or  upwash  values,  w^,  and 
then  the  y-  or  sidewash  values,  w  . 


w 


(2.2-9) 


2.2.2  Numerical  Input 

In  many  instances  it  is  inconvenient  to  determine  polynomial  coefficients 
from  modal  data.  The  values  of  w  and  f  could  then  be  supplied  directly. 

An  alternate  scheme  is  to  supply  only  h  and  require  the  program  to  take 
the  necessary  derivatives  numerically.  Various  fitting  techniques  could  be 
used  including  the  spline  fit  of  Harder25. 

A  second  alternate  is  to  supply  only  the  values  of  df../d(x/F)  and 
integrate  for  the  values  of  f. .  Aiso  needed  is  one  value  of  f.  at  eaui 
spanwise  station  to  establish  the  level  of  f^. 

Interpolation  schemes  could  be  used  to  reduce  the  number  of  input  values. 
Instead  of  supplying  deflection  data  at  each  spanw'se  strip,  it  could  be 
supplied  at  intervals  along  the  span  and  interpolated  at  the  intermediate 
spanwise  locations. 
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2.2.3  Aerodynamic  Influence  Coefficients 


There  are  many  schemes  for  generating  aerodynamic  influence  coefficients 
( AlCfe )  and  some  of  these  are  discussed  in  Part  I  of  this  report  (page  28). 

The  basic  idea  of  the  AIC  approach  is  to  derive  a  set  of  elementary  modes 
which  can  be  used  to  build  up  any  desired  mode  shape.  These  elementary  mode 
shapes  are  independent  of  the  mass  or  stiffness  properties  of  the  aircraft 
and,  thus,  can  be  used  to  build  up  a  solution  for  various  fuel  conditions, 
for  example.  One  way  to  derive  elementary  mode  shapes  is  to  deflect  a  series 
of  modal  deflection  points  on  the  structure  one  at  a  time  while  holding  all 
the  others  fixed.  Once  the  generalized  forces  for  this  set  of  elementary 
modes  are  known  for  a  particular  aircraft  planform,  Mach  number  and  frequency, 
aerodynamic  solutions  for  all  mass  and  stiffness  distributions  can  be  found 
simply  by  a  matrix  multiplication  of  the  modal  deflections  and  the  AIC 
matrix. 


Qi(j  =  (fj)1  [AIC]  {f.}  (2.2-10) 

where  f .  is  a  set  of  modal  deflections  for  mode  j  and  f.  is  the  same 

J  * 

for  inode  i . 

The  type  of  curve  fit  used  for  the  generation  of  the  elementary  mode 
shapes  may  vary.  Part  I  of  this  report  offers  several  possibilities  (see 
page  31).  Recent  work  by  Harder,  et.al.25  already  referred  to,  on  two- 
dimensional  surface  spline  fitting  of  deflected  surfaces  may  offer  a  very 
accurate  and  versatile  technique  for  generating  elementary  mode  shapes. 

For  bodies  the  usual  one-dimensional  spline  could  be  used. 

2.3  The  Lifting  Surface 

The  basic  technique  employed  to  represent  the  lifting  surfaces  is  the 
Doublet  Lattice  Method  (DLM).  The  DLM  is  an  extension  of  the  steady  flow 
Vortex-Lattice  Method  to  oscillatory  flow.  Developments  of  the  Vortex-Lattice 
Method  include  those  of  Rubber*4,  Dulmovits5,  Hedman6,  Belotserkovskii7,and 
Giesing8  who  extended  the  method  to  wing-body  cc  ibinations.  The  original 
extension  of  the  steady  Vortex-Lattice  Method  to  unsteady  flow  was  made  by 
Albano  and  Rodden10.  Adaptation  of  the  DLM  to  AIC  generation  was  done  by 
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Stahl,  et  al.12  Extensive  correlations  for  nonplanar  configurations  were 
made  by  Kalman,  et  al.13  A  refinement  of  the  method  for  nearly  coplanar 
configurations  is  reported  in  Part  I  of  this  report  (section  2.1,  appendices 
A,  B  and  C)  and  also  by  Rodden,  et  al."8  The  use  of  the  method  for  stability 
derivatives  is  given  by  Rodden  and  Giesing15,  and  its  use  in  gust  analysis  is 
given  by  Giesing  and  Rodden16.  An  extension  to  induced  drag  distribution  was 
made  by  Kalman,  et  al.17 

The  DLM  has  proven  to  be  a  simple,  versatile  and  accurate  method  for 
the  solution  of  unsteady  nonplanar  lifting  surface  problems.  The  general 
nature  of  the  solution  need  not  be  inown  in  advance  as  with  the  kernel  func¬ 
tion  techniques  which  utilize  a  standard  set  of  pressure  loading  functions. 
Elimination  of  the  loading  functions  simplifies  the  analysis  since  the 
complicated  and  time  consuming  quadrature  integrations  of  the  loading  functions 
and  kernel  are  eliminated.  The  method  is  versatile  since  there  are  essentially 
no  restrictions  on  the  configurations  that  can  be  handled.  It  has  been  shown 
to  be  accurate  by  James9  who  performed  analytic  studies  on  the  steady  two- 
dimensional  version, and  through  the  many  correlations  with  experimental  data 
and  other  analytical  methods. 

The  flow  singularities  used  to  model  the  lifting  surface  are  steady 
horseshoe  vortices  and  oscillatory  doublets  along  the  bound  vortex.  The  doublet 
line  is  equivalent, at  zero  frequency,  to  the  horseshoe  vortex,  and  thus  the 
horseshoe  vortex  need  not  be  used.  However,  since  the  effects  of  the  vortex 
system  can  be  analyzed  exactly,  while  the  effects  of  a  doublet  line  can  only 
be  approximated,  improved  accuracy  is  obtained  by  using  both  the  vortex  and 
doublet  systems.  In  this  way,  the  vortices  represent  the  steady-flow  effects, 
and  the  doublets  represent  the  incremental  effects  of  oscillatory  motion. 

The  configuration  is  idealized  by  dividing  the  surface(s)  into  small 
trapezoidal  elements  (ooxes)  arranged  in  strips  parallel  to  the  freestream 
so  that  surface  edges,  fold  lines,  and  hinge  lines  lie  on  box  boundaries 
(sketch  2.3-1).  Then,  to  '-epresent  the  steady-flow  effects,  a  horseshoe 
vortex  is  placed  on  each  of  the  boxes  such  that  the  bound  vortex  of  the  horse¬ 
shoe  system  coincides  with  the  quarter-chord  line  of  the  box.  To  represent 
the  oscillatory  increment,  a  distribution  of  acceleration  potential  doublets 
(which  have  the  steady-flow  ocublet  strength  subtracted)  of  uniform  strength 
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Sketch  2.3-1.  Surface  Idealization  into  Boxes  and  Location  of  Vortices, 
Doublets  and  Collocation  Points 

is  superimposed  on  the  bound  vortex.  The  surface  boundary  condition  is  a 
prescribed  normalwash  applied  at  the  control  point  of  each  box.  The  control 
point  is  centered  spanwise  on  the  three-quarter-chord  line  of  the  box 
(sketch  2.3-1).  The  influences  of  all  vortices  and  doublets  are  summed  for 
each  control  point  to  obtain  the  total  dimensionless  normalwash,  w,  at  the 
control  point. 

The  expressions  for  Drs,  the  elements  of  the  influence  matrix,  are 
well  documented  in  Part  I  and  will  not  be  repeated  here.  A  general  description 
of  the  method  of  integration  will  suffice.  The  expression  for  Dr$  is  given 
in  Eq.  (2.1-13). 


D  =  ff  Kds 
rs  8^  JJ 

ELEMENT 


(2.3-1) 
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The  kernel  K 


is  split  up  into  the  following  components 
,(s) 


K  = 


Vi 


T2K2S) 


(2.3-2) 


where  the  sending  and  receiving  subscripts  (s,r)  have  been  dropped  for 
convenience.  The  first  term  in  braces  is  the  steady  part  of  the  kernel  while 
the  second  and  third  terms  are  the  planar  and  nonplanar  parts,  respectively, 
of  the  unsteady  increment.  Integration  in  the  ^-direction  (parallel  to  the 
freestream)  is  done  by  lumping  the  value  of  the  integrand  at  the  one-quarter- 
chord  point  of  the  box.  The  integration  of  the  steady  term  in  the  spanwise 
direction  is  done  using  the  law  of  Biot  and  Savart.  Spanwise  integration  of 
the  planar  and  nonplanar  terms  is  done  by  approximating  the  numerators  by 
parabolas.  The  resulting  expressions  are  analytically  integrable. 


D  = 

f  (VlS)  ,  V2S)|d- 

*-7H* 

/ 1 '  a.,  r  Ain2  +  Bin  +  Cn 

n  V  /  =  I  _! _ ! _ L 


8 


-e 

e 


.(2)  _  ax  f  kf  '  +  B2n  +  C2 


8tt 


/ 

-e 


dn 


where 


A^2  +  Bf,  +  C]  «  T-j (K]  -  Kjs)) 


A^i  2  +  B2n  +  C2  «  T2(K2  - 


(s)> 


(2.3-3) 


(2.3-4) 


(2.3-5) 


(2.3-6) 


(2.3-7) 
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the  coefficients  A,B,  and  C  are  found  by  evaluating  the  kernel  at  the  center 
and  both  edges  of  the  element.  Here  n~  is  a  spanwise  coordinate  in  the  plane 
of  the  element  centered  on  the  sending  element,  e  is  the  semi -width  of  the 
sending  element  and  the  subscripts  (r,s)  have  been  dropped  again  for  con¬ 
venience.  The  result  of  inteyr  ting  Eq.  (2.3-4),  i.e.,  an  expression  for  the 
steady  horseshoe  vortex,  is  found  in  Par:  I,  Appendix  C,  Vol .  I,  page  68.  The 
integrations  indicated  in  Eqs.  (2.3-5)  and  (2.3-6)  are  found  also  in  Part  I, 

Vol.  I,  Appendix  B,  page  58.  Two  expressions  for  are  found  on  pages  59 

and  60  of  Part  T,  Vol.  1,  Appendix  B. 

The  matrix  formulation  for  an  isolated  lifting  surface  without  ground  effect, 
symmetry  or  images,  given  in  terms  of  the  influence  coefficient  matrix  [D],  is 
then 

{w}  =  [D]  {ACp}  (2.3-8) 


2.4  The  Isolated  Body 

It  is  desirable  to  use  an  axial  system  of  singularities  to  simulate  the 
effect  of  bodies  in  the  fluid  since  they  require  much  less  computational 
effort  than  do  surface  singularity  distributions.  Generally, axial  distribu¬ 
tions  are  associated  with  slender  body  theory. 

Slender  body  theory  has  had  a  long  history  and  it  is  not  the  purpose  of 
this  report  to  detail  its  development.  Only  some  of  the  more  important  works 
will  be  mentioned.  The  doctoral  thesis  by  Revel  1  (see  reference  26)  gives  a 
very  elegant  and  detailed  description  and  cites  over  200  references. 

Munk27  is  generally  the  first  of  any  references  relating  to  steady 
slender  body  theory.  Kantian20  applied  a  modified  version  of  Munk's  theory  to 
airship  hulls  (Karmafi's  method  will  be  discussed  later  in  this  section). 
Lighthi 11 29  may  be  mentioned  next  along  with  Jones30  for  extensions  to  slender 
wings  and  Miles31  for  extension  to  unsteady  flow.  Many  attempts  have  been 
made  to  extend  the  slender  body  theory  to  higher  order.  Ward32  has  given 
second-order  effects  for  the  nonlifting  component  of  the  pressure  for  super¬ 
sonic  flow.  Adams  and  Sears33  have  extended  this  method  to  subsonic  flow 
using  the  Fourier  transform.  Van  Dyke34  attempted  a  second-order  "lifting" 
solution  and  was  only  partially  successful.  Van  Dyke  suggests  using  the 
first-order  solution  for  the  lifting  components  (doublets)  but  suggests  using 
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the  higher  order  axial  flow  solution  along  with  a  higher-order  pressure  formu¬ 
lation  to  find  the  pressures  and  forces.  Such  a  formulation  allows  the  coupling 
of  the  cross  flow  with  the  axial  flow.  This  coupling  does  have  an  effect  on 
the  axial  loading;  however,  for  many  bodies  (shown  in  reference  34)  the  first- 
order  solution  is  very  good  over  a  wide  range  of  Mach  number. 

The  only  completely  consistent  second-order  theory  has  been  developed  by 
Revel!26.  Revel!  attacks  the  second-order  solution  by  a  complicated  iteration 
scheme. 

The  methods  of  Hess  and  Smith35  and  Landweber36  may  be  mentioned  in 
passing.  These  methods  are  exact,  surface  singularity  methods.  They  are, 
however,  only  valid  for  steady,  compressible  flow  or  unsteady  incompressible 
flow  and  cannot  be  used  here. 

In  all  of  the  axial  singularity  methods  described  (except  for  Revell),  the 
lifting  part  of  the  solution  is  given  by  the  first-order  slender  body  theory. 
Various  improvements  have  been  suggested  for  determining  the  loading  associated 
with  this  distribution;  however,  the  distribution  itself  is  determined  using 
first-order  theory. 

The  basic  approach  of  slender  body  theory  is  to  consider  only  the  near 
field  effects  of  the  axial  singularities  of  the  bodies  on  themselves.  In  the 
near  field  the  flow  is  both  two-dimensional  and  q^as , -steady.  When  the  body 
surface  lies  close  to  the  axis,  the  surface  sees  only  the  local  axial  singu¬ 
larity  distribution.  Also,  when  the  body  surface  lies  close  to  the  axis, 
the  characteristic  wave  length  associated  with  the  solution  of  the  wave  equa¬ 
tion  is  long  compared  to  the  distance  from  the  body  axis  to  its  surface;  this 
renders  the  problem  quasi-steady.  In  all  of  the  axial  singularity  methods 
described,  only  one  isolated  S'dy  is  considered.  Any  method  developed  would 
have  to  properly  generate  a  flow  field  away  from  the  body  surface  where  the 
slender  body  assumptions  are  no  longer  valid  (so  that  the  mutual  interference 
of  bodies  can  be  accounted  for). 

A  possible  solution  to  this  problem  would  be  to  use  the  exact  three- 
dimensional  unsteady  compressible  flow  solution  for  the  effect  of  the  axial 
singularity  system.  The  axial  integral  could  be  discretized  (see  Appendix  B) 
giving  a  piecewise  constant  singularity  distribution.  For  such  a  distribution 
the  axis  is  di /i ded  up  into  a  series  of  short  axial  elements.  The  singularity 
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strengths  are  held  constant  over  each  element  but  varied  from  element  to 
element.  The  boundary  conditions  would  be  made  to  hold  a^ng  a  line  on  the 
surface  of  each  of  the  bodies.  Multiple  bodies  of  arbitrary  thickness  should 
be  handled  easily  and  the  flow  fields  could  be  calculated  easily. 

Karmen  used  this  approach  to  predict  pressures  and  loads  on  airship  hulls 
at  angle  of  attack.  He  used  sources  for  the  volume  effect  and  doublets  for 
the  crossflow  lifting  effect.  Karman's  method  was  a  velocity  potential 
approach  (Appendix  B,  Eq.  (B-2)).  Other  approaches  could  be  used,  however, 
for  instance,  the  pressure  potential  method  given  by  Eq.  (B-4)  (Appendix  B) 
or  the  modified  acceleration  potential  approach  given  by  Eq.(B-15)  (Appen¬ 
dix  B).  For  steady  flow  the  modified  acceleration  potential  method  is  the  same 
as  the  velocity  potential  method. 

The  unsteady  pressure  potential  and  modified  acceleration  potential 
methods,  discretized  as  described  above,  were  implemented  on  the  computer  so 
that  they  could  be  investigated.  The  formulas  for  the  flow  field  due  to  a 
short  element  of  pressure  potential  doublet  and  modified  acceleration  potential 
doublet  are  as  fallows: 


Ers  ~  8tt  AfjsK^xr  ?l/4’  yr 


V  Zr 


Ca  1  Yy.>  k  ^  ^ 


Pressure  Potential 


(2.4-1) 


i  wAC./2U 

e  "  K(xr  -  C,  ,  yr  -  n,  . 


u:AC_/2U 

5  "  K(xr 


Modified  Acceleration  Potential  (2.4-2) 


The  terms  q  ,  ^  and  5]/4  are  the  leading  edge,  trailing  edge  and 
1/4-chord  point  of  the  sending  axial  element  (indexed  by  the  subscript  s). 
The  index  r  indicates  the  receiving  point.  The  boundary  condition  was 
enforced  at  one  point  per  element.  For  the  pressure  potential  method  this 
point  was  located  as  follows: 
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Pressure  Potential 


xr  =  x3/4  (3/4-chord  point  of  element) 

yr  =  nd  (2.4-3) 

r  ^a  r 

where  a^  is  the  radius  at  the  point  . 

Modified  Acceleration  Potential 

xr  =  xl/2  (mid-chord  point  of  element) 

yr  =  na  (2.4-4) 

z  =  t  +  a 
r  sa  r 

Unlike  slender  body  theory  each  element  affects  all  control  points. 

The  pressure  point  doublet  in  steady  flow  is  just  an  ordinary  doublet  line 
of  constant  strength  originating  at  the  pressure  point  and  terminating  at 
downstream  infinity.  The  pressure  point  doublet  in  steady  flow  then  is  just 
a  semi -infinite  doublet  line. 


z 


Sketch  2.4-1 
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The  two  approaches  previously  outlined  are  essentially  the  same  for 
bodies  that  close  {eliminating  the  vake).  The  only  difference  is  that  the 
elements  for  the  pressure  doublet  method  are  shifted  aft  by  1/4  of  an  element 
length.  This  difference  disappears  with  increasing  numbers  of  elements.  The 
numerical  results  for  either  method  were  found  to  be  almost  identical  for  the 
cases  considered.  Figure  1  shows  one  such  case.  A  pointed  body  of  parabolic 
shape  is  made  to  pitch  about  its  leading  edge  with  a  reduced  frequency 
(k  -  wL/Uj  of  0.1.  The  calculated  axial  loading*  is  compared  to  that  obtained 
using  the  method  of  J.  W.  Miles31.  For  this  case  the  axis  was  divided  up  into 
ten  equal  elements.  A  pressure  doublet  was  placed  on  each  element.  The  results 
for  a  thickness  ratio,  (Rmax/L)}  of  0.1  are  essentially  in  perfect  agreement 
with  these  of  Miles.  As  the  thickness  increases,  however,  the  calculated  results 
stare  to  oscillate  about  the  correct  solution.  F'"  very  thick  bodies  the 
results  are  mtaningless.  It  has  been  found  that  the  ratio  of  diameter  to 
element  length  plays  an  important  role  in  determining  the  accuracy  of  the 
results.  It  has  been  found  that  the  diameter  should  not  be  much  greater  than 
2  to  4  times  the  element  length.  Von  Karman  used  a  diameter-to-element  length 
ratio  of  about  3,  and  thus  obtained  very  satisfactory  results  for  the-  airship 
hulls  he  considered. 

The  reason  for  the  failure  of  this  method  is  not  known.  The  supersonic 
version  of  this  method,  as  developed  by  Ktfrmctn  and  Moore37  and  Tsien38,  seems 
to  work  very  well.  It  is  suspected  that  the  reason  the  supersonic  method  works 
while  the  subsonic  method  runs  into  trouble,  is  the  fact  that  in  supersonic 
flow,  singularities  or  discontinuities  are  propagated  to  the  body  surface  from 
the  axis  while  they  are  not  in  subsonic  flow.  It  is  obviously  impossible  to 
produce  a  discontinuity  in  upwash  on  the  body  surface  due  to  axial  singularities 
unless  discontinuities  can  propagate  to  the  surface  of  the  body  from  the  axis. 

It  is  then  obvious  that  in  subsonic  flow  the  closer  the  singularities  are  to 
the  surface  ( i . e . ,  the  more  slender  the  body)  the  more  accurate  the  solution 
will  be.  It  seems  logical  that  there  is  a  restricted  family  of  upwash  distri¬ 
butions  or  boundary  conditions  that  a  subsonic  axial  singularity  system  can 
satisfy  on  the  surface  of  the  body.  The  method  by  which  one  determines  whether 
a  specific  upwash  boundary  condition  qualifies  for  a  solution  is  not  yet  known. 

*The  assumption  was  made  that  the  axial  doublet  strength  could  be  used  directly 
to  obtain  the  loading  in  the  calculation. 
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Karman28  commented  that  the  doublet  distribution  obtained  by  his  method 
agrees  very  well  with  the  doublet  distribution  obtained  using  Hunk's  slender 
body  theory.  He  went  on  to  suggest  that  very  accurate  solutions  for  the 
pressure  could  be  found  if  the  exact  non-linearized  pressure  formula  were 
used  in  conjunction  with  the  doublet  distribution  obtained  using  either  his 
or  Munk's  theory. 

Karman  has  thus  suggested  a  solution  to  our  problem.  For  an  isolated 
body,  slender  body  theory  should  be  used  to  determine  the  axial  doublet 
strength.  Classical  slender  body  theory  does  not  account  for  interference 
between  the  isolated  body  and  lifting  surfaces  or  other  bodies.  A  solution 
to  this  problem  would  be  to  keep  the  finite  element  or  KctVman  approach  for 
intereference  purposes  but  to  uce  slender  body  theory  for  the  effect  of  a 
body  on  itself.  Stated  in  terms  of  influence  coefficients  this  would  mean 
that  slerder  body  theory  should  be  used  for  the  calculation  of  influence  co¬ 
efficients  of  the  bodies  on  themselves,  producing  diagonal  matrices,  and 
general  three-dimensional  wave  equation  solutions  should  be  used  for  the 
calculation  of  influence  coefficients  of  the  bodies  on  1)  lifting  surfaces 
and  2)  other  bodies.  The  interference  matrices  would  be  full  matrices. 

2.4.1  Axial  Singularity  Strength  Using  Slender  Body  Theory 

Slender  body  theory  states  that  the  flow  field  very  near  the  body  is  two- 
dimensional  and  quasi-steady.  At  each  longitudinal  station  the  flow  field 
depends  only  on  the  local  upwash  (or  si  dewash  )  and  local  body  cross-sectional 
shape.  The  axial  singularity  can  then  be  determined  from  this  two-dimensional 
quasi-steady  flow.  For  the  case  of  a  circular  cross  section  a  simple  doublet 
suffices  to  divert  the  flow  around  the  cross  section. 

us  =  27rwa^  (2.4-5) 

where  ao  is  the  local  radius  and  w  is  the  local  up-  or  si  dewash. 

For  noncircular  cross  sections  a  simple  doublet  will  not  suffice.  In  gen¬ 
eral,  a  multipole  expansion  v/ould  be  required  to  simulate  the  cross  section  in 
the  far  field.  An  alternative  approach  which  requires  no  new  basic  singulari¬ 
ties  is  to  use  doublets  spaced  laterally  a  distance  from  the  axis.  This  produces 
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the  same  effect  in  the  far  Meld  as  do  higher  order  singularities.  This 
approach  can  be  matched  to  the  proper  expansion  through  the  second  term 
(quadruple)  to  give  the  value  of  the  distance  of  the  doublet  from  the  origin. 
An  integrated  doublet,  i.e.,  two  vortices,  can  also  be  used  in  the  same 
manner. 


Consider  the  flow  about  a  body  of  elliptic  cross 
semi -major  axis  and  bQ  is  the  semi -minor  axis.  The 
gives  the  near  field  solution  as  the  two-dimensional 
Let  F  be  the  complex  potential  in  the  complex  plane. 

section  when;  aQ  is  the 
slender  body  theory 
flow  about  an  el 1  ipse. 

F 

-  -  i  w  Z  +  1  w 

(2.4-6) 

where  Z  is  the  complex  coordinate  in  the  circle  plane 
to  the  1  or  ellipse  plane  =  y  +  iz)  is 

.  The  transformation 

r 

'j 

.  Z  r  K2/Z 

C 

=  y  +  iz 

r 

-  (a0  +  b0)/2 

(2.4-7) 

K 

-  +  bo>/4 

w 

=  upwash  at  the  cross  section 

In  the  usual  procedure  for  inner-outer  matching, the  inner  solution  will  be 
expanded  in  terms  of  the  outer  variable.  Basical ly^what  is  wanted  is  the 
multipole  expansion  of  the  inner  or  slender  body  solution.  Expanding  the 
expression  for  F  and  5  gives: 
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or  retaining  one  term  past  the  doublet  gives: 


onset 

flow  doublet  quadrupole 


The  appropriate  doublet  and  quadrupole  strengths,  needed  to  represent  the 
elliptic  cross  section,  are: 

=  D2Dw  doublet 

D2D  =  2ira0(a0  +  bQ)/2  (2.4-9) 

u q  =  w27iap  {^~2 — "j  (~^~7 — ~)  quadrupole 

The  flow  about  the  ellipse  is  replaced  by  the  flow  about  a  doublet  cf  strength 
u  and  a  quadrupole  of  strength  u  .  The  approach  is  to  replace  the  doublet 

b  H 

and  quadrupole  with  two  doublets  each  spaced  a  distance  d  from  the  axis  of 
the  cross  section. 

Two  cases  will  be  considered:  1)  b  >  a,  i.e.,  major  axis  smaller  than 
minor  axis,  and  2)  a  >  b.  For  the  first  :ase,  b  >  a,  cwo  doublets  will  be 
used  (see  sketch  2.4-2).  One  doublet  placed  at  +id  and  the  other  at  -id. 


a->j 

Sketch  2.4-2 
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The  complex  potential  due  to  these  doublets  is: 


us/4tt 

F  =  f  ”  <  -?  +  -w 


Ws/4n 

c+id 


Expanding  Equation  (2.4-10)  gives: 


F  =  i  w  {  +  (ps/2n)r  ^  -  ( p^d ^/2tt ) C  ^  + 


The  term  -  usd2  is  the  quadrupole  strength.  Equating  this  to 
(2.2-4)  gives  an  expression  for  d. 


=  K 


The  second  case,  where  the  major  axis  is  larger  than  the  mi  no 
is  now  described.  In  this  case  an  integrated  doublet  is  used 
grated  doublet  gives  two  vortices;  one  positive  vortex  at  d 
vortex  at  d.  Sketch  2.4-3  gives  the  geometry. 


The  complex  potential  for  the  two  vortices  is: 

F  -  ’■»  {-< +  k 1 "  (f-Hr) 
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(2.4-10) 


..)j  (2.4-11) 

given  in 


(2.4-12) 

axis,  a  >  b. 
This  inte- 
and  one  negative 


(2.4-13) 


f 


Expanding  this  equation  gives: 


F  =  i  w  <  -  t  + 


js  *~y  *  {-¥-) 


(2.4-14) 


Equating  Equations  (2.4-14)  and  (2.4-8)  gives: 


2dr 


^s 

D2D 

d 


D2D 


(2.4-15) 


2,1  ao(ao  +  b0)/2 


^3 


T2  ^ 

~  DO 


=  \f3  K 


Equations  (2.4-15),  (2.4-12)  and  (2.4-9)  give  the  appropriate  doublet 
strength  or  vortex  strength  and  spacing  d  to  properly  represent  the  doublet 
and  quadrupole  strengths  of  the  elliptic  cross  section.  The  far  field  may 
now  be  determined  by  using  the  exact  wave  equation  solution  for  these  singu¬ 
larities  coupled  with  the  appropriate  singularity  strength  and  spacing  as 
derived  above. 


The  multipole  expansions  used  are  terminated  after  the  quadrupole.  Thus 
the  cross  section  will  not  exactly  be  an  ellipse.  Specifically,  the  ratio 
a/b  or  b/a  may  not  exactly  be  the  aspect  ratio  of  the  body.  Figures  2  and  3 
give  examples  of  body  shapes  obtained  using  two-doublets  and  two  vortices. 

In  many  instances  it  is  not  important  to  maintain  an  elliptic  cross  section 
exactly  since  the  ellipse  itself  is  an  approximation  to  body  shapes  with 
unequal  semi -major  and  semi -mi nor  axes.  The  main  idea  is  to  match  the  ratio 
of  a/b  or  b/a  of  the  ellipse  or  approximate  ellipse  to  the  physical  air¬ 
craft  cross  section.  In  some  instances  the  body  cross-sectional  shapes 
shown  in  Figures  2  and  3  are  more  representative  of  the  physical  shape 
than  is  the  ellipse.  In  any  case  it  may  be  more  important  to  correctly 
associate  the  values  of  doublet  strength  and  spacing  d  with  the  ratio  a/b 
or  b/a  for  the  approximate  ellipse. 
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2.4.2  Far  Field  Representation 


As  stated  in  the  foregoing  discussion,  exact  formulas  are  to  be  used  for 
the  effect  of  the  body  on  lifting  surfaces  and  other  bodies.  The  cross- 
sectional  area  and  shape  usually  varies  along  the  length  of  practical  bodies. 
Thus,  the  strength  and  spacing  d  must  vary  along  the  length  of  a  body. 
Numerically,  it  seems  reasonable  to  assume  that  the  cross  section  is  piecewise 
constant  over  short  lengths.  The  body  is  then  divic  d  up  into  short  elements 
(see  sketch  2.4-4).  The  two-dimensional  doublet  becomes  a  line  doublet  the 
length  of  the  element.  One  line  doublet  lies  at  n  =  +d  and  one  at  i,  =  -d. 

Each  line  doublet  or  vortex,  will,  in  general,  lie  at  a  different  value  of  d  . 
The  exact  solution  for  a  constant  strength  doublethas  not  been  derived.  How¬ 
ever,  an  expression  for  a  doublet  varying  like  e  1  ^  over  the  ele¬ 
ment  (from  f,  to  t2  where  ?c  =  U-j  +  has  been  derived  in  Appendix  B. 

This  is  the  modified  acceleration  potential  method.  The  influence  function  at 
a  field  point  then  is  the  sum  of  two  expressions,  one  for  the  doublet  element 


lyinq  at  n  -  n  +  d  and  one  for  the  doublet  eleim'nt  lyinq  at 

jo  a 


n  —  d. 
a 


1  I  ’•  Ji.  /2U 

Ers  *'  16',,  j e  b  [K(xr  "  ?is*  yr  ~  K  "  ds)’  zr  “  ?a$’ 


+  k(xr  c-j  ,  yr  (na  +  dg),  zr  ,  ...)] 
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in  Appendix 
term  Er$  g 
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S  b  j 


+  K(x  -  f.9  ,  y  -  (n,  +  d  ) ,  z 


r. 


(2.4-16) 

is  the  classic  kernel  of  lifting  surface  theory  and  is  derived 
\.  The  terms  and  t,2  are  the  element  end  points.  The 
ives  the  velocity  at  the  point  xp,  yr,  zf  normal  to  a  surface 
due  to  the  doublet  1 i ne  segment  of  unit  doublet  strength * 


If  the  short  segment  of  line  doublets  is  integrated  laterally,  as  outlined 
in  Appendix  B,  then  the  result  is  an  unsteady  trapezoidal  vortex;,  see  Sketch 
2.4.5. 

*For  convenience  the  ,<udified  doublet  (:-  -  2m)i‘s  used  in  this  analysis,  see 
Eqs.  2.4-26,  -27. 
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(2.4-17) 
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dn 


The  division  by  2d  appearing  in  Eq.  (2.4-17)  reduces  the  quadrilateral 
vortex  strength  per  unit  area  to  a  doublet  strength  (modified)  per  unit  length. 


Expanding  K  into  its  three  component  parts,  as  given  by  Eq.  (2.3-2), 

(2.4-18) 

»  ^  '  j 

d 

e(S)  =  iksf  k(S)  (x  “  "v  y  -  •••)T1/rZdn 
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-d  .  r 

-1  coA§/2U 


K-j (x  -  ••’)  +  K|S^(x  -  ?2’  •••) 


dr, 


(2.4-20) 


c(2)  -jL-  f  T* 
1  "  16^ J 


i  wAC/2U 


K2(x  -  ...)  -  K^(x  -  51 ,  ...) 


-d  l  r 

-i  uA?/2U 


K2(x  -  52»  •••)  +  K<S>(x  -  ?2,  ...)'| 

^  |  dn 


(2.4-21) 
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where  the  subscripts  r  and  s  have  been  dropped  in  Eq.  (2.4-19)  through 
(2.4-21)  for  convenience.  The  subscript  r  is  attached  to  the  receiving 
point  quantities  xr,  y  ,  zr,  yr>  etc.,  whereas  the  sending  point  subscript 
s  is  attached  to  the  sending  element  quantities,  A?s>  ^  »  ds>  etc* 

The  steady  term  E^  is  associated  with  a  trapezoidal  vortex  ot  constant 
strength  located  around  the  edge  of  the  element  and  can  be  determined  using 
the  law  of  Biot  and  Savart.  The  second  and  third  terms,  E^  and  E^  can 
be  evaluated  in  a  manner  analogous  to  that  for  and  given  in 

Eqs .  (2.3-5)  and  (2.3-6). 


A]n2  +  V  +  C1 


dn 


(2.4-22) 


(2)  1  fVitVtC2 

E  ~  7 OT  /  - 1  dn 

-d 


(2.4-23) 


+  B^n  +  C.j  =  Numerator  of  (2.4-20) 


A2n2  +  B2n  +  C?  =  Numerator  of  (2.4-21) 

This  formulation  of  1  and  E^  becomes  very  inaccurate  when  r  =  0  and 
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the  receiving  point  is  downstream  of  the  element.  The  1/r  and  1/r 
singularities  exist  only  on  the  length  of  the  element;  no  ^uch  singularities 
exist  downstream.  These  singularities  are  eliminated  in  the  numerator  in  a 
limiting  process.  However-,  when  the  numerator  is  approximated  by  a  parabola 
the  singularities  are  not  cancelled  properly.  A  way  around  this  difficulty 
is  to  f 1 c  the  entire  integral  with  parabolas. 

E^  =  ihdj  ^A3n2  +  B3n  +  C3^dri  =  TBird  (l  A3d  +  2dC3)  \2.4-24) 

-d 

d 

£l2)  =  1 ikdf  (V2  +  B4n  +  C4)dn  =  ifcd 
-d 


A4d3  +  2dC4)  (2.4-25) 
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+  ^3n  +  ^3  =  ^nte9rai')c*  °f  (2.4-20) 

A^  +  ^4n  +  C4  =  Integrand  of  (2.4-21) 

The  coefficients  A^,  Bg,  Cg,  A^,  B^,  are  calculated  by  evaluating  the 
integrands  at  three  points:  n  =  -d,  0,  d,  that  is,  the  -eft  edge,  center 
and  right  edge  of  the  element.  When  a  fieu  point  lies  close  to  one  of  these 
locations  (downstream  of  the  element)  then  r  becomes  very  small  and  accuracy 
can  be  lost  in  the  calculation.  To  avoid  this  difficulty,  a  two-point  formula 
is  used  when  r  ^O.Old.  The  offending  point  is  simply  eliminated  from  the 
calculation. 


Bgn  +  C3  =  Integrand  of  (2.4-20) 

B4n  +  C4  =  Integrand  of  (2.4-21) 

If  a  field  point  lies  close  tc  an  element  edge,  for  instance,  B  and  C  are 

determined  using  the  value  of  the  integrand  as  evaluated  at  the  other  edge 

and  the  center.  This  simplification  of  the  integration  scheme  gives  rise  to 

irregularities  in  the  normalwash  flow  field  of  the  trapezoidal  vortex  in  a 

region  downstream  of  the  element.  These  irregularities  do  not,  however, 
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invalidate  the  calculation.  Problems  involving  1/r  and  1/r  could  be 
eliminated  if  the  velocity  potential  approach  were  taken  and  new  formulas 
developed  for  the  line  and  area  integrals.  The  general  form  of  these  inte¬ 
grals  are  given  in  Appendix  B,  Eqs.  (B-l)  and  (B-2).  Details  of  the  integra¬ 
tion  over  a  small  element  have  not  yet  been  developed  but  no  problems  are 
envi s i oned . 


2.4.3  Matrix  Representation 

In  the  last  two  subsections,  expressions  have  been  derived  for  the  singu¬ 
larity  strengths  and  far  field  effects  of  short  element  lengths.  The  total 
effect,  aw,  is  obtained  by  summing  the  effects  of  all  elements.  In  matrix 
form 

vs  =  2.0  D2Dw  (2.4-26) 

(AW)  =  [E(z)]  {^z)}  +  [E(y)J  {us(y)}  (2.4-27) 
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where  aw  is  the  flow  field  caused  by  the  bodies.  Specifically,  Aw  is  the 
nornialwash  at:  i)  lifting  surface  elements,  and  2)  body  interference  elements 
(these  will  be  discussed  later)  that  lie  off  the  surface  of  the  sending  body. 


The  superscripts  z  and  y  indicate  the  doublet  orientation  in  response 
to  velocities  in  the  z-  and  y* directions.  Various  expressions  are  used  for 
E^  and  E'^';  depending  on  the  ratio  a/b  of  the  cross  section  and 
depending  on  the  location  of  the  receiving  point.  The  far  field  formulas 
(2.4-16)  and  (2.4-17)  are  used  when  the  receiving  points  fall  on  lifting 
surfaces  and  on  bodies  other  than  the  sending  body.  The  effect  of  a  body  on 
itself  is  determined  simply  by  the  inverse  of  the  expression  for  ^  given 
in  Eq.  (2.4-26).  This  means  that  the  value  of  w  used  as  a  boundary  condi¬ 
tion  is  obtained  as  the  effect  of  the  body  on  itself.  The  final  boundary 
condition  on  lifting  surfaces  and  interference  body  elements  is  the  difference 
of  w  -  aw.  The  result  of  this  difference  is  zero  for  interference  body 
elements  lying  on  the  sending  body.  Rather  than  go  through  this  operation, 
the  term  w  (excluding  the  effect  of  other  slender  bodies)  is  set  to  zero 
there.  Thus,  the  effect  of  a  body  on  itself  is  set  to  zero. 


E^z)  =  E'y-  =  0  wnen  receiving  elements  lie  on  sending  body 
w  =  0  for  interference  body  elements  on  the  sending  body 


(2.4-28) 


when  b/a  >  1 


-  Equation  (2.4-16)  with  ys  =  0 

E^  =  Equation  (2.4-17)  with  ys  =  -90° 
wr.en  b/a  <  1 


(2.4-29) 


E^  =  Equation  (2.4-17)  with  ys  -  0 
E^  =  Equation  (2.4-16)  with  ys  =  -90° 


(2.4-30) 


The  receiving  points  for  lifting  surfaces  are,  as  usual,  the  3/4-chord 
point  centered  spanwise.  The  receiving  points  for  body  interference  elements 
are  on  the  receiving  body  axis  yr  =  nar»  \  =  n3r  and  centered  longitudinally. 
For  geometrical  purposes,  bodies  are  represented  by  constant  section  tubes  to 
which  lifting  surfaces  are  attached.  In  some  instances  (e.g.,  a  tail)  this 
idealization  moves  the  surface  so  that  the  proper  aw  is  not  calculated.  A 
provision  has  been  made  to  shift  surfaces  (receiving  points)  back  to  their 
proper  position. 
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2.5  Interference 


The  first  step  in  the  detemii nation  of  interference  has  been  taken  in 
section  2.4.  Specifically,  the  incremental  normalwash  hw  due  to  the 
isolated  bodies  has  been  found,  Eq.  (2.4-23).  The  resulting  normalwash  is 

Wj  =  w  -  aw  (2.5-1 ) 

and  is  known  at  all  lifting  surfaces  and  bodies.  This  may  be  viewed  as  a  new 
normalwash  distribution  to  be  satisfied  by  the  lifting  surfuces  and  bodies  and 
need  not  be  considered  further.  Simply  replace  w  with  w~. 

The  basic  approach  to  be  taken  in  the  solution  of  this  new  modified 
boundary  value  problem  is  to:  1)  generate  an  approximate  Green's  function  for 
lifting  surfaces  in  the  presence  of  several  bodies,  and  2)  generate  a  residual 
flow  used  to  render  the  Green's  function  exact.  In  simpler  terms,  an  image 
system  is  generated  within  each  body  to  divert  the  flow  around  that  body  when 
it  is  in  the  presence  of  the  lifting  surface.  The  image  system  is  not 
completely  effective  in  doing  this,  however,  and  a  residual  flow  must  be  added. 
This  residual  potential  is  a  simple  axial  singularity  distribution,  very  similar 
to  the  axial  system  discussed  in  section  2.4. 

2.5.1  The  Method  of  Images 

The  method  of  images  is  not  new.  lennertz39  in  1927  and  later  Koeing40  were 
two  of  the  first  to  use  the  method  for  steady  flow.  The  basic  idea  of  the 
method  is  to  match  each  singularity  external  to  the  body  with  one  internal  to 
the  body  at  the  "image"  point.  The  strength  of  tne  internal  or  image  singu¬ 
larity  is  directly  related  to  that  of  the  external  singularity  strength  so  that 
no  new  unknown  distributions  are  introduced.  The  image  singularity  exists  to 
negate  the  flow  through  the  body  surface  generated  by  the  external  singularity. 

The  method  of  images  has  been  put  to  use  in  different  ways.  In  most  of 
the  aoproaches  the  residual  flow  is  ignored.  Exceptions  are  Rehorst41  and 
Wu  and  Talmadge42  who  generated  complicated  expressions  for  the  residual  flow 
fields. 
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Zlotnick  and  Robinson43  used  images  of  unswept  horseshoe  vortices  that 
are  placed  along  the  1/4-chord  point  of  a  swept  wing.  The  calculation  of 
induced  velocity  due  to  a  particular  bound  horseshoe  vortex  includes  the 
image  as  well.  The  resulting  integral  equation  for  the  span  load  is  solved 
using  a  discretized  Weissinger  approach.  Some  attempts  were  made  to  account 
for  the  effect  of  body  thickness  ratio  on  the  lift  of  the  body. 

Gray  and  Schenk44  used  a  more  approximate  method.  The  basic  approach 
was  to:  1)  determine  the  span  load  excluding  the  fuselage;  2)  form  the  image 
system  of  the  known  loading;  3)  calculate  the  resultant  upwash  distribution 
generated  by  the  image  system  on  the  wing,  and  4)  repeat  the  span  load  calcu¬ 
lation  using  the  new  upwash  distribution. 


Multhopp's  method45  and  extensions  of  it  by  Weber,  kirby  and  Kettel46 
require  a  mapping  of  thefuselage  to  a  vertical  or  horizontal  slit.  The  fuse¬ 
lage  is  effectively  reduced  to  either  a  plane  of  symmetry  or  a  segment  of  the 
lifting  surface.  In  either  case  the  resulting  simplified  problem  can  be 
solved  using  standard  methods.  Such  an  approach  ignores  entirely  the  residual 
flow  field.  This  method  cannot  be  extended  much  further  and  generalization  to 
more  complicated  configurations  seems  unlikely. 


Giesing8  recently  has  incorporated  the  method  of  images  into  the  lifting 
surface  theory.  All  horseshoe  vortices,  both  external  and  image,  possess 
sweep.  Non-midplane  configurations  are  considered.  In  addition,  the  residual 
potential  is  accounted  for  using  an  axial  doublet  distribution.  The  approach 
described  in  the  present  report  is  an  extension  of  this  method.  Spanqler  and 
Mendenhall22  have  developed  a  similar  method  for  more  general  configurations. 


Extensions  of  the  method  of  Gray  and  Schenk  have  recently  been  made. 
Borland23  has  extended  the  method  to  fuselage  cross  sections  of  elliptic  shape. 
Chou24  has  generalized  the  image  procedure  to  include  nacelles.  Also,  an 
attempt  was  made  to  generalize  the  image  approach  to  account  for  longitudinal 
variations  in  cross-sectional  area. 


Other  approaches  have  been  taken  into  account  for  wing-fuselage  interfer¬ 
ence.  For  instance,  slender  body  theory  has  been  generalized,  refined  and 
extended  to  determine  the  flow  about  very  slender  wing-body  combinations.  No 
attempt  will  be  made  to  review  all  of  the  slender  configuration  approaches. 
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Instead,  reference  will  be  made  to  an  excellent  survey  paper  written  by 
Lawrence  and  Flax47.  This  approach  has  not  been  considered  for  the  present 
method  since  more  than  just  slender  configurations  are  considered. 

2 . 5 . 2  The  Method  of  Images  for  Circular  Cross-Sections 

The  image  point  and  singularity  strength  are  determined  using  two- 
dimensional  theory.  Appendix  C  shows  how  the  Thcmpson  Circle  Theorem  is 
applied  to  obtain  the  image  point  location  and  strengths  for  three  types  of 
external  singularities  in  the  presence  of  a  body  of  circular  cross  section. 

The  three  types  of  singularities  are:  1)  the  vortex,  2)  the  source,  and  3) 
the  doublet.  The  results  are  given  in  Appendix  C  in  Equations  (C-2),  (C-3) 
and  (C-5),  respectively.  For  an  external  vortex  of  strength  r  located  at 
n,^,  the  image  strength  r j  and  location  r.j,  are  given  as: 

rI  =  "r 

nj  =  (a2/o2)n  (2.5-2) 

Cj  =  (a2/p2k 

where 

2  _  2.2 

p  =  n  +  C 

Here  it  is  assumed  that  the  axis  of  the  circular  body  is  at  the  origin  of 
coordinates.  Sketch  (2.5-1)  gives  a  graphical  description  of  these  results. 
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For  u.e  source  of  strength  a,  there  are  two  images 


Or  =0 

nT  =  (a2/p2)n 
=  (a2/p2)c 

1  (2.5-3) 

a 

r  -  -o 
A2 


nT  =  0 
x2 


For  a  doublet  of  strength  \Sy^  and  in  the  y  and  z  directions, 

respectively,  the  image  strengths  and  positions  are  given  by: 


"ly>  ■  -7  Lw(y><n2  -  ?2)  +  y^z)2?n] 

0 

=  T-  Cu(z^{n2  -  C2)  +  u(y)2?n] 

1  (2. 

nj  =  (a2/p2)n 

Cj  =  (a2/p2)c 

The  location  of  the  image  in  the  x-direction  is  identical  with  that  of  the 
external  singularity.  That  is,  the  image  singularity  matches  the  external 
singularity  both  in  length  and  position  in  the  longitudinal  direction. 


The  use  of  a  multipole  expansion  or  singularity  distribution  along  the 
axis  of  the  body  requires  that  the  surface  boundary  condition  (the  residual 
normalwash  generated  by  the  external  singularity  and  its  image)  must  be  finite 
and  continuous  on  the  body  surface.  One  of  the  basic  requirements  of  the  image 
then  is  to  render  the  boundary  condition  regular  even  when  the  external  singu¬ 
larity  lies  very  close  to  the  body  surface.  It  can  be  shown  that  such  is  the 
case.  In  the  limit  as  the  singularity  approaches  the  surface  from  the  outside, 


the  image  singularity  approaches  the  surface  from  the  inside.  In  the  limit  the 
singularity  and  its  image  render  that  portion  of  the  body  surface  lying  between 
them  a  plane  of  symmetry.  The  flow  normal  to  the  body  surface  in  that  region 
is  then  zero.  This  can  be  shown  from  the  expression  for  the  imuge  coordinates 
n j »  Ct.  Consider  the  case  ?  =  0  and  n  =  a  +  e.  Equation  (2.5-4)  then  gives: 

n,  =  a  -  e  +  0(e2) 

=  -„<*>  +  0(£2) 

4Z)  -  u<2>  +  0(c2) 

when 

n  =  a  +  e 

This  equation  shows  that  in  the  limit  as  e  ->  0  the  surface  n  =  a  is  a 
symmetry  plane.  This  is  a  two-dimensional  result.  It  can  be  shown,  however, 
that  it  also  holds  in  three-dimensions. 

As  a  summary  of  the  above  discussion  it  may  be  stated  that: 

1.  The  image  singularity  furnishes  the  major  part  of  the  disturbance 
flow  necessary  to  satisfy  the  boundary  condition  on  a  body  in  the 
presence  of  an  external  singularity. 

2,  The  boundary  conditions  for  the  residual  potential  (needed  in 
addition  to  the  image  potential)  is  everywhere  regular  even  when  the 
external  singularity  approaches  the  body  surface. 


2.5.3  The  Method  of  Images  for  Moncircular  Cross  Sections 

Borland23  has  developed  an  image  system  for  elliptic  cross  sections.  The 
image  point  is  obtained  by  transforming  the  circle  image  points.  The  transfor¬ 
mation  is  the  one  that  carries  the  circle  to  the  ellipse.  The  image  coordinates 
given  in  complex  form  are: 


nI 


Ui 


~d 
r ' 


+  (1/4)K2(A  +  iB)2 
(1/2) r2  (A  +  iB) 


(2.5-5) 
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where 


The  terms  a  and  b  are  the  semi-majoi  and  semi -mi nor  axes,  respectively. 

This  formula  is  not  valid  for  singularities  lying  outside  of  the  region 
bounded  by  the  ellipse 

2  2 

- £ - -  +  _ D _ _  =  i 

( r2/ 1 K |  -  K2|K|/r2)  (r2/ 1 K j  -  K2|K|/r2) 

since  the  image  point  then  lies  on  the  wrong  Riemann  sheet.  This  is  not  a 
serious  deficiency  since  the  image  system  is  intended  mainly  to  eliminate  the 
nonuniform  flow  field  associated  with  singularities  that  lie  close  to  the  body 
surface.  Under  this  circumstance  the  behavior  of  the  image  is  identical  to  that 
described  in  the  last  subsection.  In  a  small  region  near  the  singularity  the 
surface  of  the  body  becomes  a  plane  of  symmetry  •Jn  the  limit  as  the  singularity 
approaches  the  surface. 

As  anticipated,  the  image  approach  is  not  as  effective  for  noncircular 
cross  section  as  it  is  for  circular  cross  section.  The  image  does  not  render 
the  ellipse  a  streamline  even  in  the  two-dimensional  cross  section  plane.  The 
residual  potential  then  will  be  larger  for  noncircular  cross  sections. 

A  second  method  of  determining  the  image  point,  for  noncircular  cross 
sections  is  given  in  Appendix  C.  The  method  is  basec:  on  the  concept  of  local 
center  and  radius  of  curvature.  If  the  radius  and  center  of  curvature  are 
known  on  that  part  of  the  body  cross  section  which  lies  closest  to  the 
external  singularity,  then  an  image  point  may  be  found  using  the  circle  produced 
by  that  center  and  radius.  The  image  is  calculated  using  the  formulas  for  the 
circular  cross  section  where  the  circle  center  and  radius  vary  depending  on  the 
position  of  the  external  singularity  relative  to  the  cross  section.  Sketch  C-2 
in  Appendix  C  gives  an  illustration. 
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The  formulas  for  the  radius  and  center  of  curvature  have  been  developed 
for  ar,  elliptic  cross  section  in  Appendix  C.  If  t  is  the  vector  from  the 
center  of  the  ellipse  to  the  center  of  curvature  and  a  is  the  local  radius 
of  curvature, then: 

jr.j^^cos3e  +  ^bi^|sin3e 

-  =  (a2  sin2e  +  b2  cos2e)3/2 
ab 

where  the  parameter  6  is  related  to  y  and  z  as  follows: 

y  =  a  cos  e 
z  =  b  sin  e 

The  unit  vectors  j  and  k  are  in  the  y-  and  z-directions ,  respectively. 
Sketch  2.5-2  presents  a  plot  of  the  locus  of  C"  which  is  the  evolute  of  the 
ellipse  and  a  plot  of  a"  for  an  ellipse  of  b/a  ratio  of  0.75. 

As  the  ratio  b/s  is  reduced,  the  center  of  curvature  eventually  passes 
outside  of  the  ellipse,  thus  allowing  the  possibility  that  the  image  may  pass 
outside  of  the  ellipse.  This  method, like  the  last,  works  best  for  singular¬ 
ities  that  lie  close  to  the  body.  The  determination  of  which  method  is  best  - 
either  Borland's  or  the  method  just  described  -  remains  to  be  seen.  In  any 
case,  it  is  anticipated  that  there  will  be  restrictions  on  the  ratio  b/a  and 
on  the  position  of  tne  external  singularities  as  far  as  images  are  concerned. 

The  present  method  uses  the  local  center  of  curvature  approach  with  a 
further  restriction:  the  image  point  must  lie  in  the  same  quadrant  as  the 
external  singularity.  If  it  does  not,  it  is  ignored.  This  restriction 
eliminates  the  overlapping  of  image  sheets  within  the  body.  That  part  of  the 
image  surface  indicated  by  a  dashed  line  in  the  example  of  Sketch  2.5-3  is 
the  part  that  is  ignored. 


(2.5-6) 

(2.5-7) 
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2.5.4  The  Method  of  Images  for  Unsteady  Flow 

The  method  of  images  can  be  taken  over  into  the  unsteady  regime  essentially 
unchanged.  Attempts  to  refine  the  image  to  account  for  unsteady  effects  have 
led  to  the  conclusion  that  unsteady  effects  can  be  lumped  into  the  residual 
potential . 

Consider  the  cross  section  given  in  Sketch  2.5-4.  Attempts  to  refine  the 
image  system  so  that  unsteady  effects  are  accounted  for  include: 

1.  An  adjustment  of  the  strength  and  location  of  the  image  such  that 
the  boundary  conditions  at  A,  B  and  C  are  satisfied  in  unsteady 
flow; 

2.  Adjustment  of  the  strength  of  a  second  doublet  located  at  "0"  and 
the  position  of  the  image  doublet  such  that  the  boundary  conditions 
at  A,  B  and  C  are  satisfied  in  unsteady  flov>;  and 

3.  Adjustment  of  the  strength  of  the  second  doublet  (at  "0")  and  the 
strength  of  the  image  doublet  so  that  the  boundary  conditions  at 
A,  B  and  C  are  satisfied  in  unsteady  flow. 
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Sketch  2.5-4 

All  of  the  approaches  are  based  on  the  assumption  that  the  similarity 

variable  k  =  kr  r/c  is  small.  Consistent  with  this  approximation  is  the  use 

of  Laschka's  simplified  expression48  for  the  kernel  downstream  of  a  pressure 

doublet.  (The  expression  is  valid  if  (er/(x  -  5)  is  small.)  This  expres  ■ 

_  9  _ 0 

sion  is  expanded  in  terms  of  k  and  terms  up  to  ¥  In  If  and  1 r  d”e 

retained.  In  all  of  the  attempts  outlined,  the  changes  in  strength  anti  posi- 
tion  of  the  image  are  small  (proportional  to  k  In  ¥  and  k  )  so  tnat  pro¬ 
ducts  of  high-order  tenns  are  ignored. 

The  basic  idea  of  the  first  attempt  is  to  adjust  the  image  strength  so 
that  the  boundary  condition  at  A  (and  by  symmetry  at  C)  is  satisfied  for 
any  image  position.  The  boundary  condition  at  A  and  C  requ.-es  that  there 
be  a  stagnation  point  at  A  and  C.  The  next  step  is  to  calculate  the  value 
of  the  stream  function  on  the  body  surface  at  A  and  C.  The  image  position 
is  then  adjusted  so  that  the  value  of  the  stream  function  at  B  matches  the 
value  calculated  at  A  and  C. 

It  was  found  that  the  boundary  condition  at  A  and  C  could  not  be 
maintained  when  the  external  and  image  singularities  approached  the  circular 
surface. 
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The  second  attempt  was  meant  to  correct  this  condition.  Instead  of  using 
the  image  strength  to  satisfy  the  boundary  condition  at  A  and  C,  a  second 
doublet  was  introduced  at  the  origin  to  do  so.  Again  the  stream  function  was 
calculated  at  A  and  C  and  the  image  location  adjusted  until  the  stream 
function  at  B  matched  that  calculated  at  A.  This  method  works  except  when 
the  external  singularity  is  at  a  large  distance  from  the  body.  This  is  not  a 
serious  drawback  since  the  basic  expansion  in  terms  of  kr  r/c  fails  there 
also.  The  basic  drawback  with  this  method  as  with  the  first  one  is  that  an 
adjustment  of  distance  is  involved.  The  image  position  is  adjusted.  This 
adjustment,  however,  must  vary  in  the  longitudinal  direction;  the  correction 
to  the  position  of  the  image  varies  like  e~^  uX//U<*>.  Changing  the  position 
of  the  image  with  x  is  practical  only  when  the  frequency  is  small.  This 
type  of  correction  is  complicated  and  not  worth  the  effort. 

The  third  attempt  was  directed  toward  the  development  of  a  correction  that 
does  not  involve  adjustment  of  any  distance.  The  idea  is  to  adjust  the  strength 
of  the  doublet  at  the  origin  for  any  strength  image  so  that  the  boundary  condi¬ 
tion  at  A  and  C  is  satisfied.  The  st  ength  of  the  image  was  then  adjusted 
until  the  stream  function  at  B  was  equal  to  that  calculated  at  A  and  C. 

This  approach  has  no  real  deficiencies  except  that  the  status  of  the  boundary 
condition  at  D  and  in  between  A,  B,  C,  and  D  is  unknown.  Also,  this 
,:;erhod  like  the  second  method  fails  for  singularities  located  at  large  distances 
from  the  body. 

The  conclusions  drawn  after  analyzing  the  above  th-'cr  attempts  to  correct 
the  irraae  is  that  it  is  better  to  lump  unsteady  effects  into  the  residual 
potential.  The  last  tv/o  attempts  described  actually  use  a  correction  flow 
(the  new  doublet  at  the  origin)  to  correct  for  unsteadiness.  Rather  than 
generate  a  residual  solution  ror  each  singularity  and  image  separately,  it  seems 
logical  to  generate  the  residual  potential  only  once.  The  description  of  the 
residual  potential  is  given  in  section  2.5.6. 

2.5,5  Formation  of  tne  Influence  Matrix  with  Images 

Each  lifting  surface  element  has  its  imaqe.  For  each  unsteady  horseshoe 
vortex  exterior  to  tne  body,  there  may  be  an  image  unsteady  horseshoe  vortex 
within  the  body.  The  expression  for  the  nonnalwash  influence  matrix  for 
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elements  exterior  to  the  body  may  be  rewritten  from  Eq.  (2.3-1)  as 


D  =  l f/Kd» 

-e 

The  influence  matrix  from  the  image  element  is: 

°i  =  i /’  Kid" 

-ei 


(2.5-8) 


(2.5-9) 


The  required  image  quantities  are  determined  from  the  edge  locations  of  the 
image  element.  If  npti  are  the  coordinates  of  the  inboard  edge  of  the 
exterior  surface  element,  and  are  the  outboard  coordinates, and 

and  ~  are  the  coordinates  of  the  local  center  of  curvature,  then  the 

d 

corresponding  image  locations  are: 


2 

P1 


—2 


(n 


1 


T  - ^ 

P1 

a2  , 

+  t  U] 

Pi 


=  (n1  -na)2  +  U-, 


(2.5-10) 


Similarly,  for  nj  and  These  expressions  are  slightly  different  from 
those  of  (2.5-2)  because  it  is  no  longer  assumed  that  the  body  on  r,  is  at 
the  origin. 


eI 


COS  Y  j 

sin  Yj 


(nT  -  nT  ) 
l2  *1 

eI 

(?I  -  ?  T  ) 

l2  -1! 

eI 


(2.5-11) 
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The  sweep  angle  of  the  image  element,  Xj  is  given  by 

sin  Xj  = 

(2.5-12) 

COS  = 

Each  body  in  the  flow  may  possess  an  image  of  the  external  lifting  surface 
image.  For  the  sake  of  efficiency,  it  is  not  required  that  each  body  have  a  complete 
image  system.  That  is,  every  lifting  surface  element  need  not  have  an  image 
in  every  body.  Only  those  bodies  that  lie  close  to  the  lifting  surface  element 
are  required  to  possess  an  image  of  the  element.  Also,  if  a  lifting  surface 
element  lies  out  of  the  range  of  the  body  longitudinally,  then  no  image 
exists.  Taking  into  account  all  image  elements,  the  normalwash  becomes: 


It  is  convenient  at  this  point  to  include  the  effects  of  symmetry  and 
ground  effect.  The  effect  of  symmetry  is  to  introduce  a  system  of  lifting 
surface  elements,  plus  their  images,  on  the  left-hand  side  of  the  y  =  0 
plane.  The  n-coordi nates  of  the  sending  elements  are  changed  from  n  to 
-n  .  Similarly,  the  dihedral  angle  changes  from  ys  to  -y  .  Elements 
lying  on  the  symmetry  plane  (in  the  y  =  0  plane)  have  no  image.  The  effect 
of  the  ground  plane  is  to  introduce  a  system  of  lifting  surface  elements, 
plus  images,  below  the  z  =  0  plane.  The  ^-coordinates  of  the  sending 
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elements  are  changed  from  s  to  The  dihedral  angle  changes  from  ys 
to  -y  .  The  interaction  of  ground  effect  and  symmetry  is  to  change  n  to 
-n,  5  to  and  leave  ys  unchanged. 

D  =  0(n,t,Ys)  +  S0(-n,4,-Ys)  +  e:D(n,-c  >~ys)  +  $eD(-n,-t  .Y^) 

Note  that  the  effect  of  images  is  included  in  Eq.  (2.5-13)  because  of  its 
inclusion  in  0.  Eq.  (2.5-13)  may  be  rewritten  as 

(wR)  =  [B]  {ACp}  (2.5-14) 

2.5.6  The  Residual  Interference  How 

The  image  singularity,  of  course,  is  not  completely  effective  in 
diverting  the  flow,  generated  by  the  external  singularity,  around  the  body. 

The  reason  is  that  two-dimensional  theory  was  used  to  develop  the  image.  It 
is  surprising,  however,  how  well  it  does  work.  Wu  and  Talmadge42  have 
attempted  to  solve  the  full  three-dimensional  problem  of  a  semi -infinite  line 
doublet  in  the  presence  of  a  shape  that  has  a  circular  cross  section  and  is 
infinite  in  length*.  Their  approach  is  to  break  the  potential  into  two  parts: 
<t>i  due  to  a  doublet  and  its  image  (as  discussed  above)  and  <j>2  an  incre¬ 
mental  potential  used  to  render  the  solution  exact.  The  basic  approach  is  to 
calculate  the  flow  normal  to  the  surface  as  generated  by  the  singularity  and 
its  image,  i.e.,  as  generated  by  4>-j .  Then  a  Pourier-Bessel  series  is  used 
to  negate  this  residual  normalwash.  Basically,  what  this  series  represents  is 
a  distribution  of  singularities  along  the  axis  of  the  body.  Near  the  axis  the 
Bessel  function  can  be  approximated  by  the  functions  l/rn.  Thus,  the  series 
sum  over  n  represents  a  wave  equation  multipole  expansion.  Wu  and  Talmadge 
have  shown  that  the  residual  upwash  and  thus  the  incremental  potential  is 
everywhere  small.  The  variation  with  x  is  given  as: 

lim  <jw  ->■  0(x) 
x  -  0  c 

p 

lim  <j,2  ■*  °(x"  ) 
x  -*■  <*> 

*Wu  and  Talmadge  actually  consider  a  jpt  ar.H  not  a  body.  However,  the  basic 
solution  for  the  jet  and  body  are  very  similar. 
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Although  the  image  system  makes  the  major  contribution  towards  satisfying  the 
body  boundary  condition,  there  are  still  some  effects  not  accounted  for. 

These  effects  can  be  classified  loosely  as  follows:  1)  effects  of  bound 
vorticity,  2)  effects  of  ncncircular  cross  sections,  and  3)  effects  of 
unsteadiness. 

These  effects  produce  residual  velocities  normal  to  the  body  surface. 
These  residual  velocities,  however,  are  free  from  singularities  since  it  was 
one  of  the  main  requirements  of  the  image  system  to  eliminate  such 
singularities.  Since  no  singularities  exist,  a  solution  of  the  axial  singu¬ 
larity  type  may  be  used.  Again,  slender  body  theory  will  be  used  to  find  the 
strength  of  the  axial  singularity  required  to  eliminate  the  residual  velocity 
normal  to  the  body  surface.  In  order  to  use  the  axial  singularity  system,  an 
effective  up-or  sidewash  velocity  must  be  determined  from  the  nonuniform 
residual  onset  flow  obtained  from  the  lifting  surface  and  its  image  surface. 
Slender  body  theory  allows  two-dimensional  quasi-steady  methods  to  be  employed 
to  find  this  effective  velocity.  The  basic  problem  is  shown  ir,  Sketch  2.5-5. 


ELLIPSE  PLANE 


CIRCLE  PLANE 


z. 
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The  basic  steps  in  the  solution  are  outlined  as  follows: 

1.  Transform  the  normal  velocity  boundary  conditions,  which  are  caused 
by  nenuniform  residual  onset  flow,  to  the  circle  plane; 

2.  Solve  the  potential  flow  problem  in  the  circle  plane; 

3.  Transform  the  solution  back  to  the  ellipse  plane; 

4.  Knowing  the  doublet  strength,  find  the  effective  velocity  necessary  to 
produce  it. 

The  normal  component  of  the  residual  velocity,  w,  is 


w  =  n  •  ]T wR(y)  +  twR(z)j  (2.5-15) 

where  wR  is  the  residual  onset  flow.  This  residual  onset  flow  is  caused  by 
the  lifting  surface  elements  plus  their  images  inside  the  body  (or  bodies). 
That  is 


WR  =  ws  +  WI 


For  an  ellipse 


^  _  jb  cos  9  +  ka  sin  8 


{ 


~2 - 2 - 2 - T~ 

b  cos  e  +  a  sin  e 


The  normal  velocity  transformed  into  the  circle  plane,  w,  is 


w  =  w  y +^b/a  V1  “  D~-  (b/a)2]  cos2e  (2.5-16) 


Tn  the  circle  plane  this  normal  velocity  may  be  negated  by  a  multipole  expan¬ 
sion  centered  on  the  body  axis.  In  the  circle  plane  the  multipoles  are 
designated  by  u 

v* 


-v- -hi 


f  cos  v0 

v _ 

p'+l 


sin  v9 
v 


(2.5-17) 


where  the  subscript  v  indicates  the  order  of  the  singularity  considered. 

For  example,  v  =  1  indicates  a  doublet,  etc.  If  both  sides  of  Eq.  (2.5-17) 
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are  multiplied  by  cos  ve  and  integrated  from  e 

2u  «> 

y  w  cos  Vede  =  ^  ^  u^/r"^ 

0  7=1 

Similarly,  if  cos  ~e  is  replaced  with  sin  v9, 


=  0  ->-  2v  the  result  is: 

(2.5-18) 

the  result  is 


2ir 

j  w  si n  ~e 
0 


(2.5-19) 


Eqs.  (2.5-18,-19)  give  the  values  of  the  doublet,  quadrupole,  etc.,  strengths 
to  satisfy  the  normal  velocity  boundary  conditions.  These  singularity  strengths 
may  be  transfomed  back  to  give  the  appropriate  singularities  in  the  ellipse 
plane.  Formulas  developed  previously  (Section  2.4)  for  the  far  field  properly 
represent  the  multipole  expansion  through  the  second  term  v  =  2.  The  distance 
d  is  selected  such  that  the  proper  quadripole  strength  is  generated.  In 
Section  2.4  the  distance  d  was  a  function  of  the  cross  section  only,  i.e., 
b/a,  however,  here  d  is  a  function  of  both  cross  section  and  nonuni formity 
in  the  residual  onset  flow  field.  It  is  impractical  to  use  an  influence 
formula  that  is  a  function  of  onset  flow, thus  the  effect  of  ncnuni formity  on 
the  value  of  d  will  be  ignored.  Only  the  cross-sectional  shape  will  be  used 
to  determine  d.  Since  d  has  been  determined  previously, only  the  doublets 
need  be  considered  in  Eqs.  (2.5-18,  -19).  The  doublet  transforms  unchanged 
from  the  circle  plane  to  the  ellipse  plane.  After  substituting  m  for  w 
(Eq.  (2.5-16)) , the  expressions  for  y^  and  y^  (remember  y  =  2y)  become: 

2ir  _ _ 

y(y)  =  2a 2( l  +  b/a)  j  w  cos  0  ^/l  -  [1  -  (b/a)2]  cos2e  de  (2.5-20) 

0 


u 


(z) 


w  sin  e 


-  (b/a)2]  cos2e  de 


(2.5-21) 


Numerically, the  integrals  are  replaced  by  sums  in  the  present  method.  These 
integrals  may  be  translated  into  an  equivalent  velocity  as  follows: 
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■w . wfiifm 


'"'  *CT=T*^  '-  f^wj-  ■x  •  ^KjyjUUjW. 


’™,”™W 


„(y)  =  D2D(y)w(y) 
p(z>  =  D2D(z)w(2> 


(2.5-22) 

(2.5-23) 


where 


is  simply  the  right-hand-side  of  Eq.  (2.5-20)  divided  by  D2D^, 


The  term  w^  may  be  found  in  a  similar  manner.  The  superscripts  (y)  and  (z) 
on  the  term  D2D  are  to  indicate  that  there  is  a  difference  between  them 
for  noncircular  cross  sections. 


D2D 


D2D 


<z)  = 
(y)  _ 


2 it  a  (a„  +  bj 
o  o  o' 

2tt  b  (a  +  b  ) 

o  o  o 


(2.5-24) 

(2.5-25) 


In  Eqs.  (2.5-20)  and  (2.5-21)  the  terms  a  and  b  are  those  used  for  the 
image  calculation.  That  is,  they  are  constant  for  the  entire  cross  section. 
Thus,  the  velocity  normal  to  the  surface  is  calculated  on  the  constant  section 
tube  used  for  the  images.  The  values  of  aQ  and  bQ  found  in  Eqs.  (2.5-24) 
and  (2.5-25)  vary  along  the  length  of  the  body.  When  Eqs,  (2.5-24)  and  (2.5-25) 
are  placed  into  Eqs.  (2.5-22,  -23)  and  then  into  (2.5-20,  -21),  the  results 
are: 


,(y)  _  1  a(a  +  b) 


a  (a  +  b 
oo  o 


2ir 

J  f  w  cosede 
0 


,(z)  = 


a(a  +  b) 

a  ’(  a  +  b 
0K  0  c 


2tt 

/ 


f  w  sinode 


(2.5-26) 


(2.5-27) 


f  =  ^  -  [1  -  (b/a)2]  cos2e 

The  assumption  that  is  made  for  the  calculation  is  that  the  image  cross 
section  approximates  the  actual  one.  In  order  to  more  closely  transfer  this 
onset  flow  to  the  actual  surface,  the  value  of  a(a  +  b)/aQ(a0  +  bQ)  is  set 
to  unity.  The  velocity  normal  to  the  su 'face  of  the  ellipse  due  to  lifting 
surface  elements  and  images  is  calculated  using  D  and  Dj.  If  the  sendina 
element  is  located  along  a  body  axis  and  not  on  a  lifting  surface,  then,  E^ 
and  E^  are  used.  The  dihedral  angle  at  the  receiving  point  to  be  used  is  the 
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loeal  surface  slope  on  the  elliptic  cross  section.  The  integrals  are  evaluated 
by  calculating  w  at  various  points  around  the  circumference  of  the  cross  section 

No 

w(y)  =  1  S  w  cos  8  f  A0  (2.5-28) 

IT  i—t  V  \)  V  V 

V=1 

A  similar  expression  exists  for  w^.  The  subscript  v  indicate-  that  the 
receiving  point  is  located  at  as  follows: 


y  =  n  +  a  cos  0 
a  v 

2  =  c  +  b  sin  0 
a  v 

cos  y  =  —  sin  0  (2.5-29) 

r  p  v 
v 

sin  y  =  —  cos  0 
r  p  v 

v 

p  =  ^b2  cos20  +  a2  sin2e 

Eqs.  (2.5-26,  -27),  are  applied  to  the  influence  matrix  relating  the  normalwash 
due  to  a  lifting  surface  element  and  its  image  at  the  body  surface. 


As  it  stands  in  Eq.  (2.5-14)  {w^}  is  partitioned  into  two  parts:  1) 

the  normalwash  at  the  liftina  surface  elements,  wD,  and  2)  the  normalwash 

K  (e) 

at  the  body  surface  at  various  meridian  angles  (e^),  w^  ' 


(  O 

V** 


\ 

; 


In  order  to  obtain  the  average  z-  and  y-velocities  at  the  bodies.  Eqs.  (2.5-26,-27) 
must  be  applied 

2tt 

D^z)  =  I  f  V  cos  edo  (2.5-31) 

0 

2^ 

D<y)  =  i  J  DQf  sin  ado  (2.5-32) 

0 
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Then  the  matrix  {  wg}  becomes: 


ACpj  (2.5-33) 

It  is  important  to  notice  that  the  cross-sectional  shape  (and  size)  used  in  the 
determination  of  the  effective  up-  and  sidewash  (Eqs.  (2.5-31,  -32))  must  be 
the  same  as  that  used  to  determine  the  images.  This  cross  section  is  constant 
for  the  entire  body  length.  The  effective  velocity  obtained  from  these 
calculations,  however,  is  applied  to  the  actual  body  cross-sectional  shape  and 
size  through  the  use  of  the  local  influence  function  D2D  and  the  setting  of 
the  term  a(a  +  b)/a0(aQ  +  b  )  to  unity. 

Thus  far  the  effect  of  lifting  surface  elements  and  their  images  on  bodies 
has  been  considered.  Attention  is  now  turned  to  the  interference  effects  of 
bodies  on  all  parts  of  the  configuration,  i.e.,  other  bodies  and  lifting 
surfaces. 

The  residual  interference  flow  is  to  be  generated  using  the  same  singular¬ 
ities  as  those  used  for  the  isolated  body  (section  2.4,  Eqs.  (2.4-16,  -17)). 

The  only  differences  are:  1)  that  the  distribution  of  elements  will  be  differ¬ 
ent,  2)  the  normalwash  is  determined  at  points  on  the  cross-sectional  surface 
rather  than  at  the  body  axis  (even  though  the  results  are  applied  to  the  exact 
surface) . 

The  distribution  of  elements  for  interference  could  be  identical  to  that 
for  the  isolated  body  for  simplicity.  However,  this  is  inefficient.  The 
singularity  strengths  for  the  isolated  body  are  known  and  do  not  add  to  the 
number  of  unknowns  in  the  problem  as  do  the  interference  elements.  For  this 
reason  a  generous  number  of  elements  may  be  allocated,  for  the  isolated  body, 
to  properly  describe  the  body  radius  distribution,  the  greatest  variation  of 
which  is  near  the  leading  and  trailing  edges.  The  greatest,  variation  in  the 
residual  flow  may  be  elsewhere,  specifically  near  the  body/lifting-surface 
intersection. 


(_\ ) 
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i 
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The  body  radius  must  equal  the  constant  value  used  in  the  image  calcula¬ 
tion.  The  reason  is  that  the  image  singularity  surfaces  exist  within  this 
constant  section  tube.  Thus,  two  idealizations  for  the  body  are  used:  1)  tc 
generate  the  flow  field  of  the  isolated  body,  the  exact  radius  distribution  is 
used,  and  2)  to  account  for  interference,  a  tube  of  constant  cross  section  is 
used  along  with  the  exact  body  shape. 

The  residual  flow  field  due  to  lifting  surfaces  and  images  has  nonuniform¬ 
ities  across  the  body  surface.  Thus,  the  averaging  of  this  nonuniform  velocity, 
given  by  Eqs.  (2.5-26,  -27)  is  required.  For  the  effect  of  bodies  on  other 
bodies,  however,  the  nonuni ^ormi ties  are  small  and  the  normalwash  may  be 
evaluated  on  the  axis. 


The  influence  matrix  relating  the  normalwash  at  a  field  point  due  to  an 
interference  body  element  is  formally  the  same  as  that  for  an  isolated  body 
except  that  the  values  of  y^,  y^  are  not  known. 

twn)  =  C'£(Z)]  {y<z)}  +  [E(y)]  {y^}  (2.5-34) 


E(z)  -  E(y) 

rs  rs 


=  0 


when  the  receiving  element  lies  on  the  sending  body 
except  when  r  =  s 


E^  =  D2D(z)  when  r  =  s 
E(y)  =  020^  When  r  =  s 


(2.5-35) 


b/a  >_  1 

EjZ)  =  Eq.  (2.4-16)  with  y$  =  0 
E^y)  =  Eq.  (2.4-17)  with  y$  =  -90° 


(2.5-36) 


b/a  <  1 

E^  =  Eq.  (2.4-17)  with  ys  «  0 
E"'  =  Eq.  (2.4-16)  v.'ith  y$  =  -90° 


(2.5-37) 


The  receiving  points  for  the  lifting  surface  elements  are  the  3/4-chord 

points.  The  receiving  points  for  the  bodies  are  on  the  body  axis  at  , 

r  a 

zr  =  c  and  centered  longitudinally  on  the  element. 

In  summary  then,  the  residual  interference  is  handled  in  the  following 
manner: 

1.  The  up-  or  sidewash  at  bodies  due  to  lifting  surfaces  and  their 
images  is  determined  using  an  averaging  technique  on  the  normalwash 
at  the  body  surface.  The  body  surface  used  for  this  purpose  is  a 
constant  cross-section  tube. 

2.  Interference  body  elements  are  used  to  negate  the  up-  and  sidewash 
given  in  (1)  and  also  the  up-  and  sidewash  generated  by  other  bodies. 
The  interference  body  elements  are  distributed  differently  from 
isolated  or  slender  body  elements.  This  is  done  to  reduce  the 
number  of  unknowns  in  the  problem  and  to  place  interference  elements 
where  they  are  most  needed. 

2.5.7  The  Use  of  Surface  Singularities, for  the  Residual  Interference  Plow 

In  the  last  section  axial  singularities  were  used  for  the  determination 
of  the  residual  potential.  An  alternate  method  is  presented  here.  Specifically, 
the  normalwash  boundary  condition  can  be  satisfied  using  a  distribution  of 
either  unsteady  horseshoe  vortices  or  quadrilateral  vortices  on  the  body 
surface.  Iri  Part  I  of  this  report,  unsteady  horseshoe  vortices  were  used 
on  the  body  surface  to  account  for  all  of  the  body/lifting-surface  interference. 
(Slender  body  elements  were  used  to  determine  aw.)  For  the  present  method, 
however,  surface  singularities  on  the  body  surface  may  be  used  to  account  only 
for  the  well-behaved  residual  flow.  Thus,  both  an  image  system  and  surface 
singularities  can  be  used.  The  number  of  elements  on  the  body  surface  need 
not  be  as  large  as  for  the  method  of  Part  I  since  the  residual  flow  is  well 
behaved.  However,  the  introduction  of  surface  singularities  greatly  increases 
the  number  of  unknowns  in  the  problem  and  the  computational  effort.  It  is 
anticipated  that  this  method  will  be  used  only  for  difficult  configurations  or 
if  a  very  high  degree  of  accuracy  is  required.  It  may  also  be  desirable  to 
have  such  an  approach  for  cross-checking  purposes. 
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The  normalwash,  Wy,  to  lifting  surface  elements  lying  on  body  surfaces 
is  determined  in  the  same  manner  as  the  Wy  for  axial  interference  elements. 
The  only  difference  is  that  the  velocity  normal  to  the  surface  elements  is 
determined  instead  of  z-  and  y-wash  at  the  body  axis.  For  a  single  body,  Wy 
at  these  surfaces  is  zero  since  it  is  assumed  that  the  slender  body  axial 
distribution  has  diverted  the  onset  flow  around  the  body  surface  (exceDt  for 
interference  flow).  If  two  or  more  bodies  exist  (including  symmetry  and 
ground  effect  bodies),  then  one  body  may  induce  a  nonzero  Wy  on  another 
body. 


Eventually  it  may  be  desirable  to  use  quadrilateral  vortices  that  are 
inclined  to  the  flow  field  so  that  more  details  of  the  body  surfaces  may  be 
treated.  In  such  a  case  it  is  necessary  to  use  a  more  general  formula  than 
the  one  developed  in  Appendix  B.  An  expression  for  a  pressure  doublet  inclined 
to  the  flow  was  developed  by  Berman'*9.  Such  an  expression  integrated  over  an 
element  would  give  an  inclined  horseshoe  vortex.  This  type  of  singularity  is 
unsatisfactory  for  the  body  surface  since  a  wake  trails  back  from  each  point 
on  the  surface.  Such  wakes  would  then  thread  in  and  out  of  the  body  surface. 

A  similar  expression  for  a  velocity  potential  doublet  would  be  desirable  to 
produce  an  inclined  quadrilateral  vortex  without  a  wake.  If  any  wake  did 
exist  it  could  be  added  at  the  end  of  the  body  using  unsteady  horeseshoe 
vortices. 

2.5.8  The  Image  System  for  Wing-Tail  Configurations 

The  method  of  images  requires  at  of  constant  cross  section.  The 
image  system  trails  back  on  the  inside  of  the  body  without  any  changes  in 
lateral  position.  The  application  of  the  image  system  concept  to  bodies 
of  variable  cross  section  with  special  refere.  .e  to  the  wing-tail  problem  is 
the  subject  of  this  subsection. 

If  a  lifting  surface  system  is  attached  to  a  body  at  only  one  location, 
then  a  fairly  simple  solution  to  this  problem  exists.  At  the  intersection  of 
the  body  lifting-surface  system  an  average  body  cross  section  is  selected  for 
use  with  the  image  system.  The  effects  of  the  lifting  surface  on  the  actual 
body  c3n  be  determined  on  this  average  surfac  In  other  words,  the  normal - 

wash  to  the  actual  body  surface  generated  by  the  lifting  surface  and  image 
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system  can  be  calculated  on  the  average  surface.  The  results  calculated  on  the 
average  surface  are  then  applied  to  the  actual  surface  as  if  they  were  calcu¬ 
lated  on  the  actual  surface.  If  the  residual  flow  field  is  reasonably  uniform 
(laterally),  then  very  little  error  should  be  introduced  with  this  method. 


This  approach  is  very  similar  to  that  of  Woodward  for  steady  flow.  His 
bas'ic  approach  is  to  use  slender  body  theory  for  the  exact  body  shape  in  a 
uniform  flow  to  determine  the  axial  doublet  and  source  distribution.  These 
singularities  are  then  used  to  calculate  an  o.'set  flow  to  the  lifting  surfaces 
on  the  wing.  Lifting  surface  elements  are  then  placed  on  an  idealized  body 
shape  of  constant  cross  section  to  account  for  the  rest  of  the  interference. 

If  a  wing-tail  configuration  is  considered,  then  the  tail  must  be  attached  to 
the  idealized  body  shape.  This  idealized  shape  usually  takes  on  the  radius  of 
the  body  as  it  is  at  the  wing-body  intersection.  This  may  cause  errors  in  the 
tail  loads  for  two  reasons:  1)  the  relative  positions  of  the  left  and  right 
horizontal  stabilizers  as  well  as  the  fin  may  be  changed  considerably  to 
accommodate  the  idealized  fuselage  shape,  and  2)  the  loads  in  the  presence  of  a 
very  large  diameter  fuselage  are  different  from  those  in  the  presence  of  a 
small  diameter  fuselage.  Reference  21  giv^s  examples  of  the  Woodward  approach 
for  the  B-58  and  other  aircraft  configurations. 


A  possible  solution  for  the  wing-tail  problem  is  to  use  one  idealization 
for  wing-body  interference  and  second  idealization  for  the  tail-body  interfer¬ 
ence.  This  cannot  be  done  with  the  Woodward  approach,  but  it  may  be  possible 
with  the  present  approach.  For  instance,  the  body  idealization  suitable  for 
the  wing-body  intersection  could  be  used  to  determine  wing- on- wing ,  wing- on- 

hnrlv  anH  unnn  nn  ‘♦’ail  i  nfov'fov'onrn  T  tin  Krxrlw  i  70  mi  f  a  4-l-*rv 
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tail-body  intersection  could  be  used  to  determine  tail-on- tail ,  tail-on-boay 
and  tail-on-wing  interference.  Sketch  2.5-6  gives  a  graphical  example  of  this 
approach. 


The  relative  positions  of  the  left  and  right,  lifting  surface  systems  are 
now  correct  for  the  local  influence  of  these  surfaces  on  themselves.  The 
effect  of  the  wing  on  the  wing  is  very  accurate  as  well  as  the  effect  of  the 
tail  on  the  tail.  The  effect  of  the  wing  on  tue  tail  deserves  some  spec1'-! 
consideration.  The  approach  outlined  accounts,  in  a  crude  way,  for  the  face 
that  the  wake  follows  the  flow  field  around  the  body.  Consider  for  instance, 
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IDEALIZATION  FOR  DETERMINING 
THE  EFFECT  OF  THE  WING  ON 
ITSELF  AND  OTHER  PARTS  OF  THE 
AIRCRAFT 


IDEALIZATION  FOR  DETERMINING 
THE  EFFECT  OF  THE  TAIL  ON 
ITSELF  AND  OTHER  PARTS  OF 
THE  AIRCRAFT 


Sketch  2.5-6 

the  trailing  vortex  at  the  wing  root.  This  vortex  follows  the  fuselage  and 
crosses  the  tail  plane  at  the  fuselage-tail  intersection,  i.e.,  at  the 
horizontal  tail  root.  With  the  approach  just  outlined,  the  same  effect  is 
achieved  by  placing  the  wing  and  tail -root  sections  at  the  same  points 
laterally.  If  more  information  is  known  about  the  wake  location,  either 
laterally  or  vertically,  the  tail  position  could  be  adjusted  to  more  accu¬ 
rately  reflect  the  relative  geometries  of  the  wing  wake  and  tail.  The  idea 
is  to  select  an  idealization  that  places  the  wing  wake  properly  relative  to 
the  tail  since  it  is  the  wake  that  generates  the  strong  wing-tail  coupling. 

Adams  and  Sears33  presented  a  fairly  simple  formula  for  calculating  the 
lateral  streamline  displacement  about  a  body  of  revolution  neglecting  vortex 
sheet  roll-up.  It  is 

r2U)  =  r2(?1)  -a2(5-,)  +  a2(c)  (2.5-38) 

Here  r(c)  is  the  lateral  position  of  a  streamline  (where  the  radius  is  aQ(0) 
that  had  its  origin  at  rU-j)  (where  the  radius  is  aQ(c1)).  Such  a  formula 
could  be  used  to  determine  where  wing  trailing  vortices  emanating  at  rU-j) 
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will  strike  the  tail  plane.  The  use  of  this  approach  would  cause  practical 
difficulties  for  the  lattice  method*.  In  the  lattice  method  the  trailing 
vorticity  is  idealized  in  discrete  vortex  lines.  These  vortices  must 
impinge  on  the  tail  plane  at  strip  edges  (if  wing  and  tail  are  copianar  or 
near  copianar).  The  strips  at  the  horizontal  tail  must  then  be  adjusted  so 
that  they  align  with  the  wing  edges  after  the  transformation  (2.5-38)  is 
applied. 

By  far,  the  largest  factor  in  determining  the  position  of  the  wing  wake  at 
the  tail  is  the  vortex  sheet  convection  and  roll-up.  Although  the  present 
method  does  not  determine  this  effect,  there  is  no  fundamental  reason  why  it 
could  not  be  included  at  a  later  date.  It  seems  appropriate,  at  the  present 
time,  to  delay  the  inclusion  of  the  second-order  wing-tail  idealization 
improvements  (outlined  in  this  section)  until  a  full  study  of  all  of  the 
second-order  effects  can  be  made. 

2.6  Body  Force  and  Moment  Distributions 


Heretofore  this  report  has  concerned  itself  with  the  determination  of  the 
singularity  system  that  has  been  substituted  in  place  of  the  aircraft.  This 
section  deals  with  the  loading  obtained  from  the  singularity  distribution  and 
subsequently  the  generalized  forces  from  the  loading. 

On  all  lifting  surfaces  the  loading  is  equal  to  the  singularity  strength 
and  thus  no  further  discussion  is  needed.  However,  for  the  bodies  involved, 
the  situation  is  much  more  complicated. 

The  objective  is  to  develop  a  simple  set  of  expressions  for  the  load 
distribution  on  a  body,  in  the  presence  of  a  singularity  distribution.  Both 
the  singularity  distributions  interior  to  the  body  (axial  doublets,  vortex 
quadrilaterals  and  images)  and  those  exterior  to  the  body  (lifting  surface 
elements  and  axial  singularities  on  other  bodies)  generate  lifting  forces  on 
bodies.  Unlike  lifting  surfaces,  bodi»  ^  have  nonzero  thickness.  A  change 
in  potential,  say  from  the  bottom  to  the  top  side  of  a  body,  causes  a  net 

*Such  a  method  would  be  fairly  easy  to  apply  to  loading  function  methods.  The 
spanwise  distribution  of  vorticity  would  simply  be  modified  by  Eq.  (2.5-38). 
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force.  The  flow  fields  of  all  singularity  distributions  give  rise  to  a 
variation  in  <V,  +  Wp  across  the  body  and  thus  generate  a  force  or.  the  body. 
Such  a  variation  is  not  experienced  by  lifting  surfaces  since  they  are  of  zero 
thickness.  The  only  variation  is  caused  by  the  local  potential  jump,  A<f>, 
which  generates  the  local  ACp.  The  contribution  to  the  loading  from  the 
slender  body  singularities  is  easily  obtained  using  the  theory  of  Miles.  The 
contributions  due  to  the  rest  of  the  s ingulat  i ties  will  be  determined 
numeri cally. 

All  of  the  singularities  used  to  model  the  aircraft  can  be  expressed  in 
terms  of  unsteady  pressure  doublets.  It  is  sufficient,  then,  to  develop  an 
expression  for  the  force  and  moment  on  a  body  due  to  a  point  pressure  doublet. 
This  expression  is  then  applied  to  the  equivalent  pressure  doublet  distribu¬ 
tion  that  has  been  substituted  for  the  actual  singularity  distribution. 


2.6.1  Reduction  of  All  Singularity  Distributions  to  Pressure  Doublet  Distributions 


The  reduction  of  the  unsteady  horseshoe  vortices  to  doublets  is  straight¬ 
forward  since  the  horseshoe  vortex  is  simply  an  integration  of  the  pressure 
doublet  in  the  spanwise  direction  over  a  short  element.  The  unsteady  horseshoe 
vortex  is  simply  replaced  by  one  or  more  point  doublets  depending  on  the 
accuracy  required  and  the  distance  between  the  horseshoe  vortex  and  body. 
Currently, in  the  present  method,  only  one  point  pressure  doublet  is  used. 

The  total  doublet  strength  of  an  unsteady  horseshoe  vortex  of  strength  ACp 
is  £Cp*A  where  a  A  >s  the  area  of  the  element  associated  with  the  horseshoe 
vortex.  In  the  analysis  rf  Appendix  D  it  is  assumed  that  all  point  pressure 


doublet  strengths  are  given  t,y  .  CpsA  even  chough  the  doublet  may  be  due 
an  axial  element  ano  not  a  lifting  surface  element 


to 


Consider  now  the  axial  doublet  of  strength  u  associated  with  the 
modified  acceleration  potential  (Appendix  A,  Eqs.  (A-19,  -20)).  The  points 
of  concentrated  prc^-.-rc  (delta  functions),  which  correspond  to  the  origin  of 
an  equivalent  pressure  .  mo'iec, are  obtained  from  Lqs.  (A-21)  and  (A-27)  as 


-1' 

e 


(2.6-1) 


68 


where 


AQ  =  <5(n  —  na,  ?_  ca)u(c)  (2.6-2) 

Here  6  is  a  two-dimensional  Dirac  delta  function,  and 


Thus  s  is  the  center  of  a  doublet  element  running  from  to  tg.  When 

£,  is  outside  of  the  range  shown,  jumps  to  the  center  of  the  next  element. 

Thus,  5c(c)  is  a  step  function  at  ?  =  c-j  and  c2-  For  an  isolated  element 

Sc  =  0  C  '  ^  or  r  >  f_2 

Also  for  an  isolated  element 

A(T  =  0  5  <  or  6  >  r>2 

that  is,  the  doublet  strength  drops  to  zero  off  of  the  element  as  a  step 
function  in  5. 

_  i»S  /u„ 

The  product  AQ  e  then  is  also  a  step  function  at  c  =  £]  and  ?2. 

The  derivative  of  this  step  function  is  zero  everywhere  except  at  £  =  ^ 
and  at  which  point  it  is  infinite.  The  results  are  two  delta  functions. 

^  i  'x'A5/2Um 

ACp6A  -  6(n  —  na>  5—  ta)<$(?  -  5-j)  ue  ”  (2.6-3) 

%  -i  wA^/2u\ 

-<5(r,  -  na,  c  “  ?a)<5(5  “  C2)  y  e 

where 

f  =  «c  -  «1  =  «2  -  ?c 

These  expressions  represent  two-point  pressure  doublets:  one  at  6  =  ^  of 
strength  v^e1  and  one  at  5  =  ?2  of  strength  une~"' 

Sketch  2.6-1  presents  a  graphical  illustration  of  the  singularity  (the 
imaginary  part  is  omitted  for  clarity). 
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Z,i 


Sketch  2.6-1 


Usually  the  doublet  elements  are  not  isolated  but  are  linked  end-to-end 
down  the  center  of  a  body.  In  this  case  each  end  point  possesses  two  point 
pressure  doublets  whose  strengths  may  be  summed. 


ACp6A  =  Uj  = 


He 


i  wACs/2Uo 


—  u 


s-r 


-i  w&S$ml/2Ua 


(2.6-4) 


where  s  indicates  the  element  subscript,  i.e.,  element  s  lies  directly 
aft  of  element  (s  —  1)  on  the  same  body. 

The  same  analysis  is  valid  for  quadrilateral  vortex  rings  except  that  horse¬ 
shoe  vortices  are  used  in  place  of  pressure  doublets.  At  a  particular  end  point 
of  an  element  there  exist  two  horseshoe  vortices:  one  emanating  from  the  element 
upstream  cf  the  point  and  one  from  the  downstream  element.  These  two  horseshoe 
vortices  are  each  replaced  by  one  or  more  pressure  doublets  just  as  the  lifting 
surface  horseshoe  vortices  were  replaced  by  pressure  doublets.  The  fact  that 
the  upstream  element  may  be  larger  in  span  than  the  downstream  element  ooses 
no  problem.  The  element  width  only  enters  in  the  determination  of  the  element 
area.  The  expression  for  the  equivalent  pressure  doublet  strength  for  quadri¬ 
lateral  vortices  is, 
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(2.6-5) 


AC  6A  =  AQ  2d  e 
p  s  s 


i  0)A £  /2U 

S 


-  A^S-l2ct5-le 


-i  uA?s_1/2Uo 


where  A^2d  may  be  considered  a  doublet  strength,  u,  since  2d  is  the  width 
of  the  trapezoidal  vortex.  In  fact,  for  use  with  bodies,  the  doublet  strength 
un  is  used  instead  of  a£[.  Thus,  Eq.  (2.6-3)  may  be  used  for  both  doublets 
and  trapezoidal  elements. 


The  loads  due  to  slender  body  elements  can  be  obtained  directly  using 
Miles’ slender  body  theory. 


AC 


2 


w^L  + 


(2.6-6) 


(2.6-7) 


L  =  it(2  aaQ/3x  +  i2kraQ/c) 


(2.6-8) 


M  =  n  a^/c 

where  a  is  the  local  body  width  and  b/a  is  the  ratio  of  the  semi- 
o 

minor  to  semimajor  axis.  The  terms  w2  and  wv  indicate  upwash  and  sidewasb, 
respectively.  The  relation  between  the  axial  loading  aF/ax,  and  the  pres¬ 
sure,  ac2,  ACy  is: 


aiyq 

ax 


=  46z2a0 


(2.6-9) 


sx 


Ac  2a 
y 


b 

o  a 


(2.6-10) 

\ 


To  simplify  the  logic  of  the  computer  program,  these  loads  are  converted  to 
point  pressure  doublets  and  treated  like  all  of  the  other  pressure  doublets. 

9F  Jq 

(aC  6A)  =  u  =  Ax  (1  +  b/a)  (2.6-11) 

p  z  sz  3x 
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ax  {1  +  a/b) 


(2.6-12) 


<VA)y 


3F  /„ 

-  y/q 


The  term  Ax  is  the  slender  body  element  length.  The  factors  (1  4  b/a)  for 
us  and  (1  +  a/b)  for  ps  are  applied  because  of  the  cross-sectional 
effect! vity.  For  example,  aycircle  (a/b  =  b/a  =  1)  transmits  only  half 
the  point  pressure  doublet  strength  to  the  body;  thus,  the  doublet  strength 
has  to  be  multiplied  by  (1  +  a/b)  =  (1  +  b/a)  =  2. 


2.6.2  Calculation  of  Force  and  Moment  on  a  Body  in  the  Presence  of  a  Point 
Pressure  Doublet 

The  expression  for  the  pressure  at  all  points  in  a  flow  field  due  to  a 
point  pressure  doublet  is  given  in  Appendix  D  as: 


C  (X  y  z)  -  ACpU,n,?)6A  ixM(x-g)  9  )e"UR) 
VX’y,ZJ  4tt  e  W  ( — R — j 


(2.6-13) 


where  ACpSA  is  the  strength  of  the  point  pressure  doublet  Lawrence  and 
Flax1'7  have  integrated  the  steady  version  of  this  equation  making  use  of  the 
following  assumptions.  1)  doublet  lies  at  a  large  distance  from  body,  2)  body 
is  circular,  3)  the  pressure  doublet  is  oriented  in  the  z-direction  and  also 
lies  in  the  z  =  0  plane,  and  4)  only  the  force  in  the  z-direction  is 
desired.  The  bas^'c  idea  of  the  Lawrence  and  Flax  method  is  to  integrate  the 
steady  version  of  Eq.  (2.6-13)  around  the  cross  section  for  an  arbitrary  value 
of  x  thus  producing  a  distribution  of  force  along  the  body  axis.  A  second 
method  to  be  discussed  later,  integrates  first  in  the  x-direction  (analytically) 
and  then  integrates  numerically  around  the  cross  section  to  produce  the  total 
force  (and  moment).  This  method  does  not  produce  a  distribution  for  a  single 
doublet  since  the  total  force  and  moment  are  lumped  at  the  same  x-direction 
as  that  of  the  doublet. 


A  generalization  of  the  steady  flow  method  of  Lawrence  and  Flax  is 
oerived  in  Appendix  6.  Specifically,  the  force  distribution  in  both  the  y- 
and  z-directions  is  derived  for  a  doublet  located  at  any  point  in  the  flow 
field  which  has  an  arbitrary  orientation.  If  N*  is  a  unit  vector  in  the 
direction  of  the  doublet,  and  if  n  and  r,  are  the  coordinates  of  the 
doublet,  then: 
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2  2 

AC  6A  fa: 


>(F/q)/ax  =  -J— 


where 


(2  2 

3V 

1  _“7~ 


Rg  =  (x  -  z)2  +  e2r2 


r  > 
a  o 


(2.6-14) 


”  =  “A  +  Vr 


ra  ■  "  ”a>2  *  'z  -  «a>2 

and  where  r  is  the  radial  distance  from  the  circle  origin  to  the  doublet, 

V  ^ 

and  i„  and  i‘  are  unit  vectors  in  the  e-  and  r  -directions ,  respectively. 

y  i  a 

If  this  result  is  integrated  to  obtain  the  total  lift  on  a  body  of  constant 
cross  section,  then 


aC6A  sr 

F/q  =  _l_-°  („9r6-vr> 


r,  > 

a  o 


(2.6-15) 


The  total  moment  on  the  body  is  zero  since  the  distribution  is  symmetrical 
in  x  for  the  steady  case.  An  amazing  result  is  observed  from  Eq.  (2.6-15). 
The  total  force  obtained  by  integrating  the  approximate  formula  of  Eq. 

(2.6-14)  is  exactly  correct  as  seen  by  comparison  with  equation  F-13  of 
Appendix  F.  Equation  (2.6-14)  was  derived  for  values  of  R  which  are  large, 

a 

however  the  integrated  value  is  exact  for  all  values  of  Rg  where  rg  >  aQ. 


A  similar  expansion  for  small  rg/a  gives  the  following  result  for 


uuuuicii  onao  lie  Wi  oi  I  iff  one  uOuy. 

_  AC  6A  82a2 

3(F/q)/9X  =  -Jp-  -g—  (uel0  +  urir) 

Ra 


r.  <  a„ 
a  o 


(2.6-16) 


where 


=  (x  —  z)2  +  62a2 


Again  if  this  equation  is  integrated  over  a  constant  cross-sectional  body. 


AC  AA 


a  ■  -f-  <Uot0  -  vv> 


r  <  a 
a  o 


(2.6-17) 


This  result  corresponds  to  the  exact  values  obtained  in  Appendix  F. 
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Such  encouraging  results  for  the  steady  case  (body  of  circular  cross 
section)  indicate  that  a  similar  approximate  analysis  might  work  for  both 
the  unsteady  case  and  the  case  of  noncircular  cross  sections.  A  derivation 
for  the  unsteady  case  (circular  body)  is  given  in  Appendix  G.  For  the  case 
where  the  doublet  lies  outside  of  the  body,  i.e.,  r  >  a„,  the  result  is: 

a  0 


<r/,)/»x  -  6w«x-s>-Ra] 


sTe 


VV 


ra  >  aQ  (2.6-18) 


where 
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T  =  I _! _  +  J* 
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R. 
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Here 


X  =  =  2k 

e  c 


(2.6-19) 


For  the  case  where  the  doublet  lies  on  or  inside  the  body,  i.e.,  r,  <  a  ,  the 

a.  U 

result  is: 


S(f/q)/3X  .  ^  62a2elA[H(x-«)-Ra]  j  ^  +  n 


“a  ^ 


II  Vs  +  Vr 


(2.6-20) 


r  <  a 
a  -  o 


These  results  can  be  integrated  over  a  constant  cross-sectional  body  to 
give  the  total  force.  The  term  F/q  •  can  be  integrated  analytically  to 
give 


AC  <$A  a“  /.  \  r 

Fo/q  =  Hr-  7  Y  (r)  Hj 


(2) 


(k) 


(2.6-21) 


-  (2) 

where  k  =  u>r  M/U  and  where  H)  '  is  the  Hankie  function  of  the  second 

C.  00  J 

kind  and  zeroth  order.  The  result  for  Fr/q  must  be  integrated  numerically. 
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Just  as  in  the  steady  case  the  total  forces  integrated  over  a  body  of 
constant  cross  section  are  exact.  The  results  of  (2.6-21)  and  the  results 
of  numerically  integrating  Fr/q  are  plotted  in  Sketch  2.6-2.  The  exact 
results  are  obtained  by  the  method  presented  in  Appendices  D  and  E  (for  the 
total  force). 


The  results  for  a (F/q)/3x  are  given  above  in  terms  of  u„i_ 
UpY*  A  conversion  to  the  y-z  plane  gives: 


t  >  and 


where 
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3(F  /q)  a(F^Vq)  e(F^Vq) 
ax  z  ax  \  ax 


(2. 6-22a) 


3(F/q)  9(F<y)/q)  9(F<y)/q) 
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ex  y  ax  z  ax 


(2.6-22b) 


3(F^z)/q)  2  a(F  /q)  <  2  a(Fr/q) 
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N  -  jNy  +  kNz 


tanel  =  —7 


?  -  z. 


a(F  /q) 

6 

3X 


a(F/q) 


ax  0 


a(Fr/q) 

ax 


ilEM.i 

ax  r 


65 


The  circumferential  integration  for  a  circular  cross  section  has  been 
performed.  In  the  present  method  the  force  distribution  is  given,  not  at 
every  value  of  x,  but  at  a  set  of  points  which  correspond  to  the  midpoints 
of  a  set  of  body  elements  (the  slender  body  elements).  The  value  of  the 
circumferential  integration,  a(F/q)dX,  must  be  averaged  over  these  elements. 
If  not,  then  errors  could  be  introduced  due  to  local  variations  over  the 
elements.  The  average  of  a(F/q)ax  over  an  element  is  obtained  by  integrating 
it  and  dividing  by  Ax.  Ultimately  the  total  force,  F/q  on  the  element  is 
required,  thus  the  division  by  ax  can  be  left  out. 


F/q 


-f 


xa+(ax/2) 


xa-(ax/2) 


a(frg) 

ax 


dx 


where  xA  and  ax  are  the  element  center  and  length,  respectively.  In 
equation  (2.6-22)  the  elements  that  go  to  make  up  a(F/q)/3x  are 
">(F^/q)/3x,  a(F^yVq)ax,  and  3 (F^y Va/ax .  These  are  all  integrated  in 
the  same  manner  as  is  3 (F/q) ax  producing: 


/ 


xa+(ax/2)  3(p(z)/q) 


3x 


dx 


xa-(Ax/2) 


(2.6-23a) 


p(y)  xa+(Ax/2)  3(p(y)/q) 

-X_=  [  - 1 - dx  (2.6-23U) 

q  I  ax 

\-(h*/2) 
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xa+(ax/2)  sF(yi/q) 


xA-( ax/2) 


ax 


dx 


Also,  required  for  future  analysis  is 


/q  = 


/q- 


F^y)/q 


This  case  is 


(2.6-23C) 


(2.6-23d) 
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In  addition,  future  calculation  will  require  moments  on  the  elements.  For 
this  method  (circular  bodies)  the  moments  are  zero. 

The  above  simplified  approaches  for  the  steady  and  unsteady  cases  (which 
will  be  called  Approach  I)  are  based  on  expansions  in  terms  of  a  /R,  or 

0  3 

ra/Rg-  suc^  they  are  approximate  even  though  the  total  forces  are  exact 
(the  total  moment  is  correct  for  the  steady  case  also). 

For  the  case  of  an  elliptic  cross  section,  the  simple  expansion  techniques 
used  above  and  derived  in  Appendix  G  (Approach  I)  fail  to  converge  to  sufficient 
accuracy  for  the  total  inttjrated  values.  Therefore,  these  results  are  not 
presented.  A  different  approach  (Approach  II)  will  have  to  be  taken  for  the 
case  of  noncircular  cross  sections  . 

The  basic  idea  of  the  second  approach  is  to  reverse  the  order  of  integra¬ 
tion  in  finding  the  total  force  on  the  body.  That  is,  Eq.  (2.6-13)  is 
integrated  in  the  x-direction  first  and  then  integrated  around  the  cross 
section.  The  advantages  of  this  method  are  that:  1)  unsteady  flow  can  be 
considered  for  bodies  of  arbitrary  cross  section  since  Eq.  (2.6-13)  can  be 
integrated  analytically  along  the  body,  and  2)  only  a  single  term  arises  from 
the  expression  and  not  an  entire  distribution  for  each  pressure  doublet.  The 
disadvantage  of  this  method  is  that  although  the  force  and  moment  are  correct, 
the  detailed  distribution  is  not  given.  (See  Sketch  2.6-3  for  a  graphical 
representation  of  approach  II.)  In  a  practical  problem  there  are  many 
pressure  doublets  and  approach  II  produces  a  distribution  since  there  is  a 
distribution  of  pressure  doublets.  The  distribution,  however,  is  more  abrupt 
than  the  actual  since  the  force  aria  moment  are  concentrated  at  one  point. 

The  integration  of  Eq.  (2.6-13)  to  obtain  the  total  force  and  moment  is 
carried  out  in  Appendix  D.  The  integration  in  the  x-direction  is  performed 
analytically  in  terms  of  two  Hankel  functions.  The  results  are: 
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m  -  -i  §  kH,(2>(k) 

=  m  §  bi(2)(id 

fr' '(F)  =  Hankel  function  of  second  kind 
k  =  2krM  r/c 
k ,  =  «c/?U 

'  oo 


N  =  direction  of  pressure  doublet 


5^  -Vr  TTT 
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'fp  -  direction  of  force  F 
i  =  direction  of  radius  r 

r  =  \  ^iy  ~  n)c  +  (2  -  r)2  =  j(y  -  n)  +  k(z  -  t) 
q  =  pU2  =  dynamic  pressure 

and  where  the  contour  integral  is  to  be  taken  around  the  cross  section.  The 
cross-sectional  shape  under  consideration  is  the  one  possessed  by  the  body  at 
the  longitudinal  location  of  the  point  pressure  doublet.  Currently,  for 
simplicity,  the  constant  crof.s  section  designated  by  a  and  b  is  uced 
instead  of  the  local  values,  aQ  and  bQ.  Appendix  D  shows  that  higher- 
order  moments  diverge  and  thus  cannot  be  used. 

The  integrals  found  in  Eqs.  (2.6-24,  -25)  can  be  evaluated  analytically 
for  only  the  simplest  cases.  One  such  case  is  the  one  considered  earlier  in 
this  section,  i.e.,  a  circle  in  steady  flow.  Appendix  F  gives  the  details  of 
this  integration.  The  results  are: 


£  -  AC  <SAI 
q  P 


(2.6-26) 
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•  N  Pressure  doublet  inside  body 

Ni(a/rJ^  Pressure  doublet  outside  body 
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(2.6-27) 


where  a  is  the  circle  radius,  ra  ir  the  radial  distance  from  the  pressure 
point  to  the  circle  center,  and 

ff|  =  if  sin  ( <f>  -  2e)  -  J  cos  U  -  2e") 
where  *  and  0*  a^e  defined  implicitly  as  follows: 

N  =  k  sin  <{>  +  "{  cos  4> 

r"  -  Ir’J  {k  sin  o’  +  t  cos  F} 

0  0 
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For  other  cross-sections  the  integrals  in  Eqs.  (2.6-24,  -25)  must  be 
done  numerically.  Appendix  E  gives  the  details  of  a  successful  numerical 
scheme  for  evaluating  these  integrals. 

The  basic  idea  is  to  break  up  the  contour  integrals  into  a  sum  of 
integrals  over  short  surface  segments.  On  each  of  these  elements  the  values 
of  f(F)  and  p(k)  may  accurately  be  held  constant.  The  integrals  then 
become: 


AC  sA  v 

r/q^\  «kj)  / 


(2.6-28) 


ac  <$a 

H/p  =  1  p(kj>  /  h 

j=1  -e. 


.v  a 

N  n  •  ip  ds 


(2.6-29) 


The  integral  over  the  small  element  extending  from  -e  to  <-e  in  the 
moment  equation  (2.6-29)  is  very  simple  since  no  singularities  exist. 
However,  the  1/r  singularity  in  the  force  integral  requires  special  treat¬ 
ment  and  its  accurate  evaluation  is  difficult.  Several  procedures  were  tried 
which  ignored  the  curvature  effect.  These  proved  to  be  inaccurate  and  had  to 
be  abandoned.  The  major  part  of  Appendix  E  is  devoted  to  the  evaluation  of 
the  element  integral  in  Eq.  (2.6-11).  It  is  convenient  to  break  T/q  and 
M/q  irto  their  component  parts 
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Similar  expressions  may  be  written  for  I4/q . 
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To  test  the  accuracy  of  the  formulas  in  Appendix  E,  comparison  is  made  to 
the  analytic  solution  for  a  circle  given  in  Eqs.  (2.6-26,  -27).  Figure  4 
presents  results  for  the  cases 

(a)  TF  =  N  =  k 

(b)  Tp  =  N  =  j- 

(In  Figure  4  as  well  as  Figures  5,6,7  and  8,  ACp6A  is  assumed  to  be  unity.) 
The  numerical  results  are  accurate  even  near  the  body  surface.  It  is  near  the 
surface,  both  inside  and  outside,  where  greatest  difficulty  is  encountered  in 
obtaining  accurate  solutions.  There  is  a  tendency  for  the  internal  solution  to 
drop  below  its  correct  constant  value  near  the  surface.  Similar  calculations 
were  carried  out  for  elliptic  cross  sections.  One  such  calculation  is  given 
in  Figure  5.  The  internal  values  stay  constant  except  near  the  surface  where 
a  slight  variation  is  noticed.  In  Figure  5  the  pressure  doublet  is  located 
in  the  2=0  plane.  In  Figure  6  a  similar  plot  is  showr.  Here,  however, 
the  doublet  lies  in  the  y  -  0  plane.  No  irregularities  are  observed. 

The  reason  for  the  difference  in  accuracy  between  Figures  5  and  6  is  the  fact 
that  near  the  surface  the  curvature  is  higher  in  Figure  5.  Figure  7  is  pre¬ 
sented  to  furnish  a  ch<-  'k  on  the  unsteady  case.  The  pressure  doublet  is 
located  on  the  axis  o'  j  circle.  In  this  case  the  integrals  of  Eqs.  (2.6-24, 
-25)  are  easily  performed  since  r  is  not  a  function  of  the  surface  coordi¬ 
nate  s  and  can  be  taken  out  of  the  integral.  The  results  are: 
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M  ACn5A 

q  =  —5 —  p(k)aT 
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In  Figure  7,  T  and  ACp£A  are  assumed  to  be  unity.  The  comparison  of  the 
analytic  and  calculated  results  is  essentially  perfect.  Notice  that  for 
the  steady  case,  the  pressure  doublet  of  strength  ACpSA  produces  a  lift  on 
the  body  of  strength  ACp6A/2.  For  other  cross-sectional  shapes,  this 
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factor  changes.  For  an  ellipse  in  steady  flow  with  a  z-oriented  pressure 
doublet  located  inside  the  cross  section,  the  variation  of  the  force  with 
b/a  is 


q  =  aVa  (;ti) 

As  the  cross  section  becomes  flat  (b  -+  0),  the  bocy  becomes  more  effective 
until  b  =  0  at  which  point  the  body  becomes  like  a  lifting  surface  and  the 
force  is  equal  to  the  pressure  doublet  strength,  AC^A.  Figure  8  presents 
a  comparison  of  this  formula  with  calculated  numerical  results.  The  agree¬ 
ment  is  very  good  even  for  very  small  values  of  b/a  or  a/b. 


2.6.3  Fffects  of  Symmetry,  Ground  Effect  and  Images 

The  method  outlined  so  far  considers  only  a  sinqle  point  pressure  doublet. 
In  actuality  for  every  pressure  doublet  located  on  a  lifting  surface,  there 
may  be:  1)  one  for  symmetry,  2)  one  for  ground  effect,  3)  one  for  the  inter¬ 
action  of  symmetry  and  ground  effect,  4)  images  in  bodies  associated  with  the 
lifting  surface,  and  5)  images  due  to  symmetry  and  ground  effect.  The  z-force 
on  body,  b,  due  to  pressure  doublet  located  at  £,n>?  oriented  normal  to 
a  surface  of  dihedral  y  is: 
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•yF^U,-n,-0  -  sin  .-n.-c) 


(2.6-31) 


where 


F(z)  .  ':z:/<) 
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(2.6-32) 


and  a  similar  result  holds  for  F^,  etc.  The  term  NAB  indicates  the 

y  *- 

bodies  associated  with  the  pressure  doublet,  i.e.,  those  bodies  that  have 

0/ 

an  image  of  the  pressure  doublet  within  them.  The  terms  uz  and  u 
arise  from  the  fact  chat  the  image  point  pressure  doublet  has  a  modified 
strength  shown  in  Appendix  C ,  Eq.  (C-5).  If  this  equation  is  rewritten 
in  terms  of  real  variables  •ed  the  circle  radius  a  ir  replaced  with  the 
radius  of  curvature  a,  t,.  i, 
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(2.6-33) 
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and  where 


za  ae  i’-t  coordinates  of  the  local  center  of  curvature  (for 


a  circle  these  <•<:.-  the  coordinates  of  the  body  axis) 


An  express!  ,r  <  t  .  lar  to  Eq.  (2.6-31)  may  be  written  for  the  force  in  the 
y-direction,  i  c.,  !>  l/.  The  miy  changes  to  be  made  are  that  is 

replaoed  by  ry  ^  <mr.  f ^  ^  replaced  by  f^.  For  the  moment  calcula- 
'  1""'i  ■ 1  ■  r  r>  r  ■  I  i ' «?  «  ’ 1  witu  Ms.  Iq.  (2.6-31)  may  be  made  to  hold  for 


pome  pressure  r 

St  till.'!  t 


f‘r"  o,  "livhiiVj  or,  l.ody  axes  by:  1)  omitting  all  images, 
'•  -  oriented  doublets  and  .  =  -90°  for  y-oriented 


So  far,  two  basically  different  methods  have  been  described;  one  for 
bodies  of  circular  cross  section  and  one  for  all  other  bodies.  The  first  of 
these  methods  produces  a  distribution  of  forces  on  all  elements  of  the  body 
due  to  a  single  doublet.  If  we  characterize  all  doublet  strengths  as  AtpM, 
then  the  total  force  at  an  element  is: 

Fzjb)  =  (ACp6A)s  (2.6-34) 

where  r  represents  the  receiving  element  and  s  respresents  the  sending 

doublet.  A  similar  expression  exists  for  FyI^.  In  matrix  notation 

lr 

{FZ-[b)l  =  [FZ(b)]  (ACp6A>  (2.6-35a) 

{FY-[b)l  =  [FY(b)]  UCpSA>  ( 2 . 6-35b ) 

The  superscript  b  stands  for  body.  There  is  a  column  matric  Fz|b^  and 
FY-[b)  for  each  body. 

For  the  second  method,  each  doublet  produces  a  point  force  and  moment  on 
the  body.  As  illustrated  in  Sketch  2. 6-4,  the  longitudinal  location  of  the 
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Sketch  2.6-4 
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center  of  force  is  at  the  location  of  the  point  pressure  doublet.  Again 
the  slender  body  element  system  is  used  to  keep  track  of  the  loads  on  the 
bodies.  If  the  longitudinal  location  of  a  point  pressure  doublet  lies  between 
the  leading  and  trailing  edges  of  a  slender  body  element,  the  force  and 
moment  are  added  to  the  load  on  that  element.  To  smooth  out  this  bookkeeping 
procedure,  it  is  assumed  that  the  force  is  uniformly  distributed  over  a  short 
segment,  on  the  axis  of  the  receiving  body,  whose  length  is  equal  to  the  length 
of  the  element  upon  which  the  point  pressure  doublet  lies.  Thus,  the  force 
and  moment  due  to  one  particular  pressure  doublet  may  be  spread  out  over 
several  slender  body  elements.  Sketch  2.6-4  gives  a  graphical  example.  The 
load  and  moment  due  to  the  pressure  doublet  is  distributed  over  slender  body 
elements  3,  4  and  5.  The  center  of  load  on  these  three  elements  lies  at  a, 
b,  and  c,  respectively.  The  fact  that  the  load  is  off  center  contributes 
to  the  moment  on  the  slender  body  elements. 

2 . 6 . ^  Redistribution  of  Body  Loads  for  Alternate  II 

The  previous  sections  have  described  the  two  methods  of  determining  the 
force  and  moment  distributions  for  bodies.  For  alternate  II  (general  body 
cross  section)  the  contribution  of  each  point  pressure  doublet  emanating  from 
each  singularity  ard  image,  etc.,  is  allocated  to  the  proper  slender  body 
element.  The  final  result  for  a  body,  then,  is  a  distribution  of  forces 
among  its  slender  body  elements.  The  distribution  thus  obtained  approaches 
the  exact  distribution  as  the  body  radius  becomes  small.  For  thick  bodies, 
however,  the  resulting  distribution  is  somewhat  distorted  even  though  the 
total  force  and  moment  are  correct.  The  only  way  to  avoid  this  distortion 
is  to  consider  that  each  pint  pressure  doublet  contributes  to  all  slender 
body  elements  as  the  result  of  its  longitudinal  distribution.  To  do  this, 
however,  Eq.  (2.6-13)  would  have  to  be  integrated  first  around  the  contour. 

Such  an  integration  has  been  done  for  circular  cross  sections  but  has  not  been 
done  for  elliptic  cross  sections.  Until  such  time  as  this  integration  can  be 
performed,  alternate  II  will  be  used  for  such  cases.  There  is  a  Wa.v  to  allevi¬ 
ate  some  of  the  distortion  of  alternate  II  and  that  is  by  redistribution. 

After  <he  effects  of  all  point  pressure  doublets  have  been  allocated  among 
the  c"''  •  body  elements,  assume  that  these  forces  have  arisen  solely  from 

an  axi  .  ,nt.  pressure  doublet  distribution  within  that  Dody.  Only  a 
relatively  small  number  of  point  pressure  doublets  are  now  considered.  It  is 
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assumed  that  Eq.  (2.6-20)  gives  a  good  approximation  to  the  shape  of  the 
longitudinal  distribution  of  the  ibrce  along  the  body,  even  though  it  only 
holds  for  circular  cross  sections.  The  body  radius  used  is  (aQ  +  bQ )/2 , 
i.e.,  the  average  of  semi-major  and  semi-minor  axes.  The  redistribution  thus 
outlined  will  not,  then,  affect  the  total  lift  of  the  distribution.  This 
redistribution  is  then  applied  to  each  of  the  slender  body  elements.  Even 
though  the  results  obtained  using  this  method  show  good  agreement  with 
experimental  data,  it  is  not  proposed  here  as  the  final  answer  to  the  distri¬ 
bution  problem.  Further  work  is  required  to  integrate  Eq.  (2.6-13)  around 
arbitrary  cross  sections. 

2.7  Generalized  Forces  and  Aerodynamic  Parameters 
2.7.1  Generalized  Forces 

The  pressures,  forces  and  moments  obtained  on  lifting  surfaces  and  bodies 
along  with  their  appropriate  displacements  are  outlined  as  follows: 

1.  aC  ,  h,  (normal  to  lifting  surfaces) 

P 


2. 

3(Fz/q)/ax, 

h2 

(z-direction  on  bodies) 

3. 

3(F  /q)/3X, 

hy 

(y-direction  on  bodies) 

4. 

3(Mz/q)/3X, 

dhz/dx 

(z-direction  on  bodies) 

5. 

3(My/q)/3x, 

dhy/dx 

(y-direction  on  bodies) 

The  virtual  work,  6w-,  done  by  these  pressure  forces  and  moments  during 

J 

a  virtual  displacement  is: 


•If 


<5h  ds 


+  9 


/ 


3(Fz/qh 

— 5T^  {,,2  dE  + 


r  3<Fv/q)i 

J  -  h~ 3  <1? 


B. 


3C 


y 


T3(Mz/q)  /dh\  r  a(M  /q).  /dh\  ) 

J  —re1 6wj d5  +  J  ~ir 1 6  a? 7 dc  ( 

B.  B.  | 


(2.7-1) 
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where  q  is  the  dynamic  pressure.  The  integration  limit  R.S.  indicates  that 
only  the  lifting  surfaces  on  the  right-hand-side  of  the  aircraft  are  considered. 

[l  if  lifting  surface  lies  in  plane  of  symmetry  (e.g.,  vertical  fin) 


|2  otherwise 

The  integration  limit  B.  indicates  all  bodies  lying  on  the  right-hand  side 
of  the  aircraft. 

[l  if  body  lies  on  plane  of  symmetry  (e.g.,  a  fuselage) 


2  otherwise  (e.g.,  a  nacelle) 

The  integral  over  the  lifting  surface  L.S.  represents  a  series  of  integrals 
over  each  of  the  lifting  surfaces  which  go  to  nu-ke  up  the  total  configuration. 
The  value  of  G  for  each  of  these  surfaces  may  be  different.  A  similar 
argument  is  valid  for  the  integral  over  the  bodies  B.,  thus  the  value  of 
g  may  vary  from  body  to  body. 

The  values  of  a(F/q)/3x  and  a(M/q)/Dx  are  the  original  values  obtained 
before  redistribution. 

iwo  definitions  for  the  generalized  force  will  be  introduced,  i.e.,  Q.. 
_  *3 


and  Q 


°A^q  .  «q, 


(2.7-2) 


(2.7-3) 


Here  q"..  is  the  i-th  generalized  coordinate,  c  the  reference  chord  length, 
A  the  total  reference  area  and  s  the  reference  semi  span.  With  the  intro¬ 
duction  of  generalized  coordinates,  the  displacements  may  be  written: 


l 


h  =  cXVl 
hz  "  ‘XVz,. 

hy  =  ‘SVy,  <2-7-4) 


Introducing  the  virtual  displacements  of  Eq.  (2.7-4)  (and  its  derivative) 

into  Eq.  (2.7-1),  and  equating  (2.7-1)  to  (2.7-2,  -3)  gives  the  results  for 

the  generalized  forces  Q..  and  Q. .. 

J  1  J 


ACpfids 


p(Fz/q)  -  f  3(F  /q)  . 

+g  J  — l — 1  f  d5  +  /  — ^ i  f  d5 
B  9?  i  o  3?  yi 


/ 

B. 


U 


dS 


B. 

id?  */- 

B. 


8(V^.i  df*i  ■  r»(yq>, 1  dfyj 


3? 


d? 


d? 


(2.7-5) 


(2.7-6) 


2.7.2  Aerodynamic  Parameters 

It  is  dpsirahle  and  sometimes  necessary  to  generate  conventional  aero¬ 
dynamic  data.  Such  data,  in  addition  to  being  useful  in  itself,  provide  an 
excellent  check  for  the  computer  Drogram  and/or  specific  cases  to  be  r by 
it. 


The  local  normal  force  coefficient  and  pitching  moment  coefficient  about 
the  local  1/4-chord  point  are: 
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s*Wh-*»v.- 


l 

ii 


C  =  7-  f  AC  d? 

n  c  J  p 


(2.7-7) 


chord 


cn|  =  -^r  /ACp  (5  -  51/4)d5  (nose  up)  (2.7-8) 


where  c  is  the  local  chord  length.  The  local  center  of  pressure  is: 

-MS,) 

c-p-Re  =  "RiT^T+  0,25 

“Im(c  ) 

C,p-Irn  =  "TmlcT  +  0,25 

n 

The  total  vertical  and  side-force  coefficients  on  lifting  surfaces  are 

r  -  0  +  <s) 


(2.7-9) 


C.  = 


(1  -6) 
A 


/  ccndn 

(2.7-10) 

R.S. 

(f)  /  ccndf; 

R.S. 

(2.7-11) 

The  total  vertical  and  side-force  coefficients  on  bodies  are 

C_  =  Y./  g  /  - de 

d§ 


O  +  «)  ,  f 

3(Fz/q) 

A  9  J 

95 

B. 

O  -«)  g  f 

9(Fy/q) 

A  9  J 

95 

B. 

(2.7-12) 


(2.7-13) 


The  pitching  and  yawing  moment  coefficients  on  lifting  surfaces  taken  about 
the  point  XM  are: 
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CM  =  tLLJl  J  ic-cm  _  ccr(c1/4  -  XM)[dn  (nose  up) 
Ac  R.S. 


(2.7-14) 


■'H 


(|)  J  |c2cm-ccn(c1/4-XM)|d?  (nose  right)  (2.7-15) 
Ac  '  R.S. 


For  bouies 

C  =  ii 


if'U 


f  ( 3 (M  /q) 

v  a(p2/q)  ) 

f\  « 

(«  XM)  ^  | 

d5  (nose  up)  (2.7-16) 


JN 


=  iL^ii  g  f  j!^ 

»  *  B.<  95 


3(F  /q)  I 

(5  -  XM)  — —  d5  (nose  right)  (2./-17) 


The  total  rolling  moment  for  the  aircraft  is: 

C£  “  ~  yiT”  i  J  ccnndn  +  I  /  ccncd?  +  /  35  "a 


n  d5 


>R.S. 


R.S. 


B. 


f3(F  /35) 

J  ~b - 43d5 


35  a 


(right  wing  down) 


(2.7-18) 


B. 
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3.0  CALCULATED  RESULTS 


3. 1  Parameter  Variation  Studies 

Before  discussing  the  correlations  of  the  present  method  with  other 
theories  and  experimental  data,  a  series  of  studies  will  be  made  to  determine 
the  effect  of  varying  some  of  the  important  parameters  of  the  problem.  These 
studies  not  only  show  how  results  are  changed  by  changing  a  parameter  but  also 
indicate  the  range  of  validity  of  the  present  method. 

3.1.1  The  Effect  of  Body  Radius 

Figure  9  presents  a  comparison  of  spanwise  lift  distribution  or.  a  wing 
(exposed  portion  only)  that  is  attached  to  a  circular  fuselage  whose  radius 
has  been  varied.  The  exposed  wing  is  the  same  in  all  cases  —  it  is  attached 
to  the  fuselage  center  —  and  is  oscillating  in  pitch  at  a  reduced  frequency 
of  unity.  The  fuselage  is  at  zero  incidence  and  thus  the  upwash  generated 
by  the  fuselage  is  not  present.  Thus,  this  figure  gives  only  the  interference 
effect  of  the  fuselage.  The  effects  of  zero  radius  and  infinite  radius  are  the 
same.  Consider  the  lift  coefficient  at  the  wing  root.  When  the  diameter-to- 
chord  ratio  is  increased  from  zero  to  0.125,  then  the  lift  drops.  A  further 
increase  to  0.50  produces  an  increase  in  lift  coefficient  over  the  zero  value 
case.  Further  increases  eventually  bring  the  fuselage  effect  back  to  the 
zero  diameter  case. 

The  effect  of  the  fuselage  when  it  too  is  given  an  angle  of  attack  is 
nrch  different.  The  effect  is  monotonic  in  fuselage  diameter- to-chord  ratio 
and  in  the  limit  of  infinite  radius  the  lift  coefficient  is  twice  the  zero 

A  •>  +  ll*,l  Urt 

U  I  UIIIC  (/Cl  vuiuc, 


3.1.2  The  Effect  of  Cross-Sectional  Shape 

The  effect  of  fuselage  ellipticity  is  illustrated  in  Figure  10.  Specifically, 
the  distribution  of  lift  coefficient  for  a  wing-fuselage  combination  is  pre¬ 
sented  for  the  aspect  ratio  b/a  ranging  from  zero  to  infinity.  The  semi-width 
a  is  constant  and  equal  to  0.2s  for  all  calculations.  The  case  b/a  =  0.0 
is  spc.ial  and  is  obtained  by  replacing  the  body  by  a  flat- lifting  surface. 

The  case  b/a  =  ■*>  is  also  special  and  is  obtained  by  placing  a  plane  of 
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symmetry  at  the  wing  "oot  (this  case  is  given  the  designation  of  "wing  alone"). 
All  of  the  intermediate  cases  which  are  not  special  should  lie  within  the 
envelope  formed  by  the  two  special  cases  described  above.  The  curve  corres¬ 
ponding  to  the  lowest  value  of  b/a  (=0.25)  does  not  fall  within  the  envelope. 

The  curve  corresponding  to  the  highest  value  of  b/a  (=2.0)  still  lies  inside 
the  envelope  but  comes  very  close  to  passing  outside  of  it.  The  upper  bound 
of  validity,  for  b/a,  of  the  present  method,  is  then  somewhere  near  or 
below  2.0.  The  lower  bound  of  validity  for  b/a  is  somewhere  between  0.5 
and  0.25.  If  configurations  are  to  be  analyzed  that  fall  outside  of  this 
range,  then  it  is  recommended  that  either  the  body  be  replaced  by  a  lifting 
surface  (b/a  <  0.5)  or  a  plane  wall  for  (b/a  >  2.0).  An  additional  alter¬ 
nate  is  to  place  lifting  surface  singularities  on  the  body  surface  in  addition 
to,  or  in  place  of,  the  image  system. 

The  use  of  lifting  surface  singularities  on  the  body  surface,  as  described 
in  Part  I  of  this  report,  does  not  always  increase  the  accuracy  of  the  solution. 
Figure  11  presents  a  comparison  of  the  present  method  with  the  method  of  Part  I 
for  the  wing-fuselage  combination  (b/a  =  2.0)  of  Figure  10.  Two  separate 
idealizations  are  used  to  describe  the  fuselage.  The  first  idealization  uses 
three  panels  or  strips  to  describe  a  quarter  of  the  fuselage.  The  resulting 
calculation  falls  below  that  of  the  present  nethod.  When  the  number  of  fuselage 
elements  is  doubled,  the  agreement  is  improved.  The  actual  spanwise  distribu¬ 
tion  of  lift  coefficient  probably  lies  between  the  present  method  and  the  six 
strip  idealization  since  Figure  10  indicates  that  the  present  method  may  be 
high  by  a  slight  amount  and  Figure  11  indicates  that  the  method  of  Part  I  may 
be  low. 

A  third  calculation  which  has  not  yet  been  discussed  is  the  use  of  lifting 
surface  panels  and  images  in  the  present  method.  This  calculation  has  not  been 
done  for  the  steady  case  of  Figure  11;  however.  Figure  14  shows  such  a  calcula¬ 
tion  for  the  unsteady  case,  and  is  discussed  below. 

3 . 1 . 3  The  Effects  of  Frequency 

A  comparison  of  the  lift  coefficient,  as  calculated  by  the  present  method 
and  the  method  of  Part  I  for  a  wing-fuselage  combination,  is  given  in  Figure  12. 
The  motion  of  the  wing  may  be  described  as  a  modified  pitch  where  the  normalwash 


corresponds  to  that  of  steady  pitch.  The  effects  of  frequency  enter  through 
the  influence  matrix  [Dy].  The  fuselage  has  a  circular  cross  section  and  is 
not  in  motion.  Agreement  between  the  two  methods  is  good  for  the  two  common 
frequencies  (kr  =  0.0  and  1.0).  Intermediate  frequencies,  kr  =  0.1  and  0.5 
are  also  plotted  to  show  the  variation  with  frequency. 

3.1.4  A  High-Wing  Configuration 

A  comparison  between  the  results  of  the  present  method  and  the  method 
of  Part  I  is  shown  in  Figure  13.  The  wing-fuselage  combination  associated 
with  Figures  10  through  12  is  used  except  for  a  small  extension  of  the 
exposed  area  of  the  wing  at  the  root  section.  The  total  semispan  is  held 
constant  and  in  order  tc  join  the  wing  to  the  fuselage;  a  small  extension 
to  the  wing  surface  is  required  at  the  wing  root.  The  comparison  between  the 
two  methods  shows  the  same  trends  observed  in  Figure  11,  i.e.,  that  the 
results  of  the  method  of  Part  I  lie  below  those  of  the  present  method. 

Also  shown  in  Figure  12  are  results  obtained  for  the  midwing  case.  These 
results  lie  below  those  for  the  high-wing  case. 

Figure  14  presents  a  similar  comparison  except  that  the  frequency  is  1.0 
instead  of  0,0  as  in  Figure  13.  The  wing  is  pitching  about  its  root  leading 
edge.  The  body  is  at  rest.  In  this  case  the  real  parts  are  in  good  agreement. 
The  imaginary  parts,  however,  differ  near  the  wing-fuselage  juncture.  The 
imaginary  lift  coefficients  for  the  present  method  lie  below  that  calculated 
by  the  method  of  Part  I.  This  seems  to  be  contrary  to  what  is  expected  since 
the  results  of  the  present  method  have  been  greater  than  the  results  of  the 
method  of  Part  I  in  Figures  II  and  13.  A  third  calculation  is  presented  to 
help  decide  which  solution  is  most  accurate.  Specifically,  the  third  solution 
consists  of  placing  lifting  surface  elements  on  the  body  surface  in  addition 
to  the  images  of  the  present  method.  The  results  agree  almost  perfectly  with 
the  present  method  near  the  wing-fuselage  intersection.  However,  over  the 
outboard  half  of  the  wing  this  calculation  agrees  best  with  the  method  of 
Part  I.  However,  it  may  be  observed  that  the  difference  among  any  of  the 
solutions  is  fairly  small  over  the  outboard  half  of  the  span.  This  figure 
indicates  that  the  image  method  is  most  accurate  for  this  case. 
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3.1.5  End  Plating  Effects 


In  all  of  the  calculations  presented  thus  far,  the  configurations  have 
been  symmetrical  about  the  y  =  0  plane.  If  the  left-hand  wing  half  is 
omitted  from  the  calculation,  then  the  effect  of  body  end-plating  can  be 
ascertainea.  Such  a  calculation  is  presented  in  Figure  15.  Again  the  same 
midwing  configuration,  as  found  in  Figures  10  through  12,  is  used  except  for 
the  missing  left  wing  half.  It  is  expected  that  the  end-plating  effects,  as 
obtained  by  the  present  method,  are  slightly  smaller  than  they  should  be.  Two 
reasons  are  offered  for  this:  1)  the  image  within  the  ellipse  is  cut  off  once 
it  crosses  the  y  =  0  plane,  and  2)  the  interference  singularities  cannot 
account  for  a  flow  that  is  antisymmetric  across  the  body.  If  the  body  were 
circular,  then  the  image  would  fully  account  for  the  end-plating  effect  and 
the  onset  flow  to  the  body  would  be  symmetric  across  the  body  thus  allov'ing 
the  axial  singularities  to  properly  account  for  interference.  However,  for 
the  elliptical  rase,  some  of  the  end-plating  effect  may  be  lost  due  to  the 
approximated  and  abbreviated  image  system.  The  resulting  residual  flow  should 
make  up  for  this  loss  in  end-plating.  However,  flows  that  are  antisymmetric 
to  the  body  surface  cannot  be  accounted  for  because  of  the  truncation  of  the 
multipole  series.  If  one  further  term  were  added  to  this  series,  then  this 
effect  could  be  accounted  for  even  for  highly  elliptic  cases  (b/a  1.0). 

A  comparison  of  lift  coefficient  distribution  as  calculated  by  the  present 
method  and  the  method  of  Part  I,  is  presentee  in  Figure  15.  The  wing  is  given 
an  angle-of-attack  ot  1.0  radians.  The  body  incidence  is  zero.  The  two 
methods  appear  to  be  in  very  close  agreement.  Figure  11  presents  the  identical 
case  except  for  symmetry.  In  Figure  II  the  calculated  results  ot  the  present 
method  lie  everywhere  above  those  of  the  method  of  Part  I.  However,  ir. 

Figure  15  the  lift  distribution  as  calculated  by  the  present  method  falls 
below  that  of  the  method  of  Part  I  over  most  of  the  span.  This  shows  that, 
relative  to  the  method  of  Part  I,  the  present  method  has  lost  some  of  its 
end-plating  effect. 

The  end-plating  effect  of  a  body  whose  ellipticity  b/a  falls  below 
unity  (b/a  <  1.0)  is  a  different  matter.  In  this  case  the  important  component 
of  the  residual  flow  to  the  body  is  symmetric.  As  the  body  becomes  flatter,  the 
residual  flow  may  lose  some  accuracy,  as  already  illustrated  in  Figure  10,  but 
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the  ability  to  furnish  an  end-plating  effect  remains  intact  even  for  very 
small  b/a  ratios.  Of  course,  if  b/a  becomes  small  enough,  the  body  can 
be  replaced  by  a  lifting  surface. 

3.1.6  The  Effect  of  Pickup  Point  Location 

The  averaging  process  (applied  to  the  residual  onset  flow)  outlined  in 

Section  2.5.6  involves  finding  the  flow  normal  to  the  body  surface  at  various 

points  around  the  cross  section.  These  points  are  input  quantities  to  the 

computer  program  of  the  Present  Method.  Figure  16  illustrates  the  effect  of 

increasing  the  number  of  these  points.  In  the  lower  curve,  pz/|j*  represents 

the  2- component  of  the  interference  doublet  strength  normalised  by  the  doublet 

strength  that  would  exist  in  the  fuselage  if  it  were  pitched  up  at  the  same 

angle-of-attack  as  the  wing.  In  this  case  the  wing  alone  is  pitched  up  to 

one  degree  and  the  value  of  u*  is  0.001095.  The  upper  curve  gives  the  value 

of  M  /„*.  Of  course,  for  this  symmetric  case,  this  quantity  should  be  zero. 

If  the  values  of  were  chosen  symmetrically  (say  e  =  45°,  135,  225,  315) 

then  the  ficticious  / ,*  would  not  appear.  As  the  values  of  o-  are 

y  v 

rotated  away  from  the  symmetric  position,  the  ficticious  uy/u*  appears.  At 
15°  of  rotation  it  is  at  its  worst  (upper  curve  given  by  circles).  At  this 
worst  possible  point  the  number  of  points  is  doubled.  The  result  is  a  two- 
thirds  reduction  of  the  error.  Thus,  as  the  number  of  integration  points 
increases,  the  accuracy  increases,  iiowever,  a  judicious  choice  cf  the,e 
points  will  allow  fewer  of  them  to  be  considered  for  the  same  accuracy. 

3.2  Correlation  of  the  Present  Method  with  Other  Methods  and  Experimental  Data 

3.2.1  Wing-Fuselage  Combinations  (Emphasis  on  Wing  Loads) 

A  wing-fuselage  comparison  presented  in  Reference  8  is  represented  in 
Figure  17.  In  addition,  calculations  using  the  Present  Method  are  shown. 

The  span  load  for  a  wing  (at  4.7°  angle-of-attack)  attached  to  a  circular 
fuselage  (at  0.7°  angle-of-attack)  is  given  in  this  figure.  The  agreement 
between  the  Present  Method  and  che  method  of  Reference  8  is  excellent  as  is 
the  correlation  with  the  experimental  data50.  Since  this  is  a  steady  case, 
the  agreement  between  the  two  methods  should  be  perfect;  however,  a  change  in 
the  numoer  of  spanwise  strips  near  the  tip  has  caused  a  slight  disparity. 
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A  second  wing-fuselage  comparison  presented  originally  in  Reference  18 
and  in  Part  I  of  this  report,  is  reproduced  in  Figure  18.  Shown  in  this 
figure  are  span  loads  for  two  modes  of  motion; 


Body  Pitch,  h2 
Body  Camber,  hz 


e 


iwt 


The  first  of  these  is  pitching  (of  the  body  only)  about  the  point  x  which  is 
tie  coordinate  of  the  center  of  the  booy.  The  second  mode  is  a  cambering  of 
tie  body  only.  The  nodal  point  is  at  xc.  The  maximum  amplitude  of  oscillation 
occurs  at  the  body  leading  and  trailing  edges  and  is  equal  to  the  maximum  body 
radius  at  its  center.  The  body  length  is  L. 


The  original  version  of  Figure  18  contained  a  comparison  of  three  differ¬ 
ent  methods.  To  this  list  the  Present  Method  is  added.  Again  the  Present 
Method  agrees  with  the  method  of  Reference  8.  (Actually,  for  this  figure, 
the  method  of  Reference  8  has  been  modified  to  account  for  the  doublet  distribu¬ 
tion  on  a  body  of  varying  radius.)  The  other  methods  referred  to  are  those 
of  Woodward  19  and  the  method  of  Part  I  of  this  report.  There  are  no  volume 
effects  for  this  problem  since  axial  sources  will  not  affect  the  wing  as  it  is 
placed  on  the  fuselage  for  this  case.  Actually,  the  nonlifting  (source)  effects 
of  volume  are  not  large  in  most  cases  and  „an  usually  be  ignored.  Also,  in 
unsteady  flow  there  are  ro  non  lifting  (source)  volume  effects  since  the  volume  is 

not  changing  with  time.  Currently,  no  volume  effects  are  contained  in  the  Present 
Method . 


The  results  presented  in  Figure  18  are  for  steady  flow.  The  method 
of  Reference  8  and  Woodward's  method  are  restricted  to  steady  flow;  however, 
the  method  of  Part  I  is  not.  Figure  IP  presents  a  comparison  of  the  Present 
Method  with  the  method  of  Part  I  for  the  same  configuration  considered  in 
Figure  18  but  oscillating  at  a  frequency  of  1.0.  This  figure  is  a  duplication 
of  one  found  in  Reference  18  and  Part  I  of  this  report,  with  the  results  of 
the  Present  Method  added.  The  real  part  of  the  span  load  due  to  body  camber 
was  plotted  with  the  wrong  sign  in  the  original  figures  but  is  corrected  here. 
Also,  in  Reference  18  there  appears  an  error  in  the  labeling  cf  the  ordinate. 
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The  agreement  between  the  two  results  is  good  for  the  body  pitch  mode 
but  is  only  fair  for  tne  body  camber  mode.  The  results  oc  the  Present  Method 
look  more  realistic,  especially  for  the  imaginary  part  of  the  span  loads  due 
to  body  camber.  The  probable  reason  for  the  disparity  is  the  idealization  of 
the  body  for  the  method  of  Part  I.  If  more  panels  were  placed  on  the  body, 
an  increase  in  accuracy  would  probably  result  (refer  to  Figure  11). 

3.2.2  Wing-Fuselage  Combination  (Emphasis  on  Fuselage  Loads) 

The  four  methods  compared  in  Figure  18  are  again  compared  in  Figure  20. 

The  original  version  of  this  figure  appeared  in  Reference  13.  To  this  fig¬ 
ure  the  results  of  the  Present  Method  have  been  added.  The  configuration 
consists  of  a  simple  swept  wing  attached  to  a  very  large  diameter  circular 
fuselage.  The  wing  is  at  l.G  radian  angle  of  attack  while  the  fuselage  is 
held  to  zero  incidence.  All  methods  are  in  good  agreement.  However,  the 
Present  Method  is  in  better  agreement  with  the  method  of  Reference  8  while 
the  Woodward  method  is  in  better  agreement  with  the  method  of  Part  I.  Again, 
the  method  of  Part  I  predicts  a  lift  coefficient  distribution  that  lies  below 
that  of  the  Present  Method. 

Special  attention  should  be  paid  to  the  fuselage  span  load.  The  usual  output 
of  the  Present  Method  is  a  longitudinal  loading  along  the  fuselage  length  (see 
Figure  21).  Also  the  total  lift  is  known  for  each  body.  If  this  total  lift 
is  normalized  by  the  root-chord  length  and  the  span  over  v.'hich  it  acts,  then 
an  average  span  load  is  the  result.  This  average  is  shown  as  a  horizontal 
dashed  line  across  the  fuselage.  The  area  under  the  horizontal  dashed  line 
seems  tc  be  low  compared  with  the  areas  produced  by  the  other  methods;  however, 
this  is  deceptive  and  requires  an  explanation. 

For  steady  flow  and  circular  fuselages,  it  can  be  shown  that  an  excellent 
approximation  to  the  fuselage  span  load  is  made  up  of  a  constant  load  minus 
an  elliptic  one.  The  constant  has  the  value  of  the  lift  coefficient  at  the 
root  and  this  constant  value  extends  from  the  fuselage  center  to  the  wing  root 
and  no  further.  The  lift  loss  due  to  the  elliptic  distribution  extends  from 
the  fuselage  center  to  the  edge  of  the  fuselage.  As  in  this  case  (Figure  20), 
the  fuselage  edge  extends  out  beyond  the  wing  root  position  and  thus  there  is 
a  region  between  the  wing  root  and  fuselage  edge  where  there  exists  only  the 
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negative  elliptic  distribution.  Combining  these  two  loadings  gives  the  dashed 
curve  across  the  fuselage  in  Fee  ire  .'0.  Notice  that  the  two  curves  combine  to 
give  a  posi>  /e  span  load  across  that  part  of  the  fuselage  that  does  not 
overlap  the  wing.  Over  the  overlapping  part,  the  lift  is  negative.  The  areas 
under  the  two  dashed  curves  (one  constant)  are  equal.  The  agreement  between 
the  method  of  Reference  8  and  the  Present  Method  is  excellent  over  that  part 
of  the  fuselage  not  overlapping  the  wing.  It  is  conventional  to  take  the 
negative  lift  distribution  that  exists  in  the  overlapping  region  and  add  it  to 
the  wing  span  load.  If  this  is  done,  all  curves  are  in  close  agreemert.  The 
deceptive  aspect  of  this  plot  is  the  fact  that  the  span  load  in  the  overlapping 
region  of  the  fuselage  and  wing  is  attributed  to  the  fuselage  alone.  The 
fact  is,  however,  that  the  scan  load  on  the  fuselage  in  this  region  is  actually 
negative.  This  explains  why  the  area  under  the  dashed  horizontal  curve  (which 
applies  only  to  the  fuselage)  is  equal  to  that  under  the  curved  dashed  line. 

The  spanwise  distribution  of  fuselage  lead  is  given  in  Figure  20. 

Figure  21  presents  the  longitudinal  variation  of  the  same  fuselage  lift. 

Four  curves  are  shown:  one  due  to  the  method  of  Part  I  and  three  due  to  the 
Present  Method.  The  loading  on  the  fuselage  is  due  entirely  to  the  wing 
(lift  carryover)  since  the  fuselage  is  at  zero  incidence.  expected,  the 
Present  Method  Alternate  II  does  not  compare  very  well  since  the  body  has  a 
large  diameter  compared  to  the  wing  dimensions.  As  noted  in  Figure  20,  how¬ 
ever,  the  total  lift  is  in  excellent  agreement  with  the  other  methods.  The 
moment  (not  shown)  is  also  very  good;  however,  the  distribution  of  force  is 
skewed  due  to  the  assumptions  made  in  alternate  II.  Some  of  this  skewness 
can  be  eliminated  if  the  approximate  redistribution  technique,  described  in 
Section  2.6.4,  is  employed.  The  curve  marked  Present  Method  Alternate  II 
(redistributed)  is  in  closer  agreement  with  the  method  of  Part  I.  Since  the 
body  is  circular,  Alternate  I  may  also  be  used.  The  curve  marked  Present 
Method  Alternate  I  shows  excellent  agreement  with  the  results  of  Part  I.  The 
results  of  the  Present  Method  Alternate  I  lie  slightly  above  the  results  of 
Part  I  but  this  is  expected  since  the  wing  loading  shows  the  same  trend. 

Figures  22  through  32  present  a  correlation  of  the  Present  Method  with 
the  experimental  data  of  Kornersl  and  the  theory  of  Labrujere20  et  al .  There 
are  a  total  of  four  configurations  considered:  1)  a  swept  wing  alone  (see 
Figure  23),  2)  a  swept  wing  and  a  fuselage  (D/c  =  1.0)  (see  Figure  24), 
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3)  a  swept  wing  and  a  fuselage  (D/c  =  0.68)  (see  Figure  26),  4)  and  a  strain!  t 
wing-fuselage  combination  (D/c  =  1.0)  (see  Figure  30).  The  swept  wingsin  as', 
of  the  configurations  are  similar.  Each  has  a  constant  chord  and  is  swept  45  . 
The  span  is  the  same  for  all  configurations  and  the  only  changes  to  the 
wings  are  those  necessary  lo  extend  the  wing  root  to  the  fuselage  surface. 

In  the  case  of  the  wing  alone,  the  wing  is  extended  to  the  y  =  0  plane. 

Figure  22  presents  the  idealization  of  configuration  2,  i.e.,  swept  wing  with 
D/c  =1.0.  This  idealization  is  typical  of  all  the  others.  The  solid  lines 
represent  lifting  surface  or  body  interference  element  boundaries.  The  dashed 
lines  represent  slender  body  element  boundaries. 

The  experimental  data  of  Reference  51  are  given  for  various  angles  of 
attack.  The  value  of  C£a  (per  degree)  may  be  obtained  from  this  data. 

Figure  23  presents  a  comparison  of  the  distribution  of  lift-cu^ve-slope  as 
calculated  for  a  wing  a  one  using  the  Present  Method  and  as  obtained  from 
experimental  data.  Experimentally,  there  is  an  obvious  loss  of  lift  due  to 

viscosity.  !otioe  that  c,  is  lower  at  a  =  6°  than  at  3°.  This  is 

*-n 

a  clear  indication  of  a  nonlinear  viscous  effect  in  the  data. 

Figure  24  presents  a  similar  comparison  except  that  the  configuration 
consists  of  a  wing  and  fuselage.  Trie  fuselage  is  at  zero  incidence  while  the 
wing  is  pitched  up.  Again  c  across  both  the  wing  and  fuselage  shows 

the  characteristic  loss  cf  lift  due  to  viscosity. 

This  loss  of  lift  due  to  viscosity  is>  accentuated  when  the  fuselage  is 
given  the  same  angle  of  attack  as  the  wing.  Figure  25  presents  such  a 
comparison 

Figures  26  and  27  show  comparisons  similar  to  those  given  in  Figures  24 
and  25.  The  configuration  is  almost  identical  except  for  the  fact  that  the 
fuselage  is  smaller  (D/c  =  0.68).  The  span  load  behaves  in  an  interesting 
manner  in  Figure  27.  The  comparison  seems  to  be  very  good  over  the  outboard 
naif  semispan.  The  agreement  then  deteriorates  toward  the  wing-fuselage 
intersection.  This  could  be  due  to  the  thickening  of  the  boundary  layer  in 
the  region  of  the  wing-fuselage  intersection. 

The  longitudinal  loading,  as  calculated  by  the  Present  Method,  is 
compared  to  that  dettrmined  experimentally  in  Figures  28  and  29.  The  fuselage 


incidence  is  ze"o  in  Figure  28  and  equal  to  that  of  the  wing  in  Figure-  r: . 

In  the  region  near  the  wing  Alternate  I  and  the  redistributed  Alternate  II 
are  in  good  agreement  with  the  experimental  data.  A  large  disagreement 
between  Alternate  I  arid  Alternate  II  appears  at  the  leading  and  trailing 
edge"  of  the  body  in  Figure  29.  The  basic  reason  for  this  is  the  fact  that 
the  '  ;tal  lift  of  a  slender  body  element  is  lumped  at  that  element  in  the 
method  f  Alwmate  II  where  as  it  is  cistrioutea  over  tne  body  for  Alternate 

I  and  Alternate  II  ( ^distributed) . 

Labrujere  et  a*.  0  compared  their  steady  flow  theory  with  Korner's  data 
for  an  unswept -wing/L  .ombi nation,  figures  30  and  31  present  comparisons 
of  the  spanwise  lift  coe  n'cient  distribution  for  a  wing  at  6°  angle  of  attack. 
The  fuselage  is  at  zero  degrees  incidence  in  Figure  30  and  at  6°  in  Figure  31. 

For  the  case  of  zero  fuselage  incidence  the  i-Aesent  Method  is  in  closer 
agreement  with  the  method  of  Labrujere  <  .  al.  than  with  the  experimental  data. 
However,  for  the  case  where  the  wing  anc  i.selage  are  at  the  same  incidence, 
the  distribution  claculatei  by  the  Preser  *.  Method  lies  Imost  equidistant 
between  the  datd  and  the  method  of  Labruji'--;  et  al.  The  discontinuity  in 
the  span  load  apparent  in  the  theory  of  Ladrijeie  et  al .  s  caused  by  the  fact 
that  the  fuselage  does  not  close,  due  to  t.e  v'.Ti  tunnel  s  ng,  and  thus  there 
exists  a  resultant  slender  body  lift. 

Figure  32  presents  a  comparison  of  experimen  al  and  calcinated  longitudinal 
loading  for  the  zero  incidence  fuselage  case.  The  , "distribute  loading  (Alter¬ 
nate  II)  is  in  very  good  agreement  with  the  experiment*  data  as  is  Alternate  I. 
The  original  distribution,  also  now  shown,  however,  sh.  ws  large  variations.  The 
straight  wing  case  exaggerates  these  variations  because  i  »e  contribution  of 
every  wing  lifting  surface  element  occurs  on  the  body  surf  a  -e  in  that  sho-"fc 
space  lying  between  the  root  chord  leading  and  trailing  edges.  There  is  a 
negative  dp'\  •  the  longitudinal  loading  which  is  due  to  the  fuselage  ir  re¬ 
ference  eleiii'v.  These  interference  elements  react  to  the  upwash  generated 
at  the  body  g*.c.  to  the  bound  vortices  on  the  wing. 

The  last  two  correlations  of  this  section  are  found  in  Figures  33  and  3^. 

The  experimental  and  calculated  longitudinal  loading  on  a  transport  aircraft 
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are  compared  for  two  values  of  Mach  number.  In  both  cases  ‘.he  redistributed 
loadings  are  used  for  alternate  II.  The  agreement  is  good  especially 
for  the  M  =  0.6  case.  At  the  higher  Mach  number  it  may  be  that  shock  waves 
have  interf erred  with  the  flow  at  the  point  near  the  wing  root  trailing  edge, 
where  the  experimental  data  takes  a  sharp  jump. 

3.3  T-Tail  Fuselage  Combination 

Zwadn52»S3  has  developed  a  kernel  function  technique  for  T-tails.  Of 
particular  interest  are  the  span  load  plots  of  Reference  52  for  a  T-tail 
tested  with  ooth  a  ground  plane  and  a  fuselage.  Figure  35  presents  a  compari¬ 
son  of  spanwise  loading  as  calculated  by  the  Present  Method,  by  Zv/aan  and  as 
determined  experimental iy  for  a  T-tail  in  yaw.  Results  for  the  Present  Method 
are  given  both  for  the  T-tail  with  ground  plane  and  fuselage.  The  experimental 
data  are  given  only  for  the  T-tail  with  ground  plane.  Difficulty  was  experienced 
ir  interpreting  Zwaan's  results.  In  Reference  52  a  symbol  |  k  |  was 
gi/en  as  the  absolute  value  of  the  loading.  This  symbol  was  not  defined  in 
this  r-ference.  Integration  of  experimental  pressure  data  indicates  that  | k j 
is  (cc;)/(cCy),  where  Cy  is  based  on  the  total  area  of  the  vertical  fin. 

Using  this  definition,  the  correlation  of  the  data  with  the  Present  Method 
for  the  vertical  stabilizer  is  good.  (Refer,  to  the  calculation  of  the  Present 
Method  marked  "Ground  Effect.")  Both  the  absolute  magnitude  and  phase  angle 
compare  well.  If  the  same  definition  of  |kj  holds  for  the  horizontal  stabil¬ 
izer,  ther  there  is  an  inconsistency  in  Zwaan's  calculation  as  oresented  in 
Reference  52,  Hgu^e  8.  The  jump  in  jkj  on  the  vertical  fin  at  the  intersec¬ 
tion  of  the  hoiizontal  stabilizer  should  be  twice  the  value  of  k  at  this 
intersection  on  the  horizontal  stabilizer.  If  the  area  of  the  horizontal 
stabilizer  is  used  in  the  definition  of  Cy  then  Zwaan’s  calculated  results 
become  consistent.  (The  value  of  |k|  on  the  horizontal  at  the  root  then 
becomes  approximately  0.415  while  half  the  jump  on  the  vertical  is  about 
C.425.)  The  new  definition  for  both  Zwaan's  calculation  and  the  experimental 
data  are  used  to  give  the  results  of  Figure  35  for  the  stabilizer.  The  agree¬ 
ment  is  good. 

Tne  effects  of  using  a  ground  plane  in  place  of  the  fuselage  are  shown  in 
Figure  35.  The  two  calculations,  done  using  the  Present  Method,  show  that  the 
g"ound  plane  increases  the  loading  substantially,  both  on  the  vertical  fin  and 
horizontal  stabilizer. 
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Experimental  data  have  also  been  obtained  for  the  T-tail  and  fuselage  in 
steady  flow  (yaw).  The  shape  of  the  fuselage  was  taken  from  a  small  sketch 
given  in  Reference  52,  Figure  12.  This  sketch  is  reproduced  in  Figure  35. 

Figure  36  shows  a  similar  sketch  with  the  addition  of  the  idealized  fuselage 
that  is  used  with  the  Present  Method.  The  results  correlate  well  except 
near  the  fuselage. 

3.4  Wing  and  Tip  Mounted  Nacelle  Combination 

The  only  experimental  data  that  have  been  found  for  a  wing/body  configu¬ 
ration  in  oscillatory  motion  is  that  reported  by  Cazemier  and  Bergh54 >55»56. 

A  wooden  model  of  a  low-aspect-ratio  wing,  fitted  with  a  large  wing-tip 
nacelle  was  tested  over  a  wide  range  of  frequencies.  The  wind  tunnel  Reynolds 
number  was  low,  Re  =  4.7  x  106/meter,  and  the  model  was  about  one  meter  in 
semi  span.  The  test  was  set  up  to  read  the  lifting  pressure 
[Cp  (lower)  —  Cp  (upper)]  directly.  The  pressures  were  transmitted  from 
the  model  through  calibrated  tubes.  The  tubes  were  calibrated  to  eliminate 
the  phase  shift  caused  by  the  transmission  time  required  to  communicate  the 
pressure  from  the  model  to  the  pick-up  point. 

An  extensive  correlation  of  the  Present  Method  with  this  experimental 
data  is  presented  -  since  it  is  the  only  data  of  its  kind  available.  Figures 
38  through  57  are  devoted  to  this  correlation. 

Figure  37  illustrates  the  idealization  of  the  two  configurations  con¬ 
sidered,  a  wing  alone  and  a  wing/nacelle.  The  wings  in  the  two  configurations 
are  slightly  different.  The  first  wing  strip,  in  the  wing/nacelle  configuration, 
is  missing  on  the  wing-alone  configuration.  The  last  strip  on  the  aileron  is 
filled  in  for  the  wi  mQ  d 1 Cud  configuration.  The  interference  elements  are 
shown  on  the  nacelle,-  however,  the  slender  body  elements  are  not.  On  the  wing 
the  inboard  control  surface  is  called  the  flap  and  the  outboard  control  surface 
is  called  the  aileron.  The  pitch  and  roll  axes  are  presented  in  Fi.iure  37, 
as  are  the  stations  at  which  experimental  pressures  were  taken. 

A  comparison  of  ’ifting  pressure  for  the  wing  alone  in  steady  pitch  is 
given  in  Figure  38.  The  experimental  data  lies  above  the  calculated  results 
for  all  stations  presented.  This  is  very  extraordinary  since  experimental 
data  usually  lies  below  the  calculated  results.  Normally  the  wing  will  be 
less  effective  as  a  lifting  surface  due  to  the  action  of  viscosity.  The 
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reason  for  this  peculiar  behavior  is  not  known.  A  similar  comparison  for  a 
deflected  flap  is  shown  in  I  fgure  39.  The  correlation  is  good  and  shows  the 
correct  relationship  between  calculated  results  and  experimental  data.  A 
further  comparison  for  a  deflected  aileron  is  given  in  Figure  40.  The  correla¬ 
tion  is  not  quite  as  good  as  for  the  flap. 

In  all  of  the  comparisons  the  pressures  are  normalized  as  follows: 

ACn  =  AP/qS 

Is  =  1.097m,  plunge 

a  ,  pitch 

sf  ,  flap 

6_  ,  aileron 

a 

<t>  ,  roll 

where  a,  6^,  6a,  and  <£  are  the  angle-of-attack,  flap  am  1c ,  aileron  angle 
and  roll  angle.  The  reduced  frequency  is  based  on  c  =  l.OS/r,  fcr  both  the 
wing  alone  and  wing/nacelle  cases.  The  reduced  frequenc  as  delined  in 
References  54,  55  and  56  is  uic/U^,  which  is  twice  s  large  es  that  defined 
in  this  report,  wF/21^.  Thus,  reduced  frequencies  of  0.5  and  1.0  listed  in 
the  figures  of  this  report  represent  frequencies  of  1.0  and  2.0  ir  the  nota¬ 
tion  of  References  54,  55  and  56. 

A  comparison  of  lifting  pressures  for  the  wing  a 'one  oscillating  at  a 
reduced  frequency  of  0.5  is  given  in  Figure  41.  The  r.-al  part  of  the  lifting 
pressure  for  pitching  motion  is  very  similar  to  the  steady  case.  The  correla¬ 
tion  at  this  reduced  frequency  leaves  something  to  be  desired,  especially  for 
the  plunging  case.  Pressure  plots  for  rolling  and  aileron  motion  are  given 
in  Figure  42.  The  case  of  rolling  is  very  similar  to  that  of  plunging,  as  it 
should  be.  The  aileron  correlation  is  good.  Figure  43  is  the  same  as  41 
except  the  frequency  is  1.0  instead  of  0.5.  Surprisingly,  the  agreement 
between  data  and  calculated  results  is  better  at  this  reduced  frequency  (1.0) 
than  it  is  at  0.5.  The  same  comment  can  be  made  of  the  plunge  case  shown  in 
the  same  figure.  An  extra  station  is  included  for  the  pitching  case  to 
compare  Laschka's  method57  with  the  Present  Method  and  the  experimental  data. 
Laschka's  theory  (as  read  from  a  very  small  figure)  is  given  for  strip  number 
5.  Essentially,  the  results  by  Laschka  are  the  same  as  the  Present  Method. 
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A  comparison  of  pressures  for  rolling  and  flapping  is  presented  in  Figure  44. 
Again  the  correlation  is  somewhat  better  for  the  reduced  frequency  of  1.0 
than  it  is  f,.r  0.5. 

The  above  comparisons  refer  to  the  wing-alone  calculation.  These  were 
presented  so  that  the  nacelle  effects  present  in  the  n^xt  set  of  comparisons 
could  be  properly  assessed. 

The  lifting  pressure  distribution  for  the  wing/nacelle  at  station  2  is 
presented  in  Figure  45;  both  calculated  results  and  experimental  data  are 
shown.  The  comparison  is  as  good  as  the  wing-alone  case.  A  comparison  of 
lifting  pressure  for  pitch,  plunge  and  roll  at  a  reduced  frequency  of  0.5  is 
given  in  Figure  46.  Generally,  the  agreement  between  the  Present  Method  arid 
experimental  data  is  better  than  for  the  wing-alone  case.  As  much  cannot  be 
said  of  the  agreement  at  the  reduced  frequency  of  1.0.  Figure  47  presents 
this  comparison.  The  agreement  between  calculated  results  and  experimental 
data  for  pitch  is  good,  while  that  for  plunge  and  roll  is  not  as  good.  In  all 
cf  the  results  where  the  agreement  is  not  as  good  as  it  might  be,  there  seems 
to  be  a  phase  shift  of  the  results. 

The  pressure  plots  just  discussed  are  presented  to  compare  experimental 
and  calculated  chordwise  loadings.  In  many  instances,  a  comparison  of  these 
two  for  the  span  load  or  lift  coefficient  distribution  is  more  enlightening. 

For  this  reason,  the  experimental  lifting  pressure  data  avdilaLle  in  References 
55  and  56  has  been  integrated  to  give  lift  coefficient.  The  method  of  integra¬ 
tion  was  to  connect  the  experimental  data  by  straight  lines  and  integrate  by 
pi  ammeter . 

The  fir.,-,  of  the  lift  coefficient  distribution  comparisons  is  given  in 
Figure  48  for  the  steady  case.  The  upper  curve  gives  a  comparison  of  lift 
coefficient,  as  calculated  by  the  Present  Method,  and  as  determined  experi¬ 
mentally  for  the  wing  alone  case.  The  agreement  is  generally  oood  except  for 
the  strange  property  that  the  experimental  data  lies  above  the  calculated 
results.  The  lower  curve  presents  a  similar  comparison  for  the  wing/nacelle 
case.  There  are  two  curves  marked  "Present  Method."  The  lower  curve  is  the 
result  of  using  the  u  .ual  method  where  axial  elements  in  conjunction  with 
images  are  used  to  model  the  flow.  The  upper  curve  results  were  obtained 
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after  lifting  surface  elements  were  placed  on  the  body  surface  (in  conjunction 
with  images).  This  upper  curve  is  in  better  agreement  with  the  data.  The 
exults  of  the  upper  curve  are  obtained  at  the  cost  of  doubling  the  number 
of  unknowns  in  the  problem.  Thus,  accuracy  and  efficiency  must  be  balanced 
by  the  user  of  the  Present  Method.  Also  shown  on  the  wing/nacelle  (lower)  plot 
is  trie  wing-alone  case.  There  is  a  very  large  difference  between  the  wing- 
alone  and  wing/nacelle  lift  coefficient  distributions.  The  Present  Method 
goes  a  long  way  toward  predicting  the  full  effect  of  the  nacelle  even  when 
lifting  surface  elements  are  not  placed  on  the  body  surface. 

At  a  reduced  frequency  of  0.5,  the  comparison  between  calculated  results 
and  experimental  data  shows  good  agreement  for  the  real  part  and  poor  agreement 
for  the  imaginary  part.  Figure  49  presents  such  a  comparison  for  plunging 
and  pitching.  A  comparison  for  aileron  motion  is  given  in  Figure  50.  This 
figure  shows  that  the  aileron  is  not  as  effective  predicted,  espeically  for 
the  imaginary  part.  The  comparison  of  lift  coefficient  distribution  for  a 
reduced  frequency  of  1.0  is  given  in  Figure  51  for  roll  and  pitch.  The  agree¬ 
ment  between  calculated  and  experimental  data  is  worse  for  the  real  part  and 
better  for  the  imaginary  part  when  compared  with  the  results  at  the  reduced 
frequency  of  0.5.  The  results  presented  in  Figure  52  are  of  interest  when 
compared  with  those  of  Figure  50.  In  Figure  50  the  experimental  data  lay 
below  the  calculated  results  for  both  real  and  imaginary  parts.  In  Figure  51 
the  real  pa^t  of  the  data  lies  above  the  calculated  results  while  the 
imaginary  part  lies  below. 

The  next  set  of  comparisons  is  for  the  lift  coefficient  distribution  for 
the  wing/nacelle  combination.  One  peculiarity  of  these  plots  is  that  there 
are  discontinuities  in  the  distributions  at  the  root  and  tip  of  the  wing. 

These  are  caused  by  a  discontinuous  chord  length  at  these  locations.  Refer¬ 
ence  back  to  Figure  37  will  show  immediately  these  planform  peculiarities. 

The  agreement  between  calculated  and  experimental  spanwise  distribution 
of  lift  coefficient  for  the  wing/nacelle  combination  is  about  the  same  as  that 
for  the  wing  alone.  Figure  53  presents  these  results  for  plunge  and  pitch 
motions.  Except  for  some  scatter  (rtation  2)  the  real  parts  are  in  close  agree¬ 
ment  while  the  experimental  data  lies  below  the  calculated  results  for  the 
imaginary  part.  The  comparison  at  a  reduced  frequency  of  1.0  is  shown  in 
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Figure  54  for  rolling  and  pitching.  Again  the  comparison  is  about  the  same  as 
the  wing-aione  case. 

In  all  of  the  comparisons  presented  for  the  wing-alone  and  wing/nacelle 
there  seems  to  be  a  phase  shift  between  the  calculated  results  and  the 
experimental  data.  The  phase  shift  seems  to  increase  as  the  frequency 
increases.  The  reason  for  this  is  not  known.  It  is  felt  that  the  agreement 
should  be  better,  espeically  for  the  wing-alone  calculation,  since  similar 
comparisons  have  been  made  in  References  10,  12  and  13  and  show  good  agreement. 

Further  analysis  of  the  lift  coefficient  data  has  been  done  in  an 
attempt  to  isolate  the  nacelle  effect  on  the  wing.  Specifically,  the  wing- 
alone  distributions  have  been  subtracted  from  the  wing/nacelle  distributions 
both  for  the  calculated  results  and  the  experimental  data.  Some  irregularities 
are  introduced  when  this  is  done  since  the  wing  without  the  nacelle  is  not 
exactly  the  same  as  that  with  the  nacelle.  Some  of  the  effects,  then,  will  be 
due  to  this  difference  and  not  due  to  the  nacelle. 

The  steady  flow  results  for  pitch  are  given  in  Figure  55.  The  experimental 
data  lies  generally  between  the  two  approaches  of  the  Present  Method  (one  with 
axial  body  elements,  one  with  surface  elements,  both  with  images).  The  same 
type  of  comparison  is  given  in  Figure  56  for  a  reduced  frequency  of  0.5.  The 
oata  shows  scatter  for  the  real  part  both  for  plunge  and  pitch  although  the 
agreement  is  excellent  near  the  wing/nacelle  intersection.  The  agreement  for 
the  imaginary  parts  for  pTunge  and  pitch  is  not  very  good.  When  the  reduced 
frequency  ia  increased  to  1.0,  the  correlation  for  the  real  part  improves, 
generally,  as  mown  in  Figure  57.  The  imaginary  part,  however,  is  not  as  good, 
especially  for  pitch.  Also  shown  in  Figure  57  are  results  obtained  from 
Part  I  of  this  report  (also,  Reference  18)  for  the  wing  in  pitch.  The  effect 
of  the  nacelle  in  these  results  is  very  large.  This  discrepancy  may  be  due  to 
the  fact  that  a  slightly  different  wing  planform  was  used  in  Part  I.  Specifically, 
the  gaps  at  the  wing  root  and  tip  are  not  present  in  the  calculation  presented 
in  Part  I.  Also,  the  nacelle  cross  section  was  represented  by  a  rectangle. 

This  rectangle  was  large  enough  to  envelope  the  nacelle.  Thus,  the  end-plating 
effect  was  larger. 
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3.5  Wing-Nacelle-Fuselage  Combination 

A  good  examp’e  of  a  complex  wing-nacell e-fuselage  combination  is  found 
in  the  B- 58  bomber.  Data  were  made  available  from  General  Dynamics , Fort  Worth 
Division,  through  Dr.  R.  G.  Bradley.  The  data  are  For  configurations 
with  and  without  nacelles.  A  detailed  description  of  the  configuration 
is  found  in  Reference  58.  The  idealization  of  this  configuration  for  calcula¬ 
tions  by  the  Present  Method  is  given  in  Figure  58.  The  fuselage  was  cambered 
as  well  as  the  wing.  The  wing  camber  was  restricted  to  a  conical  region  at 
the  leading  edge.  This  region  is  easily  identified  in  Figure  58  as  the  one 
containing  the  conical  rays  emanating  from  a  point  on  the  x-axis.  The  camber 
is  approximated  by  three  constant  slope  conical  wedges.  A  more  conventional 
idealization  is  developed  for  that  region  of  the  wing  without  camber  and  for 
the  pylon.  The  outboard  nylon  is  omitted  since  it  was  thought  to  be  too  small 
to  effect  the  calculations.  There  is  a  gap  between  the  wing  and  outboard 
nacelle.  The  outboard  nacelle  is  at  a  negative  incidence  and  therefore  the 
leading  edge  is  considerably  lower  than  the  trailing  edge  where  it  attaches 
to  the  wing.  The  average  gap  at  the  center  of  the  nacelle  is  the  one  shown 
in  the  figure.  The  small  arrows  on  the  fuselage  and  nacelle  axes  indicate 
the  directions  of  force  (and  doublet  orientation)  that  are  allowed  for  the 
calculation.  The  interference  elements  are  shown  on  the  fj„°lage  and  nacelles. 
Shown  on  the  fuselage  alone  are  the  slender  body  elements  indicated  by  tick 
marks  on  the  actual  fuselage  shape. 

A  comparison  of  the  Present  Method  with  experimental  data  for  the  configu¬ 
ration  with  nacelles  is  presented  in  Tigure  59.  Specifically,  tne  span  load  for 
the  wing  at  4  degrees  anqle-of-attack  (in  addition  to  the  conical  camber) 
presented.  The  experimental  data  lie  above  the  calculated  values.  Similar 
results  were  obtained  by  Bradley  and  Miller21  for  the  B-58  bomber.  The  cause 
of  tins  is  found  in  the  very  high  sweep  of  the  leading  edge.  As  is  well  known, 
wings  with  large  leading  edge  sweep  angles  develop  a  leading  edge  vortex  which 
causes  an  increase  in  the  lift.  Figure  60  presents  a  similar  comparison  for 
the  case  of  the  wing  without  the  nacelles.  Again  the  leading  edge  vortex 
increases  the  lift  over  the  calculated  value. 

To  determine  the  effect  of  the  nacelles  on  the  wing  span-loading,  the 
calculated  results  of  Figures  59  and  60  are  replotted  on  a  single  curve 
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in  Figure  61.  A  calculation  was  made  in  whicli  the  nacelle  diameters  were 
doubled.  The  results  of  this  calculation  are  also  shown  in  Figure  61.  The 
obvious  conclusion  that  is  drawn  from  this  figure  is  that  as  the  nacelle 
diameter  is  increased,  the  lift  on  the  wing  is  increased.  The  nacelles  's 
well  as  the  pylon  tend  to  prevent  spanwise  flow  and  in  this  way  increase 
the  lift. 

The  span  load  on  the  B-58  wing  in  oscillatory  motion  is  given  in  Figure  62. 
Specifically,  the  aircraft  was  made  to  pitch  about  the  wing  "apex"  (determined 
by  an  extension  of  the  wing  leading  edge  to  the  center  of  the  fuselage)  at  a 
reduced  frequency  of  0.5  { c  =  19.15  in.).  The  normalizing  chord  length  used 
for  the  span  load  is  cR  which  is  the  root-chord  length  (34.7  in.).  These 
dimensions  are  for  the  wind  tunr.el  model.  One  further  calculation  is  pre¬ 
sented  i»-  Figure  62.  The  aircraft  is  made  to  operate  in  ground  effect,  at  a 
height  of  10  in.  above  the  ground.  Tne  height/root  chord  ratio  is  0.288.  The 
effect  of  the  ground  is  appreciable  especially  for  the  real  part  of  the  span 
loading  inboard  of  the  pylon.  There  the  load  is  nearly  doubled.  Also,  nearly 
doubled  for  the  ground  effect  case  is  the  computing  time.  The  reason  for 
this  is  th3t  the  kernel  is  evaluated  twice  as  many  times  in  tne  ground  effect 
case.  Each  ^ending  point  on  the  aircraft  has  a  counterpart  in  its  ground 
effect  image. 

The  span  load  calculations  presented  in  Figure  62  are  for  unit  angle  of 
attack,  no  cambering  modes  are  present.  For  this  case  a  simple  idealization 
of  the  wing  is  used. 

The  fuselage  loading  associated  with  the  unsteady  case  (out  of  ground 
effect)  is  presented  in  Figure  63.  The  original  loading  for  Alternate  II  is 
presented.  '  this  case  the  original  loading  and  the  redistributed  loading 
should  be  very  close  because  of  the  slenderness  of  the  fuselage.  The  major- 
contribution  to  the  fuselage  lift  is  obtained  in  the  region  of  the  wing 
(indicated  in  Figure  63  by  a  heavy  line  lying  on  the  x-axis).  It  is  inter¬ 
esting  to  note  that  the  point  of  maximum  loading  for  the  imaginary  part  lies  a 
considerable  distance  aft  of  the  point  of  maximum  loading  for  the  real  part. 


99 


The  loading  on  tie  nacelle  is  considerably  less  thin  that  on  the  fuselage. 
Figure  64  shows  that  the  maximum  lift  (or  side  force)  coefficient  is  about 
6.0,  while  Figure  63  shows  the  maximum  for  the  fuselage  is  about  30.  Figure  64 
presents  the  vertical  loading  on  the  inboard  nacelle  and  the  horizontal  load¬ 
ing  on  the  outboard  nacelle.  The  load  distributions  given  in  Figure  64  have 
not  been  redistributed  and  are  the  original  ones  of  Alternate  II. 
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4.0  conclusion:  and  recommendations 


4.1  Conclusions 

A  method  for  predicting  oscillatory  loads  on  very  general  configura¬ 
tions  has  been  developed.  Configurations  may  include  a  combination  of  any 
or  all  of  the  following  components:  1)  lifting  surfaces  such  as  wings,  pylons, 
stabilizer,  fin,  etc.,  with  arbitrary  dihedral,  2)  partial  or  full  span 
control  surfaces  and,  3)  bodies  such  as  fuselages,  nacelles,  stores,  etc. 
with  elliptic  cross-sectional  shapes.  The  operating  conditions  are  also 
very  general:  1)  all  frequencies  of  practical  interest  and  all  subsonic 
Mach  numbers,  2)  symmetry  and  ground  ecfect,  3)  mutual  interference  of  lifting 
surfaces,  and  4)  multiple  modes  of  oscillation  (described  by  polynmials) . 

One  of  the  main  features  of  the  Present  Method  is  the  efficiency  wUh 
whicn  these  configurations  and  conditions  are  handled.  With  the  kernel 
function  techniques  even  simple  plane  wings  present  time  consuming 
numerical  difficulties.  If  it  were  possible  to  extend  the  kernel  function 
method  to  the  status  of  the  present  method  the  computational  effort  might 
be  excessively  large. 

The  effects  of  varying  some  of  the  important  parameters  irvolved  have 
been  studied.  Specifically  che  effects  considered  are  those  of:  s)  body 
radius,  2)  cross-sectional  aspect  ratio,  3)  frequency,  4)  wing  position 
on  the  fuselage,  and  o)  end  plating  effects  of  a  body.  The  conclusions 
drawn  from  these  studies  are  outlined  below: 

1)  The  body  radius  has  two  distinct  effects:  one  is  due  to  the 
flow  field  associated  with  a  body  at  angle  of  attack  and  the  other  is 
associated  with  interference.  The  fuselage  angle  of  attack  effect  increases 
the  wing  loading  as  the  ratio  of  fuselage  diameter  to  wing  span  is  increased. 
The  fuselage  interference  effect  shows  variations  as  the  radius  ranges 

from  zero  to  infinity  but  is  the  same  at  the  two  extremes. 

2)  The  interference  effect  of  fuselage  cross-sectirnal  aspect  ratio 

(ratio  of  b/a  where  b  =  semi -height,  a  =  semi -width)  is  monotonic, 
reaching  its  greatest  value  when  the  fuselage  is  a  vertical  slit  b/a  -  . 
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The  angle  of  attack  effect  is  just  the  opposite  reaching  a  minimum  when 
b/a  =  «>.  The  ‘image  system  for  elliptic  cross  sections  starts  to  fail 
when  b/a  is  greater  than  2.0  or  less  than  0.5. 

3)  The  effects  of  frequency  for  wing- body  combinations  are  very 
similar  to  the  wing  alone  case. 

4)  Placing  a  wing  high  on  a  fuselage  of  elliptic  cross  section 
tends  to  brinq  the  left  and  riqht  winq  halves  closer  together  This 
increases  thp  interference  effect  of  the  body.  The  angle  of  attack  effect, 
however,  is  reduced. 

5)  The  Present  Method  predicts  the  end  plating  effect  of  a  body  very 
accurately  for  a  circular  cross  section  (b/a  =  1).  However,  when  the 
ratio  of  b/a  is  near  the  value  of  2.0  (or  greater)  the  predicted  end  plating 
effe'~t  starts  to  fall  Delow  the  correct  value.  There  are  two  reasons  for 
this  behavior.  The  first  is  the  fact  that  the  image  within  the  ellipse  is 
terminated  at  its  centerline.  The  second  reason  is  that  the  interference 
elements  cannot  account  for  an  anti -symmetric  flow  across  the  body  cross 
section. 

If  conditions  or  configurations  are  simplified,  then  other  methods  may 
be  used  to  compare  with  the  present  method.  Cases  which  involve  only  lifting 
surfaces  may  be  handled  by  the  methods  of  Laschka57  and  Zwatt,^.  Comparisons 
of  these  methods  with  the  Present  Method  show  good  agreement.  For  the  steady 
case  the  theories  oi  Woodward1 ^  and  Labrujere20  may  be  used  for  comparison. 
Loads  predicted  by  the  Present  Method  fall  slightly  below  the  thick  winq  theory 
of  Labrujere  but  lie  slightly  above  those  of  Woodward. 

The  Present  Method  is  also  compared  with  the  method  of  Part  I  of 
this  report.  The  method  of  Part  I  employs  lifting  surface  elements  on 
the  body  surface  to  account  for  the  body  interference  effects.  The  results 
are  highly  dependent  on  the  idealisation.  If  the  cross  section  is 
represented  by  orly  a  few  defining  elements  then  the  results  can  be 
inaccurate  especially  for  unsteady  flow. 

Several  general  conclusions  can  be  made  about  the  comparison  of  the 
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Present  Method  and  the  method  of  Part  I.  In  steady  flow  the  loading,  as 
calcilated  by  the  method  of  Part  I,  usually  falls  s'':htiy  below  that 
calculated  by  the  Present  Method.  As  the  number  of  body  surface  elements 
is  increased  (for  the  method  of  Part  I)  the  agreement  between  the  two 
methods  improves.  For  unsteady  flow  the  situation  is  more  complicated. 

There  exist  differences  between  the  two  methods  especially  near  the  wing- 
fuselage  intersection.  If  the  representation,  in  the  method  of  Part  I, 
is  crude  then  variations  in  the  results  occur  that  do  not  exist  in  the 
results  of  the  Present  Method.  In  some  instances  a  third  calculation  was 
done  in  which  both  images  and  lifting  surface  elements  were  used.  This 
calculation  confirmed  the  results  of  the  Present  Method  in  most  cases. 

Tne  exceptions  to  this  were  cases  involving  bodies  placed  at  wing  tips. 

The  total  loads  and  moments  on  bodies,  as  calculated  by  the  Present 
Method,  are  very  accurate.  The  spariwise  variation  of  lift  coefficient  agrees 
very  well  with  the  method  of  Part  I,  Woodward's  method  and  tne  method  of 
Reference  8.  The  longitudinal  distribution  of  load  along  the  body  surface 
is  not  as  accurate  for  Alternate  II.  The  reason  for  this  is  that  the  total 
load  on  the  body  due  to  one  lifting  surface  element  is  found  and  lumped  at 
the  same  longitudinal  location  as  the  lifting  surface  elemenc  for  this  method. 
The  longitudinal  distribution  on  the  body  arises  then  only  from  the  longi¬ 
tudinal  distribution  of  lifting  surface  elements.  A  redistribution  of  the 
load,  based  on  an  approximate  procedure,  increases  the  accuracy  under  most 
conditions . 


iiany  c^mpa  i  ipwiu  vs  v'»c  i  cou  i  to  wi  me  r  i  ciciib  i'ic  u.uu  wi  fi 

experimental  data  are  presented  in  this  report.  Some  of  the  results  are 
apparently  contradictory .  For  this  reason  the  results  will  be  discussed 
for  each  of  the  experiments. 
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The  data  of  Martina  for  a  wing-fuselage  combination  agrees  very 

well  with  the  results  of  the  Present  Method.  Similar  data  were  obtained 
51 

by  K<5rner  for  low  Reynolds  Number  and  lew  Mach  Number.  The  data  show 
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a  characteristic  loss  of  effectiveness  relative  to  the  calculated  results 
which  is  the  direct  result  of  the  boundary  layer.  Specifically  tne  boundary 
layer  is  much  thicker  on  the  upper  surface  and  effectively  produces  an 
"uncanbering"  of  the  wing.  The  boundary  layer  is  also  very  thick  at  the 
wing-fuselage  intersection  and  thus  it  is  not  surprising  that  the  loss  of 
lift  is  found  to  be  greatest  there.  The  redistributed  longitudinal  loading 
or.  the  fuselage  is  in  good  agreement  with  the  experimental  data. 

The  only  experimental  data  found  for  a  wing/body  combination  in 

54  55  56 

oscilatory  flow  were  reported  by  Cazemier  and  Bergh  ’  '  ’  who  considered 
a  wing  with  a  tip-mounted  engine  nacelle.  These  data  were  obtained  at  low 
Reynolds  Number  and  low  Mach  Number.  For  reference  purposes  some  of  the 
data  were  obtained  for  the  wing  alone  case.  The  experiment  for  the  wing 
alone  in  steady  flow  was  found  to  lie  slightly  above  the  calculated  results. 
This  is  a  very  strange  result  for  this  configuration  and  casts  some  doubt 
on  all  of  the  comparisons.  Usually  the  data  are  below  the  calculated  results. 
For  the  unsteady  cases  there  seemed  to  be  a  phase  shift  between  the  data 
and  calculated  results.  Also  plots  are  presented  to  isolate  the  effects 
of  the  nacelle.  Ir  general,  the  real  r-  in-phase  components  of  the  calculated 
results  and  experimental  data  agreed  wc>l  but  the  imaginary  parts  did  not 
agree  as  well.  Generally  speaking,  thf  comparisons  of  calculated  and 
experimental  chordwise  loading  showed  good  to  fair  agreement.  The 
functional  shapes  of  the  pressure  as  determined  experimentally  and  as 
calculated  are  in  good  to  fair  agreement  even  though  the  integrated  value 
cf  lift  may,  in  some  instances,  not  be. 

Finally,  a  comparison  is  made  between  calculated  and  experimental 
span  loading  for  the  8-58  bomber.  The  leading  edge  sweep  of  the  B- 58  is 
60°.  Because  of  this  a  leading  edge  vortex  ’s  formed.  This  caus0^  the 
experimental  data  to  lie  above  the  calculated  results.  The  Present 
Method  does  not  account  for  this  nonlinear  effect. 

4.2  Recommendations  for  Further  Work 

The  recommendations  for  future  work  fall  into  three  categories: 

1)  improved  efficiency,  2)  improved  modal  representation  and,  3)  improved 

body  representation. 
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A  large  amount  of  time  is  expended  in  calculating  the  kernel  function. 

The  kernel  is  evaluated  at  three  points  on  each  lifting  surface  element*. 

For  each  element  there  may  be  several  images,  a  symmetry  plane  and  a  ground 
effects  plane.  Efficiency  would  be  improved  if  appropriate  far  field 
formulas  are  used  for  receiving  points  that  lie  at  a  large  lateral  distance 
from  the  element.  A  second  and  more  radical  method  of  improving  the  efficiency 
is  to  cover  the  lifting  surfaces  with  trapezoidal  vortices.  The  trapezoidal 
vortex  developed  in  the  Present  Method  is  simple  but  not  very  efficient. 

A  new  formula  could  be  developed  for  the  trapezoidal  vortex  that  does  not 
require  the  evaluation  of  the  complicated  kernel.  It  world  be  simpler 
because  the  trapezoidal  vortex  possesses  no  wake  and  thus  the  expression 
contains  one  less  integration.  The  only  place  where  a  wake  is  required  is 
at  the  wing  trailing  edge. 

The  second  recommendation  pertains  to  an  improvement  of  the  aircraft 
modal  representation.  The  polynomial  approach  lends  itself  to  scientific 
investigation  where  modes  are  simple.  When  the  modes  become  complicated 
however  it  may  be  desirable  to  incorporate  other  more  practical  modal  input 
methods.  Several  other  methods  have  been  outlined  in  Section  2.2  and  will 
be  repeated  here.  One  possibility  is  to  input  a  set  of  modal  deflections 
at  various  spanwise  stations  along  the  wing  and  use  a  spline  fitting  technioue 
to  interpolate  for  intermediate  values  of  deflections  and  slopes  in  both 
the  span  and  chordwise  directions.  A  second  method  would  be  to  generate 
an  aerodynamic  influence  coefficient  matrix  [AIC]  for  such  a  set  of  structural 
deflection  points. 

The  third  recommendation  deals  with  a  series  of  improvements 
of  the  body  representation.  The  longitudinal  distribution  of  load  on  the 
body  surface  could  be  improved.  Currently  the  unsteady  pressure  equation 

*The  value  of  the  kernel  at  the  outboard  edge  of  an  element  is  the  same  as 
the  kernel  at  the  inboard  edge  of  the  element  on  the  next  strip  outboard. 

This  fact  is  taken  into  account  requiring  that  the  kernel  be  evaluated 
only  twice  per  element. 
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is  integrated  first  longitudinally  and  then  around  the  cross  section  for 
Alternate  II.  If  possible,  this  integration  procedure  should  be  reversed 
for  elliptical  cross  sections  sc  that  more  accurate  distributions  could  be 
obtained.  The  computational  effort  would  be  increased  but  it  may  be  worth 
the  effort.  The  theory  for  circular  cross  sections  has  not  been  automated. 

It  is  anticipated  that  it  will  furnish  more  accurate  results  for  bodies 
with  circular  cross  sections. 

Currently,  the  flow  field  due  to  a  body  is  generated  by  one  or  two 
singularities  placed  near  the  body  axis.  The  resulting  cross-sectional 
shaper,  may  not  represent  the  actual  body  shape  accurately.  The  flow  field 
i'Iose  to  the  body  is  highly  dependent  on  the  details  of  the  shape  of  the 
cross  section.  It  may  be  worth  while  to  investigate  the  possibility  of 
improving  upon  the  current  idealization.  A  more  refined  method  of  accounting 
for  the  end  plating  effects  of  a  body  would  be  appropriate.  Specifically,  a 
more  complicated  singularity  should  be  developed  to  account  for  an  anti¬ 
symmetric  flow  to  the  body  surface. 

The  entire  problem  of  wing-tail  interference  should  be  reconsidered. 

The  recommendations  of  Section  2.5.8  could  be  implemented  along  with  a 
method  of  accounting  for  wing  wake  convection  and  roll  up.  One  of  the 
recommendations  of  Section  2.5.8  is  to  use  two  different  average  fuselage 
diameters.  A  diameter  .appropriate  for  the  wing  is  used  for  wing-on-wing 
and  tail -on-wing  interactions.  A  diameter  appropriate  for  tne  tail  would 
oe  used  for  wing-on-tail  and  tail -on- tail  interaction. 
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Figure  4.  Total  Force  on  a  Body  of  Circular  Cross  Section  due  to  a  Pressure 
Doublet  Located  on  the  y-Axis. 


tal  Force  on  a  Body  of  Elliptic  Cross  Sectio 
a  Pressure  Doublet  Located  on  the  z-Axis. 


Figure  8.  Total  Force  on  Cross  Sections  of  Various  Shapes  due  to  a  Pressure 
Doublet  Located  on  the  Body  Axis. 
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Figure  11.  Comparison  of  the  Lift  Coefficient  Distribution  as  Calculated  by 
the  Method  of  Part  I  and  the  Present  Method  for  a  Wing-Fuselage 
Combination.  (Midwing,  Elliptic  Cross  Section) 
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Figure  12.  Effect  o'  Reducea  Frequency  on  the  Lift  Coefficient  Distrioution 
for  a  Wing -Fusel ac;e  Combination. 
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Figure  14.  Comparison  of  the  Lift  Coefficient  Distribution  as  Calculated  by 
the  Method  of  Part  I  and  the  Present  Method  for  a  Wing-Fuselage 
Combination  in  Unsteady  Flow.  (High-Wing,  Elliptic  Cross  Section) 


Figure  15.  The  End-Plating  Effect  of  a  Body  of  Elliptic  Cross  Section 
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Comparison  of  Experimental  and  Calculated  Span  Load 
Fuselage  Combination. 


EXTENSION 


V'S 


Figure  19.  Span  Load  Distribution  for  a  Wing-Fuselage  Combination  in  Unsteady 
flow,  (a)  Fuselage  Alone  in  Pitch,  (b)  Parabolic  Camber  of 
Fuselage  Alone. 
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Figure  20.  Comparison  of  the  Lift-Curve-Slope  Distribution  as  Calculated  by 
the  Present  Method  and  Various  Other  Methods  for  a  Wing-Fuselage 
Combination. 
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Figure  21.  Comparison  of  the  Fuselage  Axial  Load  Distribution  as  Calculated 
by  the  Present  Method  and  the  Method  of  Part  I . 
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Figure  22.  Idealization  of  a  Wing-Fuselage  Combination. 
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Figure  23.  Comparison  of  Experimental  and  Calculated  Lift-Curve-Slope 
Distribution  for  a  Swept  Wing. 
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Figure  25.  Comparison  of  Experimental  and  Calculated  Lift-Cunr'-Slope 
Distribution  for  a  Swept-W-.  ng/Fusel age  Combination. 
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figure  27.  Comparison  of  experimental  and  Calculated  Lift-Curve-Slope 
Distribution  for  e  Swept-Wing/Fuselage  Combination. 
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Fiuure  31  Comparison  of  Experimental  and  Calculated  Lift  Coefficient  Distri 
bution  for  a  Straight-Wing/Fuselage  Combination. 
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Figure  32.  Comparison  of  Experimental  and  Calculated  Fuselage  Axial  Load 
Distrioution.  (Fuselage  at  Zero  Incidence) 
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Figure  34.  Comparison  of  Experimental  and  Calculated  Fuselage  Axial  Load 
Distribution  for  a  Mach  Number  of  0.88. 
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Comparison  of  Experimental  and  Calculated  Span  Load  for 
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Figure  38.  Comparison  of  Experimental  and  Calculated  Lifting  Pressure 
Distributions  for  the  V.J.  101-C  Wing  Alone.  (Pitch  k  - 


Figure  39.  Comparison  of  Experimental  and  Calculated  Lifting  Pressure 

Distributions  for  the  V.O.  101 -C  Wing  Alone.  (Deflected  Flap 
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Figure  41.  Comparison  of  experimental  and  Calculated  Lifting  Pressure 
Distributions  for  the  V.J.  iOl-C  Wing  Alone.  (Pitch  and 


Figure  42.  Comparison  of  Experimental  and  Calculated  Lifting  Pressure 
Distributions  for  the  V.O.  101-C  Wing  Alone.  (Oscillating 
Aileron, and  Roll,  =  0.5) 
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Figure  43.  Comparison  of  Experimental  and  Calculated  Lifting  Pressure 
Distributions  for  the  V.J.  101 -C  Wing  Alone.  (Pitch  and 
Plunge,  k  =  1 .0) 
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Fiqure  44.  Comparison  of  Experimental  and  Calculated  Lifting  Pressure 
Distributions  for  the  V.J.  101 -C  Wing  Alone.  (Roll  and 
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Figure  45.  Comparison  of  Experimental  and  Calculated  Lifting  Pressure 
Distributions  for  the  V.O.  101 -C  Wing-Nacelle  Combination 
(Pitch,  k  =  0.0) 
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Figure  48.  Comparison 
Distributio 
in  Steady  F 
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Fiqure  49  Comparison  of  Experimental  and  Calculated  Lift  Coefficient 
Distributions  for  the  V.J.  101 -C  without  Nacelle.  (Pitch, 
Plunge,  k  =  0.5) 
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Fiuure  52. 


Comparison  of  Experimental  and  Calculated  Lift  Coefficient 
Distributions  for  the  V.J.  101-C  without  Nacelle,  (Flap, 
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igure  53.  Comparison  of  Experimental  and  Calculated  Lift  Coefficient 
Distributions  for  V.J.  101-C  with  Nacelle.  (Pitch,  Plunge, 
■>) 
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rigure  'lb.  Comparison  of  Experimental  ond  Calculated  incremental  Lift  due  to 
the  Nacel'-e  for  +he  V  j.  lOi-C.  (Pitch,  k  =  0.0) 
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t-igure  57.  Comparison  of  Experimental  and  Calculated  Incremental  Lift  due  to 
the' Nacelle  for  the  V.J.  10! -C.  -  (Pitch,  Roll,  kr  =  1.0) 
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Fi-gure  61.  Comparison  of  The  Span  Loading  on  the  B-58  Aircraft  for  Various 
Nacelle  Diameters. 


OF  GROUND  EFFECT 


Figure  62.  Comparisor  of  the  Span  Loading  on  the  B-53  Aircraft  as  Calculated 
with  and  \\ithout  Ground  Effect  ir;  Unsteady  Pitch. 
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APPENDIX  A.  BASIC  SINGULARITIES 


The  basic  differential  equation  governing  th  ;  potential  m  compressibh 
lsentropic  flow  is: 


P  0  ‘  0  *  0 

vvx  yy 


-  Kj-  <P  4  2|MW  0^ 

a  a 


where  P  “  1  -  M  fo "  subsonic  flow  and  where  0  ~  0e  tor 

use illatory  flo\>  . 


The  first  step  m  the  solution  of  [A-  1 )  is  to  define  V  a ->  lollows- 


•  i  Ulx 


Thu  equation  reduces  ( A  - 1 )  tc 


( A-2) 
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7  7 


-x2* 


when  X  uM/P  ^ 


An  elementary  source  solution,  0  ,  for  (A-3)  for  outgoing  waves  is- 


(A-3V 


(A-4) 


(A-5) 


wh<-  re 


(x  -£) 


(y  -  n )  '  (/  -  C) 
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Applying  Green's  theorem  to  the  lifting  .surfaces  using  (A-5)as  the  elci:  entaix 
solution  gives: 


( A-o) 


where  s  represents  a  suiface  over  which  A--  is  distributed  end  where 
<•  v  /<? N  represents  a  point  oscillatory  compressible  doublet  oriented 
normal  to  the  surface  (N  ).  Using  the  relation  between  and  0  .  given 
m  (A_d),  in  Ecjualion  (A-o)  gives  the  final  expression  lor  0  in  terms  of  A 0, 
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Since  i  MM(x-0  —  R) 


equat'on  (A-7)  can  he  rewritten  as- 
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I  his  is  Ine  final  expression  for  the  potential  in  terms  of  the  potential  jump 
distribution  on  the  surface  "S'  .  The  term  L  rep r<  sent  >  a  point  double1, 
oriented  normal  to  the  surface  (in  ine  N  direction). 

The  expression  for  the  linearized  pressure  is 
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I  h<s  terms  0^  and  0  each  satisfy  Equation  (A-l)  The  pressure 
i  ooffu  lent  then  satisfies  the  equation  also.  To  obtain  a  solution 

for  <*p  simply  replace0  by  in  Equation  (A -3). 
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The  minus  sign  in  front  of  AC  Ciaicates  a  difference  in  the  definition  of 

P 

A,i.e.,the  difference  between  upper  and  lower  surface  quantities.  For  0 
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(A-  1  1 ) 

Equation  (A- 10)  will  give  the  pressure  field  anywhere  in  the  fluiJ 

due  to  a  pressure  loading  distribution  on  the  surface  "S”.  Normall/  AC 

is  unknown.  An  integral  equation  for  can  be  formed  if  an  expression 

relating  potential  or  velocitv  to  AC  can  be  found 

P 

Com 'dor  the  expression  for  the  pressure  (A-9).  In  terms  of  ^  the 
pressure  is 
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Define  an  aocele  ration  potential  such  that 
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The  expression  for  C  is  then: 
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Placing  (A  12)  into  (A-10),  on  the  right-hand-side  gives; 
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Integrating  both  sides  of  this  equation  with  respect  co  x  from  -  oo  to  x 
gives: 
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where  -<  is  the  dummy  integration  variable,  and  (-on)  is  zero.  It 
can  easily  be  shown  that 
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Substituting  this  equation  in  the  expression  for  t2  -  i2  (  -<o)  and  making 
the  change  of  variable  t  -  x  gives: 
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The  boundary  conditions  on  the  body  surface  can  usually  be  given  in  terms 
of  velocity  at  the  general  point  x,  y,  z  in  the  direction  normal  to  the 
surface 

w 
U  CO 

I<  =  V  K  ^  .  $ 

The  term  K  is  the  usual  kernel  function. 


V0 

U~ 


(A- 


It  is  often  stated  that  if  the  acceleration  or  pressure  potential  is 
used  (Equation  (A- 13)  or  ( A- 14)  )  then  the  wake  need  not  be  considered 
since  the  surface,  S,  over  which  the  integration  is  performed  does  not 
include  the  wake.  The  wake,  course,  is  present.  In  fact  each  pressure 
doublet  (of  strength  AC  )  contains  us  own  wake.  The  wake  is  *n  actuality 

P  - 1  to  l£-£  ) 

pj  *  extending  dowi 


a  line  doublet  of  strength  e 


Mistream  from  the 


point  |  ,  to  infinity.  This  is  easilv  proven  using  Equation  (A-8'.  The 
term  L  is  the  potential  field  due  to  a  unit  point  doublet  oriented  in  N 
direction. 

of  strength  e  y  gives- 


direction.  The  potential  due  to  a  continuous  distribution  of  pcint  doublets 


r  co 

/ 

I"  U” 


U-  £  ) 

°  (x  -  t ,  y  -  H  ,  z-c,  M,  w)d£ 


which  upon  substitution  for  L  becomes 
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Making  the  transformation  t  =  x  -  4  gives 
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This  expression  is  just  K0  (Equation (A-l 4)).  In  summary  then 


.  x>  CO 

-iTT 
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K0 


(A-l  6) 


This  proves  that  the  point  pressure  doublet  xs  in  actuality  a  line  doublet 
of  varying  strength  (e  ^°hextending  from  the  point  £Q  to  down¬ 

stream  infinity.  See  Sketch  A-l. 


SKETCH  A-l 
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Modified  Acceleration  Potential  Approach 

One  further  type  of  singularity  that  will  be  of  use  is  the  modified 
acceleration  potential.  Equation  (A-13)  gives 


_  j-ii-x 


Q  =  <±>  o  U< 


or  for  surface  distributions 


1TJ  ^ 

A Q  =  Ad>  e  ® 


If  a  modiiied  acceleration  potential  is  introduced,  then 
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(A- 1 7 ) 


where  is  a  known  function  of  £.  Solving  for  Ad)  in  terms  of  AQ  gives: 
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Placing  this  expression  into  Equation  (A-8)  gives: 
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where 
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The  pressure  AC  may  be  found  using  AQ  in  Equation  (A-12)  as  follows: 

/  Qj  \  ■  OJ  4 

AC  -  \  AQ  e  "Uco  e  U<0 


(A- 21) 
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Axial  Distributions 


Thus  far  only  surface  distributions  have  been  considered.  In  order 
to  solve  body  problems  axial  distributions  of  singularities  must  be  used. 
Consider  the  three  types  of  surface  distributions  studied  so  far:  velocity 
potential,  pressure  and  acceleration  potential. 
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In  order  to 

convert  th 

'“se  into 

axial  distributions  the  surface  distributions 
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»  lumped 

into  a  line  at  h  =  q  ,  4  =  £  .  Thus 
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where  ©(h  -  n  >  4  -  4-J  is  delta  function  acting  at  the  body  axis  h  ,  4 
a  a  a 

with  the  units  of  iength.  Placing  these  values  into  the  surface  integrals 

above  gives: 
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As  the  frequency  goes  to  zero  p(£)  becomes  p(|). 
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.APPENDIX  B.  DISCRETIZATION  OF  SURFACE 
AND  AXIAL  INTEGRALS 

Equations  (A- 22)  and  (A-25)  maybe  discretized  as  follows: 

?  1  ^  A<|>s  ff  Lfydz  ,  surface  (B-l) 

S'l  ASs 

5  =  S  I  -  V  ’  axiS  (B_2) 

s  1  A£ 

s 

The  terms  AS  and  A|  indii  ate  small  elements  of  area  and  axial  length, 
s  s 

respectively.  In  steady  flow  the  integral  over  the  small  area  ASg  produces 
a  vortex  quadrilateral,  i.  e.  ,  a  vortex  that  lies  along  the  parameter  of  the 
element.  This  type  of  element  possesses  no  wake;  therefore,  for  a  lifting 
surface,  a  wake  must  be  added.  The  integral  over  a  small  axial  clement 
produces  a  line  doublet  the  length  of  the  element. 

Equations  (A-23)  and  (A-26)  may  be  discretized  as  follows: 

^  =  8ir  E  A^p  /K*«.«1/4)di  ,  surface  (B-3) 

&=1  s  Aq 

?  =  8^  E  t8  A?  =  «,/4.  n  =  V  4  =  ;a),  axis  (B-4) 

s-  1 

The  term  Ah  indicates  the  width  of  the  surface  element  in  the  plane  of  the 
element.  The  longitudinal  integration  is  performed  by  lumping  the  value 
of  the  integrand  at  the  1  /4-chord  point  of  the  element.  The  length  of  the 
element  is  A|  for  both  the  surface  and  axial  line  integral.  The  1 /4-chord 
line  of  the  element  may  be  swept,  therefore, 

*1/4  =  ^1/4  (V) 

In  steady  flow  the  surface  incegral  produces  a  horseshoe  vortex  whose  bound 
portion  lies  along  the  1 /4-chord  line  of  each  element.  This  integral  is  the 
basis  of  the  Vortex-Lattice  Method  in  steady  flow  and  the  Doublet-Lattice 
Method  in  unsteady  flow.  The  axial  integral  produces  a  semi-infinite  line 
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doublet  whose  origin  lies  at  the  1 /4-chord  point  of  the  element.  The  dis¬ 
cretizations  given  in  Equations(B- 1 ),  (B-2),  (B-J)  and  (B-4)  are  consistent 
with  each  other  even  though  a  simpler  technique  was  employed  for  Equations 
(B-3)  and  (B-4).  A  system  of  vortex  quadrilaterals  covering  a  surface  and 
wake  can  be  made  into  a  horseshoe  vortex  system  covering  only  the  surface 
and  vice  versa.  Similarly,  a  system  of  finite  length  axial  doublets  can  be 
made  into  a  semi  -  infinite  axial  doublet  system. 

Equations  (A-24)  and  (A-27)  may  be  discretized  as  iollows- 


4>  =  ^  AQs  LQ  ds  ’  surface  (B-5) 

ASs 

5  =  8^  E  f,  J  LQ(1=  Vl^a>  d£  '  axis  ,B'6) 

'-1 


The  integrations  over  the  small  area  and  line  segments  could  be  performed 
as  they  stanc;  however,  further  analysis  will  show  the.,  these  integrals  can 
be  built  up  of  expressions  previously  derived.  Consider  first  the  axis  inte¬ 
gral  where  the  actual  limits  and  (leading  and  trailing  edges  of  the 

element  A 4  )  have  been  introduced.  Define  the  axial  integral  as  I  where 
s  rs 


s 

Using  the  expression  for  Lq  given  in  Equation  (A-20)  gives: 


( B  -  7 ) 


s 


The  term  £  is  a  known  function  of  £. 
ment  s,  where 


Hdi 

Define  £ 
b  c 


as  the  center  of  the  ele- 


*1  +  h 


?!  $  £  $2 

s  s 


(B-8) 
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The  function  £  (£)  is  then  a  series  of  step  functions.  The  integral  I_ 

C  i'S 

may  be  split  into  two  parts: 


The  K(t)  terms  poss*  ss  wakes  but  the  integral  I  must  not  possess  a  wake. 
The  wake  of  the  first  term  in  Equation  (B-ll)  is  cancelled  by  the  wake  of 
the  second  term.  This  equation  represents  a  short  segment  of  doublet 
strength,  which  /aries  like  e  Ua>  ,  that  lies  between  and 

(in  Sketch  2.  6-1  the  double  arrows  indicate  that  two  equal  but  opposite 
wake  strengths  exist  at  one  point.  )  This  expression  can  also  be  used  in 
the  double  integral  of  Equation  (B-5). 


surface  (B-12; 


axis  (B-13) 
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The  surface  integral,  given  in  Equation  (B-12)  ty  the  integral  of  I  , 
represents  an  unsteady  quadrilateral  vortex  composed  of  two  unsteady 
horseshoe  vortices  (see  Sketch  (B-l)).  If  Equation  (B - 1 1 )  is  placed 
into  Equations  (B-12)  and  (B-l  3)  and  it  is  assumed  that  is  not  a 
function  of  h  ,  then 


(B-  14) 


The  first  integral  of  Equation  (B-14)  is  an  unstea  'y  horseshoe  vortex  lying 
along  the  une  £  =  while  the  second  integral  is  an  unsteady  horseshoe 
vortex  lying  along  the  line  £  =  ^  (compar  Equations  (A-40)  and  (A-41)  with 
(B-3)  and  (B-l)).  If  the  element  length  A£  Joes  var)  with  q  it  is  an  easy 
matter  to  place  it  within  the  integrals  and  h indie  it  in  the  same  way  as  K-j, 
itself  is  handled. 


W 


APPENDIX  C.  IMAGE  POINTS 


Circular  Bodies 

The  flov.  about  a  circle  in  the  presence  of  a  singularity  may  be  solved 
using  Thompson's  Circle  Theorem  (see  Milne-Thompson, Reference  (43). 
This  theorem  states  that  if  f(Z  ),  (2  =  y  +  iz)  ,  is  the  onset  flow  complex 

potential  of  the  singularity  then  the  total  flow  potential  is  given  by  F(Z  )  as 


F  (Z)  =  fvZ)  +  f(a2/Z) 


(C-J  ) 


if  the  circle  is  located  at  the  origin  of  coordinates.  The  term  f(a  / Z)  is 

the  flow  necessary  to  render  the  circle  a  streamline  in  the  presence  of 

-  2 

the  onset  flow  f(Z).  The  term  f(a  /Z,  )  is  the  potential  due  to  the  image 
singularity.  If  f(Z  )  is  the  potential  due  to  a  vortex  located  at  C  then 


f  (Z) 


=  fz  In  (Z  -U 

c.  7T 


f  (a2/Z) 


I  2 

In  (Z  -  al  ) 


'(-z'c)i 


The  first  log  term  is  the  usual  image  and  the  second  term  is  a  vortex 
located  at  the  center  of  the  circle.  This  vortex  is  expendable  since  it 
does  not  affect  the  flow  normal  to  the  circle  boundary.  The  result  is: 


Pofentia!  f  <a~/Zl> 


~7  ,  2 


f,  m  (z  -  Z) 
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The  image  vortex  is  of  opposite  sign  of  the  exterior  vortex  and  the 

2  - 

location  of  the  image  vortex  is  at  +a  /  £ 


Image  Strength 

— 

-  r 

2 

(C-2) 

Image  Location 

= 

2  .- 

a  /;  - 

‘V  4 

0^ 

p  =  Ul 


The  sketch  shows  the  image  location 

7.,  k 


Applying  the  Circle  theorem  for  the  point  source  singularity  gives 
for  the  image  potential: 

2 

f(a2/z)  =  -t  (‘"<z  ~r 1  - inZ  -  iirl 

The  very  last  tog  term  is  a  constant  and  may  be  ignored.  The  middle  log 
term  may  not  be  ignored  s4nce  it  affects  the  flow  about  the  circle.  It 
represents  a  sink  at  the  origin.  The  final  solution  for  a  source  requires 
two  images;  one  at  the  usual  image  point  and  one  at  the  origin. 
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Image  Strength 


Image  Location 


-ii; 


first  image 


(C-3) 


Image  Strength 


Image  Location 


second  image 


The  sink  at  the  origin  simply  keeps  the  circle  from  producing  a  net  out¬ 
ward  flux  of  fluid. 

The  doublet  potential  is 


f  (Z)  = 


Ji  l 

2tt  Z  -  S 


The  image  potential  is 


f  (aZ/Z) 


Zr  (  2  - 

\\ 


(C-4) 


If  we  wish  this  potential  to  go  to  zero  at  infinity  then  the  ‘"illowing  constant 
must  be  added: 


Adding  C  to  f  (a  /Z)  gives: 


f  (aZ/Z) 


(-ji)  I  1 

2rr  ;2  \  Z  -  a 2/£ 


Thug,  for  doublets 

Image  Strength 


Image  Location 


2  ,-2 

-pa/; 

(a2/p2)  C 


(C-S) 
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Elliptical  Bodies 


For  any  body  other  than  a  circle  the  image  approach  does  not 
furnish  an  exact  solution.  It  may  however  furnish  usable  approximations 
if  suitably  modified.  Borland,  Reference  (35),  suggests  a  generalization 
of  the  image  approach,  for  vortices,  through  a  transformation.  Specifically, 
the  circle-el.’ipse  transformation  is  applied  to  the  exterior  vortex  and  image 
vortex  positions.  The  position  of  the  image  vortex.  r.  j,  with  reference  to 
the  exterior  vortex  in  the  ellipse  or  physical  plane  is 


r4  +  i-  K2  (  i  +  J  r}  -  4K2  )2 

j  r2  (  l  +Jl2  -  4K2) 
a  +  b 

where  r  =  ~ 


(C-6) 


where  a  and  b  are  the  semi-major  and  semi-minor  axes,  respectively. 


This  formula  possesses  a  restriction  on  the  location  of  the  exterior 
vortex.  If  the  exterior  vortex  point  C  lies  outside  the  ellipse 


(C-7) 


then  the  image  point  will  fall  on  the  wrong  Riemann  sheet.  Thus  vortex 
images  may  be  considered  only  for  vortices  lying  within  the  ellipse 
given  in  (C-7).  If  K  — 0,  i.  a.  ,  a  circle, the  ellipse  of  (C-7)  becomes 
infinitely  large  and  thus  all  vortex  points  may  be  considered. 

A  generalization  of  this  approach  for  doublets  must  include  a 
variation  in  image  doublet  strength  in  addition  to  the  variation  in  location. 
If  a  doublet  is  constructed  of  two  vortices  of  equal  and  opposite  strength 


located  at  £  and  £  +  At,  then  there  exists  an  image  doublet  located 
at  ^  and  ^  +  At^. 


=  -  r A  (i , 


p  =  r  a  ;  i 


=  fa 


r  r  A^1  - 

0  1  =  r~TT  AU 


+  P 


a; 


If  the  limit  is  taken  then 

cH, 


I 


(C  8) 


For  a  circie  d£j/d  4 


e.  y  2 

-  -  a  /b  and  this  is  exactly  the  solution  given 

ir.  ( C  - 3 )  Differentiating  Equation  (C  — 6 )  gives: 


lii 

dt 


4  .  1  't. 

'  r  +  4  K 


(i  4  A'2  -4K2  )' 


i'2  -«2) 


[Jz 

V  b 


4K 


(C-9) 


A  second  and  more  easily  understood  approach  is  to  define  an 
eouivalent  circle  for  each  singularity  'vortex,  doublet,  source).  The 
equivalent  circle  has  the  correct  curvature  of  the  ellipse  at  a  particular 
point.  Sketch  C-2  shows  an  example. 
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SKETCH  C-2 


The  curvature  of  the  surface  of  an  ellipse  (1/a) 


a 


4. 


,2 

a  z 


ds 


k 


where 

d2y 

-a  cos  6  b2 

ds2 

,  2  .  2  o  . 

(a  sin  8  + 

,2  2„,2 
b  cos  8) 

d2z 

-  b  sin  0  a 

ds2 

,  2  .  2  . 
(a  sin  0  + 

k2  2  n  ,2 

b  cos  0  ) 

where  the  angular  measure  8  is  defined  implicitly  as 


The  result  for  a  is 


,2.2.  Z  2.. 3/2 
-  (a  sin  9  +  b  cos  0 ) 

a  =  _ _ : _  (( 

ab 

The  vector  to  the  center  of  the  circle  with  radius  a  is: 

—3*  —2* - 

c  =  r  -  n  a 

where"?  is  the  vector  to  the  point  on  the  ellipse  and  where  the  radius  of 
curvature  is  a.  The  vector  ~n  is  the  normal  vector  at  the  same  point  and 
is  given  as 


The  result  is: 


(C-13) 


This  curve  is  called  the  "evolute  of  the  ellipse". 
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APPENDIX  D.  LONGITUDINAL  INTEGRATION  OF  THE  SURFACE 
PRESSURE,  ON  A  BODY  OF  ARBITRARY 
CROSS  SECTION,  FOR  FORCE  AND  MOMENT 


The  expression  for  the  pressure  field  due  to  a  point  pressure  doublet 
is: 


C  (x ,  y,  z)  = 


ACp(s,  1,  4)6A 
4  IT 


e  i  X.  M  x  q 

¥n 


(D-l ) 


where  6A  is  the  elemental  area  over  which  the  load  AC  acts. 


coordinates  4 ,  T1 ,  4  define  the  location  of  the  pressure  point  and 
direction.  Also, 


its 


\ 


R 

r 


£  M 

V  (X-  i)2  +  P2r2 

V  (y-n)2  +  (a  -  4)2 


this  pressure  must  be  integrated  over  the  entire  body  to  give  the  lift  and 
moment. 


6F 

q 

6  M 

q 


x  dx  n  . 


(D-2a) 


(D-  2b) 
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The  term’n  .  is  not  a  function  of  x  since  a  constant  cross  section  body  has 
£ 

been  assumed.  The  vector  n  is  the  outward  normal  to  the  body  cross  section 

.  -X 

and  has  components  in  the  y  and  z  directions  only.  The  vector  ip  denotes  the 
direction  of  F  and  is  used  to  determine  the  direction  of  M. 


M  =  M  (ip  .  k)  +  N  (ip  .  j) 

where  M  and  N  are  positive  nose  up  and  nose  right  1  espectively.  ft  is 
assumed,  without  loss  of  generality,  that  the  origin  is  located  at  the  pressure 
doublet  and  that  the  moments  are  taken  about  the  origin. 


Substituting  Eq.  (D-  1)  into  (D-2a)  and  using  the  transformation  t  ~  x/pr 


gives: 


AC  6  A 


1  n  .  ip  ds 


(D-3) 


where 


i  X  (3  r  (Mt  -  n|  t  +  1  ) 
e _ dt 

•J  t2  +  1 


Making  the  substitution  Pu  =  v  t  +1  -  Mt  and  (  \  pr)P  =  _ r 

gives:  lJco 


t  = - 03 


e  k  u  M 

J=Z=—  sign  (t  -  —  )  du 
*J  -  1  P 


/,  oo  -i  k  u 

~===Z 

,  J  u^  -  1 


du  =  tt  [Y  (k)  +  i  J  (k)  ] 
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This  is  just  the  Hankel  function  of  the  second  kind  times 


i  tt  .  Thus 


i  =  i  u  H  f2)  (k) 
o 


(D-4) 


where 


It  can  be  shown  that  Eq.  (D-4)  is  proportional  to  the  two-dimensional  source 
potential  for  the  unsteady  wave  quation.  A  direct  method  of  arriving  at  (D-4) 
is  outlined  as  follows: 


/' 


i  X.  M  x  -i  X.R 

e  e 


dx 


R 


Let  =  x  -  § .  1  hen 


f 


i  X.  M{  t  +  i)  -iKR  . 

•  C  U  T 

R 


Consider  the  two-dimensional  source  potential 


0 


2-D  ~ 


-  j.  i\j\ 


R 


de 


which  is  obtained  simply  by  integrating  the  point  source  along  a  line  from 
-  co  to  4  a)  .  Again  make  the  same  transformation  t  =  x  -  £  . 
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0 


2-D 


o 

4ir 


00 

/ 


i  \M{t  +  f i  -i  '•  R 
e  c 


R 


d  r 


0. 


2-D 


4  TT 


(D-5) 


The  two-dimensional  source  solution  is  obtained  from  the  two-dimensional 
unsteady  wave  equation: 


2 

r 


+ 


r  </)  +  k2  0 


0 


(D-o) 


where  oscillatory  motion  has  been  assumed.  The  source  solution  of  D-6) 
for  outgoing  waves  is: 

0  -  -y  H  (2)  (k)  (D-7) 

4  o 

Thus  from  (D-7)  and  (D-5)  we  obtain  the  results  of  Eq.  (D-4). 

Substituting  (D_l)  back  into  (D-3)  gives  the  proper  expression  for  the 

force. 

n  .  ip  ds  (D-8) 


where 

<9H  '2)(k) 
_ o  v 

5  k 

dk 

dN~ 


w  M 

If  co 


d  r 
5N 


N 


6  F 


AC  6  A 
P 

4ri 


5H0  ;(k) 


/3H  ' 
i  7T  — L 
dk 


dk 

dN 
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placing  these  into  (D~8)  gives: 


SF 

q 


tF  ds  (D-9) 


where  i  is  the  two-dimensional  radial  unit  vector  from  the  oressure  point 
r  • 

to  a  poml  on  the  body  cross  section  (see  sketch  D- 1)  on  whic  n  the  circuit 
integral  is  to  be  taken  In  the  steady  case  (k  —  0)  the  bracketed  term  approaches 
unity. 


z 


Sketch  D-l 
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The  expression  for  the  moment  is  now  considered. 


_6M_ 

q 


AC  6  A 


~~  f 

4tt  J 


d 

dN 


I  n  i„  ds 

m  F 


(D-10) 


where 


m 


-  / 

00 


Jt2  +  1  ) 


ikpr  (Mt 
S _ dt 


J  t2  + 


1 


The  term  can  be  obtained  from  I  as  follows: 


m 


01 


i  X.  P  r  3M 


Kpr  = 


const. 


(D-ll) 


Now  I 


in  Hq(2)  (k)  and  Xpr  =  k/p  ;  thus 


e  I 

dM 


i  n  h1(2)  (k) 


d  k 
d  M 


—  =  const. 


_k 

P 


=  const. 


(D-12) 


d  k 
d  M 


aP 

dM 


const. 


P  P 


—  =  const. 


(D-13) 


Placing  (D-13)  into  (D-12)  and  the  result  into  (D-ll)  gives: 


I 

m 


Pr 

ik/p 


H 


(2)  jr 


(k) 


r  W  n  H1(2)  (k) 


(D- 14) 
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Placing  this  result  into  Eq.  (D-10)  gives 


AC  6A  /* 

P  f>  r—  [rMiT  Hj  (k)  ] - n  .  ^  ds 

4  tt  •'  Sr 


r  J~Y  [r  MttH1(2)  (k)  ] 


m  dH  (2)  (k) 

rMtr  H.  (k)  +  k  - - - 

1  Sk 


r  M  ir  [k  H  (k)  ] 


(D-l  5) 


The  final  expression  for  the  moment  is  then: 


AC  6  A 


6A  f  _  ,9.  i  .  N 

P _ f  [rM  7  k  Hq  (k)]  -I - n  .  Tf  ds  (D-l  6) 

2tt  J  r 


AC  <5A  r  {z)  _  9H 

J  M,  [h[  +  k  IV  N  »  •  V  di 


In  the  steady  case  (k  — -  0)  the  bracketed  term  approaches  zero  like  k  In  k. 


Higher  Order  Moments 


The  ntn  order  moment  is  defin  -A  as 


A  C  6  A  r  d  ^  ^ 

- P - <t - ( I  )  n  .  i„  ds 

4,  J  6 N  n  F 
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where 


I 

n 


i  X  (MX  -  R) 


R 


dx 


il  y  in 

ik/p  /  9Mn 


k/p 


-  const. 


Consider  the  second  order  moment; 


SM2/q. 


<51 

3M 


M  T  „  (2)  ,r, 

111  71  k  H1  (k) 


dh 

dMZ 


.  k  (  1  xM£ 

x  TT  — 


p3 


H12)(1c)'  163i'  H0(2)  (k) 


-1---  {  (1  +  M2)  Hx(2)(k)  -  M2k  H0(2)  (k) 


^  aiz_ 

2  or 


i  f  r  jMkrH^^OO  +  (1  +  M2)r  HQ(2)  (k)  j 


Finally  the  result  for  6M2/q  is  : 


6  M 


AC  6 A  r  i  .  N 
p  /f  r 


•  iT 


(D-17) 

i  A-  jMkr^H^'  (k)  +  r“  (HM2)!^2'  (k)}] 


,T.2u(2)-  .  2 


0 


For  steady  flow  (k  —  0)  the  Hankel  Function  i-^  H^^k)  —  ink  —  -  °°  . 

Thus  the  second  order  moment,  is  unusable.  It  seems  obvious  that  higher 
order  moments,  where  n  is  even,  will  also  give  results  that  are  unbounded. 

The  higher  order  moments,  where  n  is  odd,  will  give  finite  results  since  at 
k  -  0  the  results  are  zero.  Moments  of  higher  order  than  one  are  of  interest 

in  cete rmimng  generalized  forces;  however,  as  shown  above  this  method  can¬ 
not  furnish  them. 


ds 
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SUMMARY 


6F 

q 


6  m 

q 

f(k)  = 

p(k)  = 

k 


F 


M 


AC  6  A  r 

P _  op 

2*  / 


f  (k) 


(D-18) 


AC  6  A 
P 


2  TT 


/ 


P  (k) 


(D-l  9) 


k  H1(2)  (k) 

M  ^  k  H0(2)  (k) 

2kj.  M  r/c 
w  c 

direction  of  pressure  doublet 

direction  of  F 

i  r  =  position  vector  from  doublet  to  point  on 
body  cross  section. 

force  in  the  direction  if  i„ 

J? 


moment  normal  to  i^,  (positive  nose  up  and  nose  right). 
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APPENDIX  E.  NUMERICAL  EVALUATION  OF  THE 
FORCE  AND  MOMENT  INTEGRALS 


.The  expressions  to  evaluate  are  derived  in  Appendix  D  and  are: 

Ac  8A  /  _  f  *  N  j, 

6F/q  =  — A -  T  f(k) -  n*  i  ds 

L-n  J  r  F 


(E-l) 


6M/q 


■Ac  6  A  f  v  .  .  . 

— — ^ -  tf  P(k)  i  •  N  n  *  i  ds 

2  tt  J  r  F 


(E-2) 


The  basic  approach  will  be  to  break  up  the  cross  section  into  a  series  of 
elements  over  which  Uk)  and  P(k)  are  assumed  constant. 


6F/q 


6M/ q 


(E-3) 


(E-4) 


There  are  no  singularities  involved  in  the  integrand  of  Equation  (E-4); 
thus,  a  simple  evaluation  is  permitted. 


6M/q 


iF 


)  2 
s  =  0 


e.  • 
J 


(E-5) 


Such  is  nc-t  the  case  for  6F/q;  As  r  —  0  the  integrand  is  singular  like  1/r 
and  thus  care  must  be  taken  in  its  evaluation.  As  a  matter  of  fact  several 
procedures  were  tried  and  methods  which  ignored  the  currative  of  the  sub¬ 
elements  were  found  to  be  inaccurate  for  field  poi.its  (x,  y,  z)  lying  close 
to  the  cross  section.  The  remainder  of  this  appendix  will  be  devoted  to  the 
evaluation  of  this  singular  integral  which  will  be  termed  T. 
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It  is  assumed  that  over  each  element  the  radius  and  center  of  curvature 

a 

is  constant.  The  first  quantity  to  evaluate  is  r.  From  Sketch  (E-l)  it 
is  easy  to  see  that: 


R  +  r 


s  t 


(E-7) 


it  is  assumed  that  the  element  length  is  small  compared  with  the  radius 
of  curvature.  Thus 


-A  -IS 

€  =  n  "T— 

e  2a 


(E-8) 


2  —2 

If  terms  of  order  s“/a  and  higher  are  neglected,  then 


-A  a  A  o 

n  -  ii  +  t  •=- 

e  e  a 


(E-9) 


Noting  from  the  sketch  that 


A 

A 

A 

n 

=  k 

cos  X  - 

j  sin  X 

e 

A 

A 

A 

t 

=  j 

cos  X  + 

k  sin  X 

e 

A 

-  + 

A 

R 

=  (y 

(z  -  f, )  K 

(E-101 


A  A 

Then  solving  for  r  in  Equation  (E-7)  and  n  in  Equation  (E-9), using  (E-10), 
gives: 


n 


-(y  -  9)  +  s  cos  X  + 

s 

2a 

sin  <\{  + 

2 

(E-il) 

-  (z  -  p)  +  s  sin  X  - 

r 

w 

2a' 

cos  Xj 

S  ) 

-=—  ccs  X  -  sin  X  + 
a  ’ 

k  1 

K  1 

sin  X  +  cos  x! 
a  ' 

(E-l  2) 

2  , 

Retaining  terms  linear  in  s  / a  gives 


_ i - + 

2  2  T 
s  +  Bs  i  R 


s2  D/ a 
?  ?  2 
(s  -l-  Bs  +  R  ) 


(E-l  3) 


where 


R  sin  (  X  -  4', 


-2R  cos  (  X  -  4>p) 


(E- 14) 


and  where  the  angle  (|>R  is  defined  implicitly  as  follows 


R  cos  4>, 


y  -  n 


R  sin  4>r  =  z  - 


(E-l  5) 


The  angle  4>R  is  shown  in  Sketch  (E-l).  If 


II 

12 

j  Ny 

+  k  N 

_A 

A 

II 

j  *F 

*•  \r 

+  k  i 

(E- 1 6) 


then  the  numerator  of  the  integrand  of  Equation  (E-6)  is  (retaining  only 
terms  linear  in  s  /a) 

.a  I  2  LR  , 

r  *  N  n  •  ip  =  -=-  s  +  (J  +  )  s  +  R  K 


(E-17) 


where 


I  =.  N  ir  j3/2cos2X-y|  +  N  i-  j3/2f;in\cosX( 

y  F  I  2  I  z  FJ  ! 

Y  y  (E-18) 

+  N  ir,  {  3/2  sin2  X  -  i  !  1  N  i„  !3/2sinXcosX< 

z  P  z  |  2  )  y  Fz  I  ) 
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J  =  N  i„  -  cos  X.  sin  X  +  N  U  -  sin  X 

y  Fy  I  I  Z  Fy  I 


(E- 19) 


+ 

NZ  V 

2 

|  cos  X  sin  X 

j  + 

N  U  j 

y  1 

2  X  ) 

cos  X 

L  = 

Ny  V 

y 

|  -  cos  4>r  cos 

+  N  i_ 

2  F 

y 

|  -  sin  4>r  cos  X  | 

(E-20) 

+ 

N  U  1 

z  Ez  i 

j  -  sin  4>R  sin 

+  N,  v. 

|  | 
j  -  cos  4>r  sin  Xj 

K  =  1J  i„  cos  (j5*-,  sin  X  +  N  i„  sin  4>r>  sin  X 
Y  Fy  j  R  )  z  Fy  |  R  ) 

+  N  i_  j  -  sin  6_  cos  X  !  +  N  i_  j  ■■  cos  4>_.  cos  X! 

z  Fz  |  R  )  y  F?  I  R  | 


(E-21) 


For  convenience  let 


J  =  J  + 


(E-22) 


Then 


r  *  N  n  •  ir, 
E 


=  l  sZ  Us  +  RK 
a 


(E-23) 


The  integral  I  then  becomes 


=  f  j-i—  s2  +  Js  t  R  K 

•U* 

L  SZ  D/a _  )  , 

T  ^  \  d  S 

(R  +  Bs  +  R  )  j 


)  2  2 
l  s  +  Bs  +  R 


(E-24) 


Expansion  for  R/s  >1.0 


When  R/s  is  larger  than  1.  0  an  expansion  for  large  R/s  is  appropriate. 

1  /  s  \  ^ 

Retaining  terms  of  order  ^  I  R  I  and  lower  gives: 


=  2K(  + 


+  K 


2  /  « 


*  f  11/ 
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4  (il3  14  .  Is.) 

3  \  R /  [  a  \  R/ 
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(t  - 


1  + 


K 
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B  y 
R/  J 


R 


+  it  j  4 

a  !  r2 


1  + 


.4 

R*  I 

D  ) 


2D  ,  D 

a  +  .2 


|-(*)  i  ♦  1 

(■-«)•  4!  •'!*(■-?) 


all 


R ' 


(E-25) 
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Specifically, Equation  (E-25)  is  used  when 


(11/ e)  >1.5 


Inner  or  Exact  Solution 

For  R/e  <  1.  5  an  inner  solution  is  required.  The  integral  E-24  can 
be  performed  analytically.  A  complicated  expression  results: 


I  i  is 
I  a 


D  K'K  2 


2  2 
'  -  2R 


) 


JB  x  D  ]  b(6R2  -  F,j 

2  a  2i4R  -  B  ) 


a  (4R2  -  B2)  I  J7rZ  -  B2 
j  -IB  +  J  +  D  7  )  ln  / 

a2  2  a  2  l  y 


tan 


V4R2 


-  B^ 


R 


?  2 

e  +  R  +  eB 
2  2 

e  +  R  -  eB 


I  2e 


D 


-  I  [  (3R2  -  B2)  Be  +  (2R2  -  p2)  R2 


([ 


(4R2  -  B2)  (e2  +  Be  +  R2) 


(E-26) 
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f  -(3R2  -  B2)  Be  -f  (2R2  -  B2)  R2 
i  (4R2  -  B2)  {e2  -  Be  +  R2) 


DRK 

a 


j  r  (B2  -  2R2)  e  +  BR2 _ | 

i  .  (4R^  -  B2)  (e2  +  Ba  +  R2)  J 


?  ?  ?  \ 
-  (B  -  2R  )  e  +  BR  I 

2  2  2  2 
(4R  -  B  )  (e  -  Be  +  R  )  ) 


(E-26) 

'  continucrl 


Simplification  of  this  expression  is  desirable  to  reduce  the  computing 
effort.  Equation  (E-26)  holds  only  when  the  field  point  is  close  to  the 
cross  section.  The  elements  can  be  arranged  so  that  the  element  closest 
to  the  field  point  is  centered  with  respect  to  the  field  point.  Specifically, 
the  center  of  the  closest  element  lies  on  a  "normal"  coordinate  line  pass¬ 
ing  through  the  field  point.  Sketch  (E-2)  shows  an  example  for  an  ellipse. 

z 


y 
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An  expression  for  this  normal  coordinate  will  be  developed  later  for  an 
elliptic  cross  section.  Under  these  circumstances  the  difference  in  angles 
X  -  4>r  is  close  to  90°.  Since  B  is  proportional  to  cos  (X  -  <j>^)  it  is  small. 
An  expansion  for  small  B  may  be  made  retaining  only  linear  terms.  The 
result  is 
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J  B 
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+  R  K 


3  D  J  B  R  K  D  ) 

4  a  +  2a  j 


tan 


-l/  2Re 


,R2- 
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2  2 
e  +  R 


e  B  J 


.  D  J  B  e  (3e2+  R2) 
'*■  2  2 
2a  (e  +  RC) 


1 
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D  K  R  e 
a 


i  n  n  \ 
t  / 


for  R/ e  <1.5 


Special  Considerations  for  an  Elliptic  Cross  Section 

The  location  of  the  element  closest  to  the  field  point  must  be  chosen 
such  that  its  center  lies  along  the  normal  coordinate  curve  (which  is  a  hyper¬ 
bola)  from  the  field  point  to  the  ellipse.  See  Sketch  (E-2).  The  first  thing 
that  must  be  done  is  to  find  the  elliptic  parameter  0^  which  is  constant  along 
this  line.  Once  0^  is  kuuwu  the  un  the  eiiipse  wnere  tne  nyperooia 

mteiscCvS  it  can  be  computed.  The  equation  for  the  hyperbola  is: 

- Z_ -  _  — —  _  a2  _  b2  (E-28) 

cos  sin  0^ 

where  y,  z  are  the  coordinates  of  the  field  point  and  a  and  b  are  the  semi¬ 
width  and  semi-height  of  the  ellipse,  respectively.  Equation  (E-28)  may  be 
2 

solved  for  sin  0^ 
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-  A  + 


(E -29) 


.  2 
sin 


Ja*  +  4(a2  -  b2)  z2 
2  (a2  -  b2) 


o  7  7  7  L 

A  =  (y^  +  z  +  b  -  a  ) 


when 


2  u2 
a  =  b 


.  2,  2,2  2 

sm  o  =  z  /  (z  +  y  ) 


(E-30) 


The  quadrants  are  assigned  as  follows 


sin  0j  = 

\  sin2  0^ 

jign  z 

COS  0j  = 

/  2 

N  1  -  COS  0J 

sign  y 

(E  -  3 1 ) 


The  equation  of  the  ellipse  is  given  again  as 


a  cos  0 


z  =  b  sin  0 


(E -32) 


The  location  of  the  center  of  the  first  element  is  then  obtained  by  placing 
the  value  of  0^  in  place  of  0  in  Equation  (E-32). 

The  first  element  has  been  located.  The  other  elements  are  found  by 
dividing  the  0  coordinate  equally  from  0  =  0  360°  in  N  equal  parts. 

The  half  width  of  an  element  centered  at  0  is 


-  H  1  / 

N  a 


2  2  2  2 
a  sin  0  +  b  ccs  0 


(E-  33) 
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The  radius  of  curvature  of  an  element  centered  at  0  is: 

.2  .  2.  ,  ,2  2 
(a  sin  0  +  b  cos  0) 
a  =  j - - - 

a  b 

The  slope  of  an  element  centered  at  G  is 
X.  -  tan  *  cot  oj 

All  quantities  arc  now  known  for  the  evaluation  of  GF/q  and  6M/q  for 
elliptic  cross  section. 


(E-34) 


(E-35) 


i. 
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APPENDIX  F.  EVALUATION  OF  THE  STEADY  LIFT  FOR 
CIRCULAR  CROSS  SECTIONS 


The  most  important  contribution  to  the  lift  is  the  one  due  to  the  steady 
part.  Setting  f(ko)  =  1  in  Equation  1,  Appendix  E,  gives: 


:  6A  /  (i  •  N  )  ^  ^ 

i  f  _L - („  .  iF)  d< 

/  rt  **  r 


(F-l) 


where  AC  acts  m  the  direction  of  N  and  6F'  in  the  direction  of  i_,. 

P  P 

The  evaluation  of  this  integral  for  various  cross  sections  wi’.l  be  the  sub¬ 
ject  of  this  appendix. 


The  Circular  Cross  Section 


For  a  circle  Eq.  (1)  becomes 


6F(S)  .  ACp6A 


where  I 


'2t7  r  .  n  n  .  i„  a  d  0 
F 


(F-2) 


and  where  r  =  i  r. 
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r  .  N  n  .  i. 


ip  •  1 G  ("ra  sinO 


i_  .  is  ap  cos  0  sin  0  + 
F  w  R 


i,^  .  in  aii  -  cos  0  sin  0  + 

f  a  n  e 


(F-3 ) 


iF  •  i©  a  ^  0  sin  6 


ip  -  iR  P  y  (apR  cos  e  -  "a  cos  0) 
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where  a  rotated  coordinate  system  is  used 


(F-4) 
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z 


a 


r 

a 


Here(ya,  z&)  are  the  coordinates  from  the  bedy  axis  to  the  pressure  doublet. 

2 

Also  needed  in  Eq.  (2)  is  r  . 


2 

r 


2^2, 
r  +  a  2ar  cos  0 
a  a 


(F-5) 


Placing  (F-3)  and  (F-5)  into  (F-2)  and  deleting  integrals  that  vanish  gives: 


i„  .  is  a  p_  I  +  i  .  ip  a  p  (-  r  I  +  al  )  (F-6) 

o  F  w  0  s  F  K  R  a  c  cc  ' 


where 
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■  iTi 


sin^  6  d  6 
A  -  B  cos  0 

cos  0  d  0 
A  -  B  cos  0 


cc 


and  A 
B 


.  2ir 


r  +  a 
a 


2  a  r 


cos  0  d  0 
A  -  B  cos  0 
2 


The  technique  of  splitting  integrals  into  odd  and  even  parts  will  be  used 
The  technique  will  be  illustrated  for  I  only. 
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2tt  2 
sin  9  d8 

A-B  cos  0 


sin  0  d  0 
A-B  cos  0 


sin2  0  d  0 
A+B  cos  8 


sin  0  d0 
A-B  cos  0 


TT  2 

sin  0  d0 
A+B  cos  8 


Combining  the  first  two  and  last  two  integrals  gives: 


11  sin  0  d  8 


' u  cos  0  sin2  0  d0 


2a  /  2  2  2q  +  2B  /  2  2  2 

A  A  -B  cos  0  /  A-B  cos  0 

o  *  o 


The  second  of  these  integrals  is  the  "odd1'  contribution  and  is  zero. 


(F-7) 


where  I 


/•  sin  0 

-fr/2  | 

1  2  2  1 

A2-B2coS29  d8 
o 

4a2  r  2  A  1 
a 

i  ra  -a  1 

2  2 
'(ra  +a  M 


2  2 
a  r  6 
a 


( F-  8 ) 


where  Q  =  - 


(  I  2 


-(ra2  +  a2  )  j 


In  a  similar  manner 


I  =  2BI 
c  n 


1  =  2Ai 

c  c.  n 


(F-9a) 

(F-9b) 
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ttA  2 

fit*  ms 


where  I 


r  ' £  cos  G  d  G  it  Q 

n  =  2  ~Z  2  27  =  ,  2  21  2  2~ 

/_  A  -B  cos  0  2a  r  r  -  a 

JO  a  1  a 


(F-9c) 


The  term  (-r  I  +  al  )  is  to  be  formed  for  Equation  (6). 

cl  C  O  C 


•r  1  +  al 

a  c  cc 


-  2a  fr  2  -  a2)  I 
a  n 


(F-  10) 


Placing  (F- 9c)  into  (F-10)  and  the  result  into  (F-6)  along  with  the  result 


for  Ig  gives: 


-  -  /  ^  Q 

V  ■  ’e  z 

\  ra 


+  lF  '  lR 


2  J  siSn  (ra2-aZ  > 

1  '  (F-ll) 


where  Q  = 


2  |  I  a 


2  2  I  .  2  ,  2  .  I 

r  -a  -  ( r  4  +  a  ( 


2  2  2 

r,>  a  (  doublet  outsid°  body)  Q  =  a  ,  sign  (r  -  a  )  =  + 

cl 

2  2  2 

r  <  a  (doublet  inside  body)  Q  =  r  ,  sign  (r  -  a  )  ~  - 

a  & 


Equation  (F-ll)  can  be  written  as  two  expressions  depending  on  whether 
the  doublet  lies  internal  or  external  to  the  body.  Placing  these  two 
expressions  into  ( F  —  1 )  gives: 


Doublet  Inside  Body 


AC  6  A 


)  V  - 


(F  —  1 2) 
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Doublet  Outside  Body 


'AC  6A 


(F- 13) 


where 


^0  l0 


•Vr 


Nj  =  P0  '0  -  pR  iR 


Equations  (F-12)  and  F-l3)  describe  a  very  peculiar  behavior.  Take  for 
instance  the  case  where  these  equations  have  a  Common  point  -  a. 


6F(s)/q 


6FlS'/q 


AC  6  A 
P 
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AC  6  A 
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iF  .  N 


•  v  •  N1 


AC  dA 
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AC  6A 

~ f—  {(iF  •  de)Pg  MiF  -iR 


fF-1  ’) 


r  -  a  -t  f 
a 


Notice  that  the  contribution  of  the  component  (one  oriented  normal  to  the 
radius  vector  from  the  circle  center,  see  sketch)  is  continuous  when  the 
doublet  passes  through  the  body  surface.  However, ruch  is  not  the  case 
for  the  component  (oriented  parallel  to  the  radius  vector;. 

When  the  doublet  passes  inside  the  body  the  forc°  on  the  body  d.ie  to  fr 

changes  sign.  ,Vhen  the  doublet  is  anywhere  inside  the  body,  the  force 

is  equal  to  ^^p^A  and  is  in  the  direction  of  the  doublet,  N  . 

2  AC  oA 

When  the  doublet  is  outside  of  the  body  the  force  is  ^qual  to  — ^ —  /  (r  /a)* 

— N 

and  is  in  the  direction  of  Nj. 
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APPENDIX  G.  CIRC UMF ERENT IAL  INTEGRATION  OF 
THE  PRESSURE  FOR  BODIES  OF  CIRCULAR 
CROSS  SECTION  IN  BOTH  STEADY 
AND  OSCILLATORY  FLOW 


Steady  Flow 

The  expression  for  the  pressure  field  due  to  a  point  pressure  doublet 
_n  steady  flow  is 


Cp  (x,  y.  z) 


ACp  5  A 
4ir 


(G  - 1 ) 


where 

R2  =  (x  -  +  £2r2 

r  2  =  (y  -  ?7)2  +  (z  -  £)2 


Lawrence  and  Flax'7  have  integrated  this  equation,  in  an  approximate 
manner,  around  a  circular  cro^s  section.  Essentially,  they  obtained  only 
the  force  in  the  z-direction  (see  Sketch  G-l)  due  to  a  pressure  doublet 
oriented  in  the  same  direction.  Their  results  will  be  rederived  and 
extended  to  obtain  the  force  in  the  y-direction  due  to  a  pressure  doublet 
oriented  in  the  same  direction.  Also,  it  will  be  shown  that  the  force  in 
the  z-direction  due  to  a  doublet  oriented  in  the  y-direction  (and  vice  versa) 
is  zero. 


For  simplicity  of  notation  the  origin  of  coordinates  is  fixed  at  the 
circle  center  and  the  pressure  doublet  is  assumed  to  lie  on  the  y-axis. 
This  results  in  no  loss  of  generality  since  the  y-  and  z-axis  may  be 
thought  of  as  simply  a  coordinate  system  that  has  been  translated  and 
rotated  so  that  these  assumptions  are  satisfied. 


Equation  (G-l)  may  be  rewritten  as:  * 


Cp  (:  ,  y,  z)  =  ~ j"-  -jrr  {(a  -  O  p0  (y  -  V)  f*r  j 


(0-2) 


whei  o 


N 


ie  Ne  ■*  ir  N; 


=:  u  and  u  are  used  to  be  consistent  with  Appendix  F. 
'  0  r 
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Introducing  the  coordinates  of  the  circle,  y  =  aQ  cos  6  ,  z  =  a  sin#,  noting 
that  £  -  0  and  T)  =  r  and  expanding  in  terms  of  a  /R  gives: 

cl  0  3- 


„  ,  ,  4Cp  iA  0  ‘  (.  3  a  r.  i32  cost?  )  .. 

Cp  (x,  y.  a)  =  —fa—  -gy  .1 - -  +  .  .  .  X 

&  *  3  )  /. 


(C.-3) 


a  sin#  ix  +  (a  cos  $  -  r  )  ix 
o  ‘  0  o  a'  1  r 


where 


and 


R2  =  R^+  yS2a 2  -  2a  r  /?2cos  # 

a  ^  o  o  a 

R2  =  (x-£)2  +  £2r2 


Equation  (G-3)  may  be  integrated  circumferentially  to  give  the  total  force 
per  unit  length. 


where 

R2  =  (x-£)2+£2a2 

Notice  that  there  are  no  c^oss  terms;  i.  e.  ,  z-force  due  to  y-doublet  or 
vice  versa.  Equations  (G-4)  and  (G-5)  cover  the  entire  range  of  r  . 

One  way  to  measure  the  accuracy  of  these  equations  is  to  integrate  tf  :m 
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longitudinally  from  -®  to  -1®,  holding  constant,  and  compare  the  results 
with  the  exact  values  obtained  in  Appendix  F. 

Integration  of  Equation  {G-4)  gives: 

r 

F/q  -  j  3  (F/q  )/?xdx 


F/q  =  AC£A4  |,,0  k  +  prT)  .  ra<aQ  (G-7) 

These  results  correspond  precisely  to  the  exact  values  even  though  approxi¬ 
mate  methods  have  been  employed  in  the  derivation.  It  can  be  shown  that 
even  though  the  total  integrated  values  are  exactly  correct  the  distributions 
are  only  approximately  valid. 


ilCp  6 A  ao 

- 2 -  72 


(G-6) 


'  ra  >  ao 


Oscillatory  Flow 

The  results  for  steady  flow  are  very  encouraging.  The  exact  total 
integrated  lift  is  obtained  from  an  aporoximate  result.  A  similar  line  of 
development  may  be  made  for  oscillatory  flow. 


The  expression  for  the  pressure  reld  due  to  a  point  pressure  doublet 
in  oscillatory  flow  is: 


Op  (x,  y,  z) 
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where 


,  cuM 

A  =  /S2Uco 


Taking  the  N-derivative  and  substituting  £=0,  r]  =  r&  and,  z  =  aQ  sin# 

and  y  =  a  cos  6  gives 
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^  ,  x  4Cp  6A  _  i  X(M  (x  -  £)  -  r) 

Cp  (x,  y,  z)  =  — f- e  ' 


(G-9) 
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X  jlp  +  R^j  0  |  ‘‘e  ao  sin 6  +  %(ao  cos«  -  ra> 

s 

If  R  is  expanded  in  terms  of  &q!R^,  as  in  the  steady  case  and  if  e  is 

expanded  in  terms  of  A  r)  a$2<.os  6  / R  then: 
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To  check  the  accuracy  of  the  expressions  given  in  Equations  (G-10) 
and  (G-ll)  an  integration  will  be  performed  in  the  x-direction  from  -°o  to 
+»  and  the  results  compared  with  the  method  derived  in  Appendices  D  and 
E  which  are  the  exact  results. 


If  the  highest  order  terms  in  equations,  Eqs.  (G-10)  are  dropped, 


then: 
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The  integral  e.ppearing  here  is  precisely  the  integral  I  appearirg  in  Appendix 
D,  Equation(D-3). 


I  =  i7rH  ^(k) 

o 


(2)  — 

where  Hq  (k)  is  the  Hankie  function  of  the  second  knd  and  zeroth  order 


then 
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Sketch  (G-2)  shows  that  the  agreement  between  (G-i2)  and  the  results  of 
Appendices  D  and  E  is  almost  perfect.  If  the  highest  order  terms  are 
drooped  from  Equation  (G-ll)  and  the  results  integrated  numerically  the 
results  again  agree  with  the  results  of  Appendices  D  and  E.  This  result 
is  also  shown  in  Sketch  (G-2).  The  simplified  formulas  are  then: 
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-  Exact  (Integration  of  E-18) 

O  Approximate 

(Eq.  G-12,  anci  Numerical  Integration 
of  G-14) 


where 

Ra2  -  (x-*)2  +  ^2ra2 

R&2  =  (x-£)2  4  /S^2 

As  in  the  steady  case  the  integrated  values  agree  with  the  exact  values. 
However,  the  distributions  are  again  only  an  approximation. 

The  coordinate  system  used  to  derive  the  formulas  (y,  z)  for  the  forces 
is  one  that  has  been  rotated  an  angle  6^.  This  rotation  was  performed  so 
that  the  pressure  doublet  will  lie  on  the  y-axis  of  the  rotated  system.  A 
transformation  back  to  the  y,  z  system  is  necessary  to  complete  the  deriva¬ 
tion. 
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The  new  variables,  a(F  ^/q)/ax,  a(F  ^yVq)/ax  and  d[F  ^yVq)/ax  are 

z  y  z 

defined  implicitly  by  comparing  Equation  (G-17a)  with  (G-17b)  and  by 
comparing  Equation  (G-18a)  with  (G-I8b). 


APPENDIX  H.  INTEGRATION  OP  d(F/q)/<3x 
OVER  A  BODY  ELEMENT 


The  circumferential  integration  of  the  pressure  for  a  circular  cross 
section  was  performed  in  Appendix  G.  In  the  present  method  the  force 
distribution  is  given,  not  at  every  value  of  x,  but  at  a  set  of  points  which 
correspond  to  the  midpoints  of  a  set  of  body  elements.  The  value  of  the 

— Jk 

circumferential  integration,  d(F/q)/^x,  must  be  aver  ged  over  these  elements 
If  it  is  not  then  errors  could  be  introduced  due  to  local  variations  over  the 
element.  As  an  example,  Sketch  (H-l)  shows  what  might  happen  if  averag¬ 
ing  is  not  performed. 


Typical 
Body  Element 


SKETCH  H-l 
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The  expressioro  to  be  integrated  may  be  written  in  the  following  form. 
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The  terms  x^  and  Ax  are  the  element  midpoint  and  length,,  respectively. 
It  x.  +  t  is  substituted  for  x  in  the  expression  for  R  and  the  result 
expanded  for  small  t/R  the  result  is 
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>  r-.A  /  Ra  . 
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Placing  Equation  (H-9)  into  (H-6)  and  expanding  for  small  At  gives: 
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Define  a  new  reduced  frequency  variable,  k. 
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If  k  is  introduced  into  Equations  (H-1I)  and  (H-12)  then: 
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Approximate  formulas,  yielding  less  than  2%  error,  are  used  when 

Ax/R  S  1/2.  Let  5  -  Ax/R  and  T  =  (x.-^)/R 
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6  £  1/7  further  approximations  may  be  made: 
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The  formulas  given  in  Equations  (H-14)  and  (H-15)  are  valid  for  r&  >  a^. 

For  r  <  a  simply  replace  R  with  R  and  set  A  =  0. 
a  0  aA  aA 


(*.-«)2+  P2a2  1 


r  <  a 
a  o 
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